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Abstract 

In this dissertation, written under the supervision of Prof. Albert Fannjiang, 
we study statistical and ergodic properties of randomly perturbed (noisy) classical 
and quantum dynamical systems. We concentrate on the discrete time dynamics 
generated by Lebesgue measure preserving maps defined on rf-dimensional torus. We 
introduce the notion of the dissipation time which enables us to test how the system 
responds to the noise and in particular to measure the speed at which an initially 
closed, conservative system converges to the equilibrium when subjected to noisy 
interactions with its environment. 

We study the asymptotics of the dissipation time in the limit of vanishing noises 
and prove that it provides a robust criterion of the chaoticity of the underlying conser- 
vative system. The results formalize in a rigorous and quantitative way the idea that 
the dissipation is fast for chaotic systems and slow for regular ones. In the classical 
setting, we show that chaotic systems, e.g., Anosov diffeomorphisms possess logarith- 
mic dissipation time while for non-chaotic maps the corresponding asymptotics is of 
a power-law type. In case of diagonalizable ergodic toral automorphisms we com- 
pute the exact value of the dissipation rate constant and show that it is equal to the 
reciprocal of the minimal dimensionally averaged KS entropy among all irreducible 
components of the rational block diagonal decomposition of the map. 

In case of quantum systems we introduce the notion of the dissipation time for 
both finite and infinite dimensional quantizations on the torus. We study the simulta- 
neous semiclassical and small noise asymptotics of the quantum dissipation time and 
relate it to the notions of the Ehrenfest time and the dynamical entropy of the quan- 
tum system. We concentrate on quantum toral symplectomorphisms (generalized cat 
maps) for which we compute the exact asymptotics of their quantum dissipation time 
and show that it coincides with a classical one in the semiclassical regime in which 
the magnitude of the Planck constant does not exceed the the size of the noise. 
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Introduction 



The main subject of this dissertation is the study of statistical and ergodic proper- 
ties of noisy dynamical systems. We investigate the problems of irreversibility and 
approach to equilibrium for randomly perturbed classical and quantum systems ex- 
hibiting various degrees of chaoticity. 

The origin of irreversibility in dynamical systems is usually modeled by small 
stochastic perturbations of the otherwise reversible evolution. These perturbations 
may be attributed to many different sources: uncontrolled interactions with the envi- 
ronment, internal stochasticity of the system or unavoidable simplifications made in 
theoretical models of real-life experiments; e.g., some internal variables neglected in 
the equations. In experimental or numerical investigations, stochasticity or noise is 
introduced respectively by finite precision of the preparation and measurement appa- 
ratus, and by rounding-off errors due to finite precision of numerical computations. 
The important common feature is that noises, intrinsic (internal stochasticity) or 
extrinsic (random influence from the environment), can, on appropriately long time 
scales, induce or emphasize effects that would be absent or difficult to discern without 
noise. 

In this work we concentrate mainly on one such effect, the effect of dissipation. 
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The term dissipation refers in our study to the loss of the energy of fluctuations of 
densities, or equivalently, observables of the system during the course of noisy evolu- 
tion. The strength of dissipation can be determined by measuring the speed at which 
an initially closed, conservative system converges to the equilibrium when subjected 
to noisy interactions with the environment. The latter task can be accomplished by 
studying an appropriate time scale on which the influence of the noise becomes no- 
ticeable, i.e., affects the dynamics on characteristic spatial scales of the whole system. 
Such time scale will be called the dissipation time and will constitute the main object 
of our study. 

Intuitively speaking, the dissipation time is a time scale on which the magnitude 
of initial density fluctuations is brought below a certain fixed threshold and hence the 
system finds itself in an intermediate state, roughly speaking, 'half-way' from its final 
equilibrium. From a physical point of view, the dissipation time captures the time 
scale on which the system, due to the action of environmental noise, achieves a certain 
fixed level of Boltzmann-Gibbs entropy (cf. Section Ti. 2. 3|) . 

The main goal of this work is to determine the relation between ergodic, and 
in particular chaotic, properties of unperturbed, conservative systems and dissipative 
properties of their noisy counterparts. The main method is to study the asymptotics 
of the dissipation time in the limit of vanishing noises. Our main task in the first 
part of the work is to characterize in a rigorous and quantitative way the rate of the 
dissipation for classical systems. The results will support and formalize an intuitive 
understanding that dissipation should be fast for chaotic dynamics and slow for a 
regular one, the difference being more and more visible as the magnitude of the noise 
decreases. The fact that we are mainly interested in the small noise limit has a direct 
physical interpretation. 

Indeed, in the experimental setting considerable effort is usually made to elim- 
inate the influence of the noise on the system by isolating it from its environment 
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(at least to some reasonable degree). It is, however, impossible to eliminate the noise 
completely. In the theoretical approach such situations are usually modeled by lim- 
iting procedures (the magnitude of the noise is positive but assumed to be arbitrary 
small) . 

The notion of the dissipation time, as described above and defined in Section 
I2.2.2[ is relatively new. It has been introduced in the context of classical, continuous- 
time systems in jlHl HH] (cf- Section 12.2.11) and later developed and extended to 
discrete-time classical and quantum systems in the following series of works jSOl EH 



One of the most important findings is that the asymptotics of the dissipation 
time provides clear and robust characteristics of chaoticity of underlying conservative 
(noiseless) systems. To explain the connection between the dissipation time and 
chaoticity we need to review briefly some basic notions from the theory of chaotic 
systems and relate them to our results. 

Chaotic behavior of classical dynamical systems has been studied with an ex- 
ceptional intensity over the period of the last fifty years. Many equivalent ways of 
defining or characterizing chaos were developed over that time. Among the most 
important and well-known approaches to chaoticity one has to mention at least the 
following 

• Algebraic approach: K-systems. In this approach chaoticity is characterized 
through the Kolmogorov property (see Definition 13. 1|) . which encodes in an alge- 
braic language the idea that the system, although deterministic, behaves effec- 
tively like a memoryless process. This approach allows for some generalization 
to quantum systems ,i.e., to the noncommutative algebraic setting with the clas- 
sical definition recovered as a particular (commutative) case (cf. Section 13.1.11 
and Chapter H}. 
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• Geometric approach: Hyperbolicity. Chaoticity is characterized here by uniform 
hyperbohcity ,i.e., strict positiveness of all Lyapunov exponents. The condition 
expresses the geometric picture of two nearby trajectories separating from each 
other at an exponential rate in time. Different, slightly weaker, formulations 
are also allowed in this approach, e.g., almost uniform hyperbolicity or quasi- 
hyperbolicity (with a typical example given by ergodic but not hyperbolic toral 
automorphisms, cf. [T^llU8j ). 

• Ergodic approach: Strong mixing. In this approach, fast, i.e., at least exponential 
mixing is required if the system is to be called chaotic. This characterization is 
especially useful if the dynamics is modeled on the level of densities or observables 
of the system (not directly on the phase space). In particular the property is 
equivalent to fast (at least exponential) decay of correlations and the approach 
is particularly useful in spectral analysis (cf. Sections l2 . 31 and 1^^ . 

• Entropic approach: Positiveness of KS Entropy. The system is called chaotic 
here if it has completely positive Kolmogorov-Sinai (KS) entropy. The term 
completely positive entropy |lU6j refers to the property that the entropy of any 
partition of the phase space is strictly positive. Mere positiveness of KS entropy 
do not guarantee chaoticity of the system, as it is not difficult to construct an 
example of a nonergodic map with positive entropy (an example will be given 
later in this Introduction). One of the advantages of this approach lies in the 
fact that it gives a clear information-theoretical interpretation of the notion of 
chaoticity. 

We note that in the classical setting some of these approaches are equivalent. 
For example, by Pinsker theorem |lU6j (see also |llUj ). the system has K-property 
iff it has completely positive KS entropy. However, it is worth noting here that 
some attempts to generalize both notions to quantum systems led to nonequivalent 
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counterparts (for more details see Chapter El or [IH])- Some of these approaches allow 
one to introduce different degrees of chaoticity. In the geometric approach, different 
levels of hyperbolicity can be specified. In the ergodic formulation, one may require 
a particular speed of decay of correlations within a particular class of observables 
determined, e.g., by some regularity properties. 

In this dissertation we introduce another characteristics of chaoticity, namely 
the asymptotics of the dissipation time. The difference between the above-mentioned 
approaches and the present one lies in the fact that chaoticity is tested here in an 
extrinsic way. We test how the system responds to the noise. Given this information, 
and not necessarily knowing all the details of the underlying conservative dynamics, 
one can distinguish chaotic from regular behavior. This makes the dissipation time 
a robust criterion of chaoticity (cf. remark before Proposition 13.91 in Section [3. 1.2|) . 
Moreover, since all real-life experimental systems are inevitably subject to noisy in- 
teractions with their environments the present approach is well suited for practical 
purposes. 

Another reason for considering the asymptotics of the dissipation time as a test 
of chaoticity is that the notion can be almost literally and quite successfully adapted 
from the classical to the quantum setting. Similarly as for the classical dynamics, 
the quantum dissipation time provides here a good criterion which enables us to 
distinguish chaotic from regular behavior in an appropriate semiclassical regime (cf. 
Section [5.4.2|) . The importance of this observation can be understood if one takes 
into account the fact that despite great progress made in the field of quantum chaos 
in the last two decades, there is still no agreement on what quantum chaoticity really 
means. The main problem lies in the fact that one cannot simply take any particular 
classical definition of chaoticity and apply it directly to a quantum system. Indeed, as 
mentioned above, in classical dynamics chaoticity is usually connected with the notion 
of a trajectory of a system and the arbitrary closeness of two nearby trajectories (the 



Chapter 1. Introduction 



geometric approach), or equivalently with theoretical abihty to resolve the details of 
the phase space to arbitrary small scales (the entropic approach) . For obvious reasons 
neither of these notions has any direct counterpart in quantum case. On the other 
hand, even if some quantum generalization can be constructed (e.g., via the algebraic 
approach), the way to obtain it is usually non- unique, and the same classical notion 
can have many nonequivalent quantum counterparts (for more detailed discussion 
and references see Chapter |3]). The theory of the quantum dissipation time developed 
in this dissertation can be viewed as one of the many possible ways of approaching 
the problem of chaoticity in quantum systems. The second part of this work will be 
entirely devoted to this subject. 

Now we pass to a more systematic discussion of our main results. We start 
with the classical setting considered in the first part of the dissertation. One general 
comment is appropriate here. Namely, the notion of the classical dissipation time is 
independent of any particular mathematical model of the noisy dynamical system one 
chooses to work with. However, in order to fix the attention and, more importantly, 
be able to derive concrete, rigorous results one needs to choose a certain class of 
models for which a uniform framework can be constructed and the results for different 
systems can be compared, provided that they belong to the selected class. In this 
dissertation we choose to work with discrete time systems (maps) defined on compact 
phase spaces (represented by (i-dimensional tori) with Lebesgue measure as a natural 
invariant measure for both conservative and noisy dynamics. 

As a matter of illustration and to build up some intuition we start by presenting 
some simple but representative examples of classical maps for which exact results re- 
garding the asymptotics of their dissipation time are available. The simplest examples 
(toy models) of that kind are as follows 

I. Translations on T"^, defined by Fx = x + v, represent the simplest examples 
of nonergodic, if l,vi,..,Cd are rationally dependent, or otherwise ergodic but 
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Map 


Ergodic Properties 


Dissipation Time 


I. Translations 


Not ergodic or not weakly-mixing 


r, = 6-2 + 0(l) 


II. Cat maps 


Exponentially mixing 


r, = ^ln(e-i)+0(l) 


III. Angle doubling 


Exponentially mixing 


T, = i^ln(6-i) + 0(l) 



Table 1.1: Asymptotics of dissipation time for typical maps 

not weakly- mixing maps. 

II. Cat maps on T^, i.e., elements F G SL{2,'Z) satisfying |TrF| > 2 projected 
on the torus (see POj), provide simple examples of uniformly hyperbolic, expo- 
nentially mixing, fully chaotic systems. 

III. Angle doubling map Fx = 2x mod 1 provides an example of a uniformly 
expanding, exponentially mixing, noninvertible chaotic map. 

Let e denote the strength of the noise and T^ the the operator representing 
the action of the noisy dynamics associated with the above conservative maps on the 
observables of the system. The classical dissipation time Tc ("c" stands for "classical") 
is defined as follows (for detailed definition see Section 1^. 2. 2|) . 



Tc = min{n G Z+ : ||r,"|| < e"^} 



where ||T"|| denotes the operator norm of T". 

The corresponding asymptotics of Tc as a function of positive, but vanishing 
magnitude of the noise e is summarized in Table 11.11 

Two observations emerge from the above picture. The first general observation 
is that two qualitatively different behaviors of the dissipation time can immediately 
be noticed: 

• If the dynamics is regular then Tc{e^^) diverges (as e vanishes) in a power-law 
fashion. Here we speak of slow or simple dissipation (long dissipation time). 
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Figure 1.1: Asymptotics of dissipation time 

• If the dynamics is chaotic then Tc{e^^) has logarithmic asymptotics. In this case 
we speak oi fast dissipation (short dissipation time). 

We note that when the rate of divergence of Tc as a function of e~^, with e — > 0, 
is fast then the dissipation is slow (dissipation time is long), and vice versa. 

Figure ITTT] illustrates both behaviors and explains this terminology. Indeed, the 
number of iterations required to keep ||T"|| at a constant level (here e~^), i.e., the 
dissipation time is plotted here against the inverse magnitude of the noise for typical 
regular (upper curve) and chaotic (lower curve) systems. If we fix some small amount 
of noise, say eo = 10~^, the reduction of the norm of T" to the prescribed level is 
achieved much faster (9 iterates) in a chaotic (logarithmic) regime than in a regular 
(power-law) one (36 iterates). 

The second observation is of a rather particular nature. Namely, we note that in 
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the case of these simple chaotic systems, the constant of the logarithmic asymptotics 
is reciprocally proportional to their Kolmogorov- Sinai entropy h{F). As we will see 
later this observation does not generalize in any obvious way to higher dimensions 
(cf. Theorem 13 .Tj] and to more complicated maps. 

In view of the above observations a natural question arises whether the power- 
law and the logarithmic are the only possible asymptotic behaviors of the dissipation 
time. 

To address this question, we first note that it is possible to provide quite sharp 
general upper and lower bounds for the dissipation time within a vast class of dy- 
namical systems. Indeed, using methods from the spectral theory of non-normal 
operators and in particular the notion of the pseudospectrum (see Section IT^ . and 
investigating basic geometric properties of conservative maps (hyperbolicity and local 
expansion rates studied in Section EH)), we arrive at the following results : 

• The dissipation time of an arbitrary dynamical system generated by a measure- 
preserving map is never finite (i.e. r^ -^ cxd, as e — > 0). 

• The rate of divergence of Tc is never faster than power-law in e, 

rc<e"",aG (0,2]. 

• If the map F is C^ then the divergence of Tc is never slower than logarithmic: 



^ In(e-i) < 



\n\\DF[^_ 

where HDFUoo = sup^jg-ipd ||(DF)(a;)|| denotes the highest expansion rate of F (if 
\\DF\\oo = 1 then Tc ~ e""). 

Almost all non weakly-mixing systems (a modicum of regularity is required, 
cf. Theorem I2.12|) undergo power-law (i.e., the slowest possible) dissipation: 

re~e-",ae (0,2]. 
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The question which arose from the first observation is now reduced to the prob- 
lem of estabhshing an logarithmic upper bound for the dissipation time of a largest 
possible class of systems exhibiting some chaotic properties, or equivalently deciding 
whether there exists a map for which an intermediate (i.e., contained strictly between 
power-law and logarithmic) asymptotics hold. 

The second observation suggests that for systems with logarithmic dissipation 
time, the value of the dissipation rate constant (i.e., the prefactor of the asymptotics) 
should provide valuable information about the underlying conservative dynamics. 

Chapter 01 is entirely devoted to the study of these two problems. 

As to the first problem, we developed two different methods which allowed us 
to establish logarithmic asymptotics respectively for linear (toral automorphisms - 
Section I3.1.2|) and nonlinear (C^ Anosov diffeomorphisms - Section 13. 3|) hyperbolic 
maps in arbitrary phase space dimension. 

Both methods rely eventually on quite advanced number theoretical or respec- 
tively spectral analysis and it seems that there is no 'short-cut' way to establish 
logarithmic asymptotics for any chaotic dynamical system except for simple 1- or 
2-dimensional toy models (e.g., cat maps). 

As far as abstract (i.e., not related to any particular map) results are concerned, 
we derive in Section 12.61 a general connection between mixing properties (the rate of 
decay of correlations) of both conservative and noisy dynamics and the rate of the 
divergence of the dissipation time. In particular we show that within a large class 
of maps, strong (exponentially fast) mixing implies logarithmic dissipation time. On 
the other hand we also prove that methods used in the computation of the dissipation 
can be used to determine in certain cases the (precise) rate of decay of correlations 
(cf. Proposition 13.91 in Section 13. 1.2|) . 

Let us also comment here briefly on the second problem. The exact solution is 
now only available in the case of diagonalizable toral automorphisms. The result is 
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established in Theorem 13 .71 where the dissipation rate constant is proved to be equal 
to the reciprocal of the minimal dimensionally averaged KS entropy among irreducible 
components of the rational block diagonal decomposition of the map. At this point 
another question arises: why does the minimal dimensionally averaged KS-entropy 
appear in the constant instead of KS-entropy itself? 

The complete answer to this question is not known. However, the following 
considerations can shed some light on it. It is known that the knowledge of KS-entropy 
itself (e.g., its positiveness) is not sufficient to determine whether the system is chaotic 
or not. Indeed, consider the following toral automorphism in 4-dim represented in 
the block-diagonal form 

""l 



1 

1 1 



1 2 

The first block is simply the identity and the second is a hyperbolic automorphism 
(the famous Arnold's cat [ID])- The entropy of F is positive and equals the entropy 
of Arnold's cat, but the system is not even ergodic (toral automorphisms are ergodic 
iff no root of unity lies in their spectrum - see Section [3.1.11) . The minimal dimen- 
sionally averaged entropy is in this case and the system undergoes slow (power-law) 
dissipation characteristic for non-chaotic systems. 

The fact that the dissipation rate constant averages the KS entropy over the 
dimension of the irreducible block is of separate importance. We will not explore it 
fully here. Let us only mention the following simple example. Consider two matrices 
Fi G SL{di,Z) and F2 G SL{d2,1j) and assume that they have the same or almost 
the same spectra, but with different degeneracies. If it happens that di ^ ^2 then 
also hKsi.Fi) ^ hKs{.F2) while their dimensionally averaged counterparts are of the 
same order h{Fi) ~ h{F2). This reflects the natural intuition that if the strengths of 
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Lyapunov exponents of two systems are comparable then the degree of their chaoticity 
should also be comparable (i.e., independent of the dimension). 

As far as nonlinear maps in the context of the second problem are concerned, 
we derive lower and upper bounds for the dissipation rate constant for C^ Anosov 
systems f Theorem 13. 24|) but the exact value (and even its existence, not to mention 
its connection to the KS entropy) remains unknown. 

To conclude the description of the first part of the dissertation we want to men- 
tion that in Section 13.21 we collect some of the many possible generalizations and 
applications of the above described results, especially the ones concerning toral au- 
tomorphisms. In particular we investigate the possibility of defining the dissipation 
time for maps with degenerate noise kernels f Section 13. 2. 3p . and we study the rela- 
tions between the asymptotics of the dissipation time and some typical time scales 
encountered in the study of the so-called kinematic dynamo problem f Section 13. 2. 4|) . 

Now we pass to the description of the second part of the present work, devoted 
to the study of the dissipation time in the quantum mechanical setting. 

The second part begins in Chapter |3J called the Interludium as it constitutes a 
separate and almost independent part of this work. It is meant as a historical overview 
and in the same time as a quick but comprehensive introduction to the specific area of 
quantum mechanics on the torus. We describe there the two most important and most 
commonly encountered in the literature quantization schemes for toral maps. We put 
a special emphasis on careful explanation of their origins. We also discuss similarities 
and differences between these two approaches. In particular, we concentrate on the 
role which semiclassical analysis of spectral properties of quantum chaotic systems 
on the one hand, and the introduction of several nonequivalent notions of quantum 
dynamical entropy, on the other, played in the development of both quantization 
methods (and vice versa). 

In the first part of Chapter El (cf. Sections 15. II and l5.2|) both quantization meth- 
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ods are presented in a systematic and rigorous way. We develop a general framework 
(based on Weyl quantization), which unifies both approaches (i.e., each can be derived 
as a special case of the general scheme). This prepares the ground for an introduc- 
tion of the quantum noise (Section 15. Hj) and the notion of quantum dissipation time 
f Section 15. 4j) . The final part of the chapter is devoted to our results on semiclassical 
analysis of the dissipation time of canonical toral maps f Sections 15.4. II and I5.4.2J) . 

Before we discuss these results in more detail a word of caution is necessary here. 
Namely, as mentioned above and explained in the Interludium, the quantization of 
any classical system and in particular a toral canonical map can be performed in 
many different ways. We would like our results to depend as little as possible on the 
particular quantization scheme and this is why certain effort was made in this work to 
present the results in most unified way possible. Nevertheless, complete independence 
is not possible and it is necessary to describe precisely the quantization principles 
and methods adopted in a given approach before the results can be stated (this non- 
uniqueness of the setting and dependence on quantization procedures constitutes one 
of the most fundamental differences between classical and quantum descriptions). We 
sketch some of the quantization principles briefly here and refer to Chapters 0] and 
for a detailed presentation. 

Quantization on the torus is usually approached from the following two, non- 
equivalent points of view: 

• Finite dimensional approach. This method (sometimes referred to as canoni- 
cal quantization) was originally introduced in jHS], and later generalized and 
developed in a number of works [131 iSl EOl I^ITU5] . 

• Infinite dimensional approach. Usually referred to as *-algebraic noncommuta- 
tive deformation of the torus, see, e.g., P^ ITHj. 

In this work we do not distinguish between these two approaches in terms of 
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whether they are 'canonical' or 'algebraic' (cf. jlOl)- The reason is that the real dif- 
ference between these two quantizations lies in the choice of the fundamental Hilbert 
space of pure quantum states of the system whose classical counterpart has the 2d- 
dimensional torus as a phase space. After the choice is made (i.e., the space is chosen 
to be either finite or infinite dimensional), both quantizations can be studied in any, 
including in particular 'abstract' *-algebraic, framework. 

It turns out, however, that the finite dimensional approach is much more suitable 
for semiclassical analysis of the dissipation time, and most of the results of Chapter 
El are stated in this setting. For completeness and to illustrate the difference we 
also consider the infinite dimensional case. In particular we prove that in this case, 
regardless of the value of the Planck constant, quantum and classical evolutions of 
toral symplectomorphisms coincide. 

In the finite dimensional setting, the geometry of the torus (the phase space) 
and standard requirements of quantum mechanics (conjugacy between the position 
and momentum representations) restrict the space of admissible wave functions (pure 
states) to quasiperiodic Dirac delta combs. Planck's constant h = 2TTh is then re- 
stricted to reciprocals of integers N E Z^, and the resulting quantum Hilbert space of 
pure states is A^'^-dimensional (for detailed explanations see Section fA. 21 in Appendix 
IX| . The unitary and noisy quantum dynamics can be implemented either on the set 
of all quantum states, i.e., density operators (the Schrodinger picture corresponding 
to the Frobenius-Perron approach in classical case) or on the A^^'^-dimensional algebra 
of quantum observables (the Heisenberg picture - quantum counterpart of classical 
the Koopman formalism). 

For a fixed finite dimensional quantum system (i.e., fixed Planck constant) and 
vanishing noise the presence of a pure point, unitary spectrum of the quantum prop- 
agator forces the dissipation time to have the same, trivial (i.e., power-law) asymp- 
totics regardless whether the underlying conservative classical map is chaotic or not 
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(cf. Proposition 15. 14|) . To recover useful information about the dynamics one needs 
to perform simultaneously the small noise and the semiclassical limit. It is thus clear 
that the notion of the quantum dissipation time is intrinsically of a semiclassical na- 
ture. It is, however, not enough to consider just 'a semiclassical limit', since the way 
the limits are taken here matters considerably. 

The reason for this is the following: the classical dissipation time becomes 
larger and larger in the small- noise limit, but the correspondence between classical 
and quantum evolutions holds only up to a certain "breaking time", which diverges 
only in the semiclassical limit (see Chapter 0]for detailed discussion). Therefore, when 
seeking traces of chaoticity in quantum systems with classical chaotic counterparts 
one needs to consider sufficiently fast semiclassical limit. 

On the other hand, one has to avoid falling into triviality waiting on another 
extreme of the problem. Namely, if the Planck constant is sent to zero too fast w.r.t. 
the size of the noise, the dissipation time can be shorter or even much shorter than the 
"breaking time" and the quantization effects disappear too quickly to be noticeable 
in the asymptotic behavior of the system. 

The above situation is common to all approaches to quantum chaoticity in 
finite dimensional systems (see, e.g., results regarding the spectrum of noisy quantum 
propagators [2HI HUUj , or the study of the decoherence rate and quantum dynamical 
entropy 13 123 IMl EZl lini)- 

The problem one really needs to solve here is to determine the relation between, 
on the one hand, spatial scales represented by the the size of the noise and the 
magnitude of the Planck constant, and on the other hand, temporal scales, such as 
dissipation and "breaking" times, on which some traces of original classical chaoticity 
are still present in quantized systems. 

The difficulty of the problem lies in the fact that in order to solve the prob- 
lem one needs to control simultaneously the behavior of four different asymptotic 
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parameters (noise, the Planck constant, dissipation and breaking times). 

It may be of some interest to remark that the problem resembles (to some 
extent) some problems occurring in numerical simulations of chaotic or turbulent dy- 
namical systems. Any chaotic system develops, in a relatively short time, extremely 
complicated structures on smaller and smaller spatial scales. The quantization as 
well as numerical discretization inevitably imposes a finite resolution of the details 
of the phase space, the quantum 'mesh spacing' being constrained by the Heisenberg 
uncertainty principle AgAp > h. The introduction of a small amount of noise corre- 
sponds to numerical instabilities due to finite precision of any numerical computation 
(rounding errors). 

Just as numerical approximations break down on sufficiently long time scales, 
one expects a similar breakdown of the quantum-classical correspondence. In the 
quantum framework, the corresponding time scale (the above "breaking time") is 
usually referred to as the Ehrenfest time te (see Chapter 0}. For chaotic systems, 
the first signs of discrepancy may appear around te ~ X^^ln{h^^), where A is the 
largest Lyapunov exponent (this is the earliest time scale on which the chaotic dynam- 
ics might develop structures beyond the reach of quantum phase-space resolution). 
Numerous works, both theoretical and numerical, have been devoted to study this 
phenomenon. Recently some rigorous results |22J have been obtained, describing the 
breakdown of the classical-quantum correspondence on such a time scale. However, 
this breakdown effectively occurs for a very particular class of maps and values of h, 
and hence does not necessarily refiect the generic behavior (cf. Chapter HJ. 

For a given noise strength e > 0, a quantum dynamical system always resembles 
its classical counterpart if Planck's constant is small enough. This is refiected in 
Proposition I5.15( which states that for any canonical map and noise strength e, the 
quantum dissipation time converges to its classical counterpart (see also Corollary 
I5.16|) . For a general map, it is more difficult to determine precisely a regime where 



Chapter 1. Introduction 17 

both e, h —>■ 0, and such that classical and quantum dissipation times have the same 
asymptotics. We address this problem in the case of quantized toral automorphisms. 

For symplectic toral automorphisms in arbitrary dimension, we prove (Theorem 
15.171 and Corollary I5.18|) this asymptotic correspondence between the classical and 
quantum dissipation times r^ ~ r^ in the regime where 1 ^ e > Ch. In this case, one 
indeed has Tg < te, which intuitively justifies the correspondence. 

Around the "boundary" of this regime, that is for e ~ Ch, the noise strength is 
comparable with the "quantum mesh" , and the dissipation and the Ehrenfest times 
are of the same order. It is important to stress that our methods allows us to prove 
that the asymptotics of the quantum and classical dissipation times coincide in an 
exact way, that is, including the prefactor in front of the logarithmic asymptotics (the 
"dissipation rate constant"). 

In case of quantum coarse graining (cf. def. (J5.28J) ) we derive similar semiclas- 
sical result but under a little bit stronger assumption on the convergence rate for h - 
the existence of a positive P such that e^h~^ > 1. 

We also investigate the opposite situation, when the classical dissipation time 
is longer than te- In the situation where e -C h^^^, we obtain (for ergodic toral 
automorphisms) : 

2 



^g 



> C ( - j > Tc ~ ln(e ^) > I In n| ~ Te, 



In this case, three different time ranges can be distinguished. Before te, classical 
and quantum evolutions are identical (and do not dissipate). In the range te ^^t <^ 
Tc, the evolutions may differ, but neither dissipates yet. On the time scale Tc -^t <^ Tg, 
the classical system dissipates, while the quantum one remains dissipationless until 

t > Tg. 

This result shows in particular that after passing the Ehrenfest time the quan- 
tum and classical evolutions do not have to part from each other in an immediate and 
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complete fashion. For example, the difference between the systems is still not visible 
on that time scale if one restricts the observations only to dissipative properties of 
both systems. 

To finish the description of our results we want to remark that a lot of important 
questions remain still open in this area and that the whole theory is rather in the 
initial stage of its development. The first problem is to generalize the above sharp 
estimates to nonlinear quantum maps (e.g. Anosov diffeomorphisms) . This requires 
the application of different techniques (cf. ^])- In particular appropriate estimates 
are needed on the constant in Egorov theorem for such maps J2H] (the problem is 
nonexistent in the linear case due to the fact that the semiclassical approximation is 
exact). 

Also recent results on quantum dynamical entropy for finite dimensional systems 
[3 HE] seem to indicate that just like in the case of classical maps, it should be 
possible to establish a link between the dissipation rate constant and the quantum 
mechanical entropy of the system. It is still too early however to formulate any 
concrete conjectures (since the relation is not yet established for a sufficiently general 
class of classical systems) . We discuss this problem briefly in Chapter HI 

We now conclude the Introduction with a few remarks regarding the organiza- 
tion of the material in this work. The dissertation is divided into five chapters and 
two appendixes (in the list below the descriptions do not necessary match the titles) 

Chapter 1. Introduction. 

Chapter 2. Definition and general properties of dissipation time. 

Chapter 3. Dissipation time of classically chaotic systems on T*^. 

- Section 3.1 and 3.2 - Linear maps: toral automorphism and generalizations. 

- Section 3.3 and 3.4 - Nonlinear maps: Anosov systems. 
Chapter 4. Interludium. 

Chapter 5. Dissipation time of quantum systems on T*^. 



Chapter 1. Introduction 19 

- Section 5.1 - Weyl quantization on the Torus. 

- Section 5.2 - Semiclassical analysis of dissipation time. 
Appendix A. The dynamics of Cat maps. 

Appendix B. Wigner transform. 

Two chapters (1 and 4) are of expository character. The remaining chapters 
contain the main results of the work and constitute the original contribution to the 
field. Appendixes contain necessary technical background regarding two specific top- 
ics. Throughout the work attention was paid to ensure mathematical correctness 
and completeness. All results are presented with full proofs. For the convenience of 
the reader, technical proofs of secondary importance are collected at the end of each 
non-expository chapter in a separate section. 



Part I. Classical Mechanics 21 



Part I 



Classical Mechanics 



22 Chapter 2. Dissipation time - definition, properties and general results 



Chapter 2 

Dissipation time - definition, 
properties and general results 



2.1 Setup and notation 

2.1.1 Evolution operators 

Let (T'^, 'B{T'^),m) denote the rf-dimensional torus, equipped witli its a— field of Borel 
sets and the Lebesgue measure m. Let F : T^ ^ T'^ be a map on the torus preserving 
the Lebesgue measure: for any set B G S(T^) we have m{F^^{B)) = m{B). In 
general, F is not supposed to be invertible. In the following we call such a map 
'volume preserving' with implicit reference to the Lebesgue measure. 
The map F generates a discrete time dynamics on T'^, which in terms of pathwise 
description can be represented by the forward trajectory {F"(£Co), n G N} of any 
initial point (particle) Xq G T'^. However, instead of looking at the evolution of a 
single particle, one can consider the statistical description of the dynamics, that is 
the evolution of a density (more generally a measure) describing the initial statistical 
configuration of the system. 
Let Ai(T'^) denote the set of all Borel measures on T'^. For any /i G Ai{T'^) and 



Section 2.1. Setup and notation 23 

/ G CO(T'^) we write 

/"(/) = / f{x)di2{x). 

The map F induces a map F* on M. (T"^) given by 

(F»(/)=M/°i^), forall/GC°(T'^). 

This map can also be defined as follows: 

{F*fi){B) = fi{F-\B)), for all B e S(T^). 

In particular if /i = 6x0 then F*{jj) = Sf(xo) and one recovers the pathwise description. 
If /i is absolutely continuous w.r.t. m, then F*{fi) preserves this property (since the 
measure-preserving map F is nonsingular w.r.t. m, see j79t p.42]). The corresponding 
densities 9 = -j^ ^ L^(T'^) are transformed by the Frobenius- Perron or transfer 
operator Pp |T^ : 

\ (irri / dm 

If the map F is invertible, Pp is given explicitly by: 

flW*'rn 

{Pf9){x) = {g o F-'){x)-—{x) = g o F'^x). 

am 

If the map F is differentiable, and the preimage set of x is finite for all x, the Perron- 

Frobenius operator is given by 

9iy) 



y\F{y)=x I ^^^^1 

where Jriy) is the Jacobian of F at y. 

On the other hand one can consider the dual of the Frobenius- Perron operator, called 
the Koopman operator, which governs the evolution of observables / G L°°(T'^) in- 
stead of that of densities g E L^. The Koopman operator Up is defined as 

UFf = foF (2.1) 
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Due to the nonseparability of the Banach space L°°(T'^), it is often more convenient 
to consider its closure in some weaker L^ norm, which yields larger (but separable) 
spaces of observables U'iT'^). Here we will be mainly concerned with the space L^(T'^) 
and its codimension-1 subspace of zero- mean functions 

Ll{r) = {/ G L\T^) : m(/) = 0}. (2.2) 

This subspace is obviously invariant under Up and Pp, due to the assumption F*m = 
m. Throughout the whole work, || ■ || will always refer to the L^-norm (and corre- 
sponding operator norm) on Lq(T'^) (any other norm will carry an explicit subscript). 
For any measure-preserving map F, the operator Up is an isometry on L^(T^) and 
Lq(T'^). When F is invertible. Up is unitary on these spaces, and satisfies Up = 
Pp' = Pp-i. 

Although just introduced operators will mostly be considered on Lq(T'^), we will need 
from time to time to act on some more specific spaces defined usually by certain 
regularity properties of functions belonging to them. Most typically these will be 
Holder and Sobolev spaces, the definitions of which we now briefly recall and use this 
opportunity to fix the appropriate notation. 

For any m G N, we denote by C™(T'^) the space of ?72-times continuously differentiable 
functions, with the norm 

\ct\i<m 

(we use the norm |q;|i = ai + . . . + a^^ for the multiindex ck G N'^). For any s = m + r] 
with m = [s] G N, r] G (0, 1), let C''^{T'^) denote the space of C™ functions for which 
the ?72-derivatives are r^-Holder continuous; this space is equipped with the norm 

c. = ll/llc.H- E-P ^3^p^ ■ 

\a\i=m 
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The Fourier transforms of functions g G L}{W^) and / G L}{T'^) are defined as follows: 



V| G M^ 


^(0 = 






(2.3) 


Vfc G Z^ 


/(fc) = 




= (efc,/). 


(2.4) 



Above we used the Fourier modes on the torus ek{x) := e^'^*'*'^. For any s > 0, 
we denote by i/*(T'^) and if*(]R^) the Sobolev spaces of s-times weakly different iable 
L^-functions equipped with the norms || ■ Wh" defined respectively by 






2 _ 



J2{i + \k\'r\f{k)\'. 



Finally, for any of these spaces, adding the subscript will mean that we consider 
the (f/i?-invariant) subspace of functions with zero average, e.g. Cq(T'^) = {/ G 
C^iT^), m{f) = 0}. 

2.1.2 Noise operator 

To construct the noise operator we first define the noise generating density i.e. an 
arbitrary probability density function g G L-'^(]R'^) symmetric w.r.t. the origin g{x) = 
g{—x). The noise width (or noise level) will be given by a single nonnegative param- 
eter, which we call e. To each e > there corresponds the noise kernel on W^: 

g,{x) = -g (^) , 

with the convention that g^ = 6q. The noise kernel on the torus is obtained by 
periodizing g^, yielding the periodic kernel 

9e{x) = Y^ geix + n). (2.5) 

We remark that the Fourier transform of g^ is related to that of g by the identities 

g,{k) = g,{k) = g{ek). 
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The action of the noise operator G^ on any function / G Lg(T'^) is defined by the 
convolution: 

GJ = ~g,*f. 

As a convolution operator with kernel from L^, G^ is compact on Lq(T'^) (if g is 
square-integrable, G^ is Hilbert-Schmidt). The Fourier modes {efc}o^fcgz<^ form an 
orthonormal basis of eigenvectors of G^, yielding the following spectral decomposition: 

V/gL2(T'^), GJ= J2 ^(efc)(efc,/)efc. (2.6) 

This formula shows that the eigenvalue associated with Cfc is g{tk). Since g is a. 
symmetric function, this eigenvalue is real, so that G^ is a self-adjoint operator. Its 
spectral radius rsp{Gf) is therefore given by 

r,p(G,) = ||GJ = sup \g{tk)\. (2.7) 

Since the density g is positive, g attains its maximum ^(0) = 1 at the origin and 
nowhere else. Besides, because g G L^(M°'), ^ is a continuous function vanishing at 
infinity. As a result, for small enough e > 0, the supremum on the RHS of ()2.7|) is 
reached at some point ek close to the origin, and this maximum is strictly smaller 
than 1. This shows that the operator G^ is strictly contracting on Lq(T'^): 

Ve>0, ||GJ = r,p(G,) < 1. (2.8) 

In the next section we study this noise operator more precisely, starting from appro- 
priate assumptions on the noise generating density. 

2.1.3 Noise kernel estimates 

In this subsection we present some estimates regarding the noise operator. These 
estimates will later play a crucial role in the derivation of the asymptotics of the 
dissipation time. 
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We will be interested in the behavior of the system in the limit of small noise, that 
is the limit e ^ 0. It will hence be useful to introduce the following asymptotic 
notation. Given two variables a^, b^ depending on e > 0, we write 

tte < K if limsup-^ < oo, (2.9) 

a, ^ b, if lim^ = 1, (2.10) 

a, ~ b, if a, < 6, and 6, < a,. (2.11) 

In order to obtain interesting estimates on the noise operator G^, it will be necessary 
to impose some additional conditions on its generating density g, regarding e.g. its 
rate of decay at infinity, or the behavior of its Fourier transform near the origin. 
The weakest condition we are going to impose is the existence of some positive moment 
of g, by which we mean that for some a G (0, 2], 

Ma= \x\'^g{x)dx < oo (2.12) 

(we take the length \x\ = (x^ + . . .+x^)^/^ on R'^). This condition implies the following 
properties of the Fourier transform g (proved in Section 12. 8|) 

Lemma 2.1 For any a G (0,2] there exists a universal constant Ca such that, if a 
normalized density g satisfies \2.1'^) . then the following inequalities hold: 



V^ G M^ < 1 - g{i) < C„M„|ir. (2.13) 

// \2.1^J\J holds with a = 2, we have the more precise information: 

1 ~ 9{i) ~ I^P ^'^ ^^6 limit ^ ^ 0. 

In the case a < 2, we will sometimes assume a stronger property than ()2.13|) . namely 
that 

1 - g{i) ~ l^r in the limit ^ ^ 0. (2.14) 
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Note that this behavior imphes a uniform bound 1 — g{^) < C'lCT for any 7 < a and 
C independent of 7. 

Typical examples of noise kernels satisfying (|2.14j) include the Gaussian kernel and 
more general symmetric a— stable kernels |119^ p. 152] defined for a G (0,2]: 

g,U^) := J2 e-^^^-'-^^^^'e^ix), (2.15) 



where Q denotes an arbitrary positive definite quadratic form. For the values of a 
indicated, the function ge,a{x) is positive on M.'^. 

In view of Eq. ()2.7|) . the properties ()2.12|) or ()2.14|) determine the rate at which G^ 
contracts on Lq(T'^). For instance, ()2.14j) implies that in the limit e — > 0, 

l-||GJ~e". (2.16) 

The following proposition describes the effect of the noise on various types of observ- 
ables, in the limit of small noise level. The proof is given in Section 12.81 

Proposition 2.2 i) For any noise generating density g G L^(M'^) and any observable 
f G Ll{T'^), one has 

\\GJ-f\\^^0. (2.17) 

To obtain information on the speed of convergence, we need to impose constraints on 

both the noise kernel and the observable. 

a) If for some a G (0,2] the kernel g satisfies \2.1^) or \2.i4 ), then for any 7 > 



there exists a constant C > such that for any observable f G H'^iT'^), 

\\GJ-f\\<Ce<^^f\\H.^., (2.18) 

where 7 A a := min{7,a}. If f E C^(T'^), the above upper bound can be replaced by 

\\GJ - f\\ < Cei^1|V/|| < Ce^^llV/IU- (2.19) 

Using the noise operator, we are now in position to define the noisy (resp. the coarse- 
grained) dynamics generated by a measure-preserving map F. 
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2.1.4 Noisy evolution operators 

The noisy evolution through the map F is constructed by successive application of 
the Koopman operator Up and the noise operator G^. The noisy dynamics is then 
generated by taking powers of the noisy propagator 

T, = GJJp. 

In general, the operator T^ is not normal, but satisfies rsp{T^) < \\T^\\ = \\G^\\. We 
will also consider a coarse-grained dynamics defined by the application of the noise 
kernel only at the beginning and at the end of the evolution. Hence we define the 
following family of operators: 

TW = G,U^G„ neN. 

In view of the contracting properties of G^, the inequalities ||T"|| < ||Ge||", ll^e II ^ 
IIGep imply that both noisy and coarse-grained operators are strictly contracting on 



2.2 General definition of dissipation time 

Once the notation has been set up, we can pass to the precise definition of the dissipa- 
tion time. We prefer to start, however, with some remarks regarding the motivation. 
In particular we briefly recall the original, continuous-time physical setting considered 
in jlHl; where the notion had been introduced for the first time. We then generalize 
the original definition to an abstract, discrete-time setting and after discussing some 
of the most basic properties of this notion we conclude the section with another phys- 
ical interpretation, this time coming from statistical physics and expressed in terms 
of Boltzmann-Gibbs Entropy. 
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2.2.1 Motivation 

In order to gain a bit of pliysical intuition it is useful to consider the original problem 
described in jlH] , where the notion of dissipation time was introduced in the context 
of continuous-time dynamical system of a passive tracer in a fluid flow. The flow is 
prescribed in terms of given periodic, incompressible velocity filed and the path of the 
tracer is randomly perturbed by collisions with fluid particles. The latter phenomena 
being modeled by standard Brownian motion. In pathwise description the evolution 
of the tracer is given by the following Langevin-type equation 

rfx^(t) = u(x^(t))dt + y/edw{t), V ■ u(x) = 0, 

where w stands for the standard Brownian motion. In Statistical description, the 
dynamics on the densities is given by the corresponding Fokker-Planck (advection- 
diffusion) equation 

|^ = u-Vp+|Ap, V-u = 0, (2.20) 

Let us note that unperturbed skew-symmetric operator u ■ V generates an unitary 
(conservative) group f/* = e*"'^, which corresponds (for t = 1) to our conservative 
Koopman operator Up for some F. 

Now, the perturbed generator £e = u-V+|A gives rise to a semi-group of contractions 
pt _ ^tCc^ which heuristically can be thought of as continuous time-1 counterpart of 
the generator T^ = GJJp of our discrete-time noisy dynamics in case of Gaussian 
noise (we say 'heuristically', because due to noncommutativity of both terms in C^ 
there is no obvious way of associating T^ with P^ for general velocity fields u). 
In order to predict long time behavior of the tracer and in particular the influence of 
the noise on its trajectory, one is interested in determining the speed of contraction 
of the semigroup P^. In |1H] (see also more recent version in |49 ) the following time 
scale tdiss defined by equation ||P^**ss || = x/2 was suggested for consideration regarding 
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this problem and termed as dissipation time. Determining the asymptotics of trfj^s 
(in continuous-time setting) proved however to be exceedingly difficult and except for 
some very special and simple cases (e.g. cellular flow considered in ^\) the problem 
has not been solved up to date. Also there is no known example of non-trivially short 
dissipation time (i.e. logarithmic in e~^) in continuous setting. The main difficulty, 
as was already observed in JIBj , lies in the fact that to get this result one would need 
to consider fully chaotic system, while it is very difficult to construct one i.e. to write 
down a simple differential equation which would exhibit fully chaotic behavior. On 
the other hand there is no problem in constructing fully chaotic discrete-time systems 
and in fact the ergodic theory literature abounds in such examples. This situation 
provided natural motivation for the generalization of the notion of the dissipation 
time to discrete-time systems. 

2.2.2 The definition 

In its general form the classical dissipation time Tc{p) for discrete-time noisy dynamics 
Te is defined in terms of the norm || • ||p^o on the space Lq(T'^) and w.r.t. a threshold 
r/e(0,l). 

Definition 2.3 Let T^ denote discrete-time noisy dynamics. We define 

T,{p, ij) := min{n G Z+ : ||T;^||p,o < ^}, 1 < p < oo. (2.21) 

The fact that T^ acts on every Lq(T'^) as a strict contraction ensures existence and 

uniqueness of Tc{p, t]), for each t] G (0, 1). 

Figure ITT] illustrates the definition. 

We need to show that the value of the threshold r] in ()2.2H) does not affect the order 

of divergence of Tc{p, t]), as e tends to zero. 

Proposition 2.4 For any < fj,r] < 1, Tc{p,fj) ~ Tc{p,ri). 
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Figure 2.1: Definition of the dissipation time 
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Proof. Assume < fj < rj < 1. Obviously Tc{p,fj) > Tc{p,ri). On the other hand let 
k he a positive integer such that t]^ < fj. Then 

IKe IIp.O ^ 'I ^ W-'-e llp,0 ^ '/ ^ '/• 

Hence kTc{p,ri) > Tc{p,fi), which implies Tc{p,ri) ~ Tc{p,fi). ■ 
Following the argument of |lllj one can use the Riesz convexity theorem to establish 
also the asymptotic equivalence of the Tc{p), for all 1 < p < oo (to alleviate the 
notation we drop rj). 

Proposition 2.5 i) For any 1 < q,p < oo, Tc{q) ~ Tc{p). 
u) For any 1 <p < oo, Tc{p) < Tc{1) and Tc{p) < rc(oo). 

We postpone a technical proof of this proposition to Section 12.81 

In view of the above results, instead of working in general setting, one can choose some 

convenient values of p and rj and perform, without any loss of generality, all necessary 

asymptotic calculations in one notationally simplified setting. We will usually choose 

p = 2 and t] = e~^ for computational convenience. Following this choice we introduce 

the convention that Tc{p) := Tc{p,e'^) and Tc := rc(2,e^^). The obvious dependence 

of the dissipation time on e will always be implicitly assumed but rarely explicitly 

denoted. 

The corresponding dissipation time for a coarse-grained dynamics is defined in fully 

analogous way and denoted respectively by nc{p,ri), hdp) and he- In particular 

fc := min{n G N : ||tj")|| < e"^}. (2.22) 

We note that the dissipation time does not depend on whether the dynamics is ap- 
plied to densities (i.e. by the Frobenius-Perron operator) or to observables (by the 
Koopman operator). Indeed, the norm of an operator equals the norm of its adjoint 
p. 195], so that 

||tJ")|| = \\G,U^G,\\ = ||(G,f/^a)*|| = liaP^GJ, 
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and similarly for the noisy operator T^. In particular, for invertible maps the dissipa- 
tion time does not depend on the direction of time. 

As mentioned in the Introduction, we will distinguish two qualitatively different 
asymptotic behaviors of dissipation time in the limit e — > 0. We say that the op- 
erator Te (or the map F associated with it) respectively has 

I) simple or power-law dissipation time if there exists /9 > such that 



Tr ~ 



1/^^ 



II) fast or logarithmic dissipation time if 

Tc ~ ln(l/e). 

We will also speak about slow dissipation time whenever there exists some (3 > s.t. 

re > 1/e^. 

In case of logarithmic dissipation time, the dissipation rate constant Re, when it 
exists, is defined as 

i?, = lim— ^. (2.23) 

e^o ln(l/e) 

A similar terminology will be applied to the coarse-grained dissipation time fc. 

2.2.3 Physical interpretation via Boltzmann-Gibbs entropy 

In this section we briefly discuss the connection between dissipation time and Boltz- 
mann-Gibbs entropy. The results formalize an intuitive physical interpretation of the 
dissipation time outlined in the Introduction. 

First we note that on the scales exceeding Tc, the Boltzmann-Gibbs entropy ap- 
proaches its maximal equilibrium value (i.e. 0) as can be seen from the following 



Section 2.2. General definition of dissipation time 



35 



simple estimate (cf. jZHl)- Let us first restrict considerations to bounded initial 
states, i.e., / > 0, / € L°° and ||/||i = 1. Let 



ri{u) 



—ulnu, u > 
0, u = 



and let -D„ = {x E T'^ : 1 < T^f}- On one hand, we have 



v{T:^f{x))dx 



D„ 



< 



Dn 



T^ficc) 



drjiu) 
du 

l + hiM 



du 



dx 



< sup 

1<«<I|T,"/IU JD 

< (l + ln||T;/|U)||T;/-l||i 



T:f{x)-i\dx 



< (1 + ln 



l^eV-llll. 



(2.24) 



On the other hand, we have 



> / ^{T:i\x))dx > [ ^{T:f{x))dx. 

JT'i Jd„ 

In view of the inclusion relation: L°°(T"') C L^(T°') C L^(T'^), we then obtain that 
for n ^ Tc 



sup 

/>0,||/||oo<C 



7^{T:f{x))dx 



Jd 



40 



0, Vc > 0. 



For unbounded initial states, we note that, by Young's inequality. 



|Tr/||oo < IITJIU < he 



1 = \\9e 



from which we have, instead of ()2.24|) . the following estimate 



v{T:fix))dx 



D„ 



<(l + ln||(7e||oo)||r;/-l||i. 



where in view of ()2.16p 



Inll^fjoo ~ ln(l/e). 
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Therefore for sufficiently fast diverging n ^ Tc(1) such that 

ln(l/e)||T;(/-l)||i,o^O (2.25) 



one obtains 

40 



sup 

/>0,||/||i=l 



r^{T^f{x))dx 

Yd 



0. 



The condition ()2.25p typically results in a slightly longer time scale than Tc{1). 

On the other hand, we can bound the Li distance between the probability density 

function / and the Lebesgue measure by their relative entropy via Csiszar's inequality 



\f{x)-g{x)\dx<J2 I f{x)\n{f{x)/g{x))dx 

ijd y jjd 

with g{x) = 1. We see immediately that the decay rate of 



sup 

/>0,||/||i=l 



v{T^f{x))dx 
provides an estimate for Tc{1) and, consequently, for Tc{p),p G (1, oo). 

2.3 Dissipation time and spectral analysis 

In this section we investigate the connection between the dissipation time of the noisy 
propagator T^ and its pseudospectrum together with some spectral properties of Up 
and Ge- All the operators considered in this section are defined on Lq{T'^). In the 
framework of continuous-time dynamics, some connections have recently been ob- 
tained between, on one side, the pseudospectrum of the (non-selfadjoint) generator 
A, and on the other side, the norm of the evolution operator e^^ jSH]. We con- 
sider complementary, discrete-time setting, which allows for generalizations and more 
transparent proofs. We start with the definition of the pseudospectrum, and then de- 
rive general abstract lower and upper bounds for the dissipation time. In the following 
sections we will apply these results to determine the asymptotics of the dissipation 
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time under some dynamical assumptions regarding the underlying conservative maps 
(e.g. lack of weak- mixing) . 

2.3.1 Pseudospectrum 

In this short subsection we define the pseudospectrum of a bounded operator |123j 
and state some of its properties. 

Definition 2.6 Let T be a bounded linear operator on a Hilbert space Ti (we note 
T G CiTi)). For any 5 > 0, the 5 -pseudospectrum of T (denoted by csiT)) can be 
defined in the following three equivalent ways: 

(/) as{T) = {XeC:\\{X-T)-'\\>S-'}, 
{II) as{T) = {XeC:3ven, \\v\\ = l, \\{T - X)v\\ < 6}, 
(III) asiT) = {A G C : 3fi G C{n), \\B\\ < 6, A G a(T + B)}. 

We will apply these definitions to the operator T^. For brevity, the resolvent of this 
operator will be denoted by Re{X) = (A— TJ^^. We call S"" the circle {A G C : |A| = r} 
in the complex plane, and define the following pseudospectrum distance function: 

d,{r) := mi{6 > : asiT,) n 5" ^ 0}. 

From the definition (I) of the pseudospectrum, one easily shows that this distance is 



d-\r) = sup ||i?.(A)||. (2.26) 

|A|=T- 



also given by 



We have the following property (proved in Section IT^ 

Proposition 2.7 For any isometry U and noise generating function g, one has 

4(1)^^0. (2.27) 

This means that for any fixed 6 > 0, the pseudospectrum as{T,) will intersect the unit 
circle for small enough e. 



38 Chapter 2. Dissipation time - definition, properties and general results 

2.3.2 General bounds for the dissipation time 

In this section we consider both fully noisy and coarse grained dynamics. We start 
with 'non-finiteness' results. 

Proposition 2.8 For any measure-preserving map F and any noise generating func- 
tion g, both fully noisy and coarse-grained dissipation times diverge in the small-noise 
limit e — > 0. 

Proof. We skip the subscript F to alleviate the notation. We only use the fact that 
U = Up is an isometry. We start with the full noisy case and prove by induction the 
following strong convergence of operators 

From Proposition 12. 2K ). this limit holds in the case n = 1. Let us assume it holds at 
the rank n — 1. Then we write 

rj^nj ^ uT^-if + (^G, - I)UT^-^f. 

From the inductive hypothesis, T^~^f >U"'~'^f, so that the first term on the RHS 

converges to f/"/. Applying Proposition I2.2fc ) to the function f/""/, we see that the 
second term vanishes in the limit e ^ 0. From the isometry of U, we obtain that 

for any n > 0, \\T"'\\ ^1, so that r^ >oo. In coarse-grained version, similarly as 

above, we have 

lir/-t/"/ll = ||G,f/;(G,-/)/ + (a-/)f/"/|| 

< ||(G,-/)/|| + ||(G,-/)t/'^/||^0. ■ 

Now we pass to abstract spectral bounds. 

Theorem 2.9 For any isometric operator U on Lq(T"') and noise operator G^, the 
dissipation time of the noisy evolution operator T^ = GJJ satisfies the following esti- 
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mates: 



^-f^ < r, < ,, ,,,\,,,., + 1, (2.28) 

4(1) llndlGJ)! 

1 / e 



re < inf T^-TTT In ^^T • (2-29) 

rsp(TO<r<i |ln(r)| \de{r)J 

We notice that the first upper bound does not depend on U at all, but only on the 
noise. Using the estimate ()2.16|) . we obtain the following obvious corollary: 

Corollary 2.10 If the noise generating density satisfies the estimate fl2.14|) for some 
a G (0,2], then for any measure-preserving map F the noisy dissipation time is 
hounded from above as follows r^ < e~°. 

Proof of Theorem VIM 

1. Lower bound 

We use the following series expansion of the resolvent |125| p. 211] valid for any |A| > 

oo 

i?^(A) =^A-"^^T;^. (2.30) 

n=0 

Considering that VspiT^) < \\G^\\ < 1, we may take |A| = 1, and cut this sum into two 
parts: 



R^{\) = Y^ X-'^-^T^ + A-"^T;=i?e(A). 



ra=0 

Taking norms and applying the triangle inequality, we get 

Tc-l 



Pe(A)|| < llJ^A— iT;|| + |Ar-||T-||p,(A)|| 

n=0 

< re + e-ii?,(A)|| 
^||i?,(A)||(l-e-i) < r,. 

Taking the supremum over A G Si yields the lower bound. 



40 Chapter 2. Dissipation time - definition, properties and general results 

2. Upper bounds 

To get both upper bounds, we use the following trivial lemma. 

Lemma 2.11 Assume that (for some value of e) the powers ofT^ satisfy 

'in G N, ||t;^|| < r(n), 

where the function T{n) is strictly decreasing, and T{n) >0. Then the dissipation 

time is hounded from above by 

where T^~^'' is the inverse function ofT. In particular, for the geometric decay T{n) = 
Cr"" with r G (0, 1), C>1, one obtains r^ < ^^^ + 1. 

The upper bound in Eq. ()2.28|) comes from the obvious estimate 

in-i7i|| ^ \\ri II n 
IMe II S ll'-^ell ■> 

on which we apply the lemma with C = l, r=||G'£||. 

To prove the second upper bound, we use the representation of T" in terms of the 

resolvent: 

T," = -^ / X'R,{\)d\ 



27Vi J gr 



valid for any r > rgpiT^). Thus for all r G {rsp{T^), 1), one has 

ll^ni < ^ / \\nR.{\)\\\d\\ < sup ||i?.(A)||r"+^ = -1 

27r Jsr IA|=r d,{T 



r"+\ 



We then apply Lemma 12.111 on the geometric decay for any radius rsp{T^) < r < 1, 



with C = :t^>1. 
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2.4 Dissipation time of not weakly-mixing maps 

In order to better control the growth of Tc-, we need more precise information on the 
noise and the dynamics. In the present section, we restrict ourselves to the dynamical 
property of weak-mixing. We recall [^2] that the map F is ergodic (resp. weakly- 
mixing) iff 1 is not an eigenvalue of Up (resp. iff Up has no eigenvalue) on Lq(T'^). 
We now use Theorem 12.91 in the case where U = Up is the Koopman operator for 
some measure-preserving map F on T'' to establish the following important result. 

Theorem 2.12 Assume that the noise generating density g satisfies the estimates 
(I2.12|) or ()2.14j) with exponent a G (0,2]. If F is not weakly- mixing and at least one 
eigenfunction of Up belongs to H'^{T'^) for some 7 > 0, then T^ has slow dissipation 
time: 



r^ I C- 



Proof. Let h G H'^{T'^) be a normalized eigenfunction of Up with eigenvalue A. 
Applying Proposition 12.21 ii), we get 

||(A-T,)/^|| = ||(/-G,)/.||<ire^^" 

for some constant K > depending on g and h. This implies that ||-Re(A)|| > 
therefore taking the supremi 
in Theorem 12.91 then implies 



therefore taking the supremum over |A| = 1 yields de{l) ^ > -jt-^. The lower bound 






< Tc. ■ (2.31) 

Remark 2.13 Recall that if g satisfies (J2.14J) with exponent a, then the dissipation 
time is also bounded from above, as shown in Corollaru \2.1(A If one eigenfunction h 
has regularity H'^ with 7 > a, then both results imply that the dissipation is simple, 
with exponent a. 
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Remark 2.14 The above results can be stated in more general form: Up does not 
need to be a Koopman operator associated with a map F. The result holds true for 
any isometric operator U on Lq with an eigenfunction of Sobolev regularity. 

The dependence of the lower bound in (j2.31|) on 7 can be intuitively explained as 
follows. In case of non-weakly-mixing maps the eigenf unctions of Up are, in general, 
responsible for slowing down the dissipation. The rate of the dissipation is affected 
by the regularity of the smoothest eigenfunction. In principle, irregular functions 
undergo faster dissipation giving rise to slower asymptotics of r^. It is not clear, 
however, whether the actual asymptotics of the dissipation time will be slower than 
power law in case when all eigenf unctions of Up on Lq{T'^) are outside any space 
H^'{T'^) with 7 > 0. 

The above theorem serves as a source of examples of 'non-chaotic' ergodic dynamical 
systems. A typical example of ergodic but not weakly mixing transformations for 
which this corollary applies is the family of 'irrational' shifts on T"^ i.e. maps Fx = 
a; + c on T*^, where c = (ci, .., q) is a constant vector such that the numbers 1, Ci, .., q 
are linearly independent over rationals. More general and less trivial examples of 
ergodic maps giving rise to a slow dissipation time will be discussed in Section 13.2.21 
(cf. Remark IT^ . 

In Corollarv 12.71 we have shown that for any map F and arbitrary small 6 > 0, the 
pseudospectrum asiT^^) intersects the unit circle for sufficiently small e > 0. If F 
is not weakly-mixing, the spectral radius of T^ (that is, the modulus of its largest 
eigenvalue) is believed to converge to 1 when e ^ 0, and the associated eigenstate 
h^ should converge to a "noiseless eigenstate" h. This "spectral stability" has been 
discussed for several cases in the continuous-time as well as for discrete-time maps on 
T^ [TTllTnn]. 

On the opposite, if F is an Anosov map on T^ (see Section IT^ . the spectrum of T^ 
does not approach the unit circle, but stays away from it uniformly: rsp{T^) is smaller 
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than some vq < 1 for any e > J2I]. Simultaneously, \\T^\\ -^ 1, so we have here 
a clear manifestation of the nonnormality of T^ for such a map. In some cases (see 
|1UUJ and the linear examples of Section EH)), the operator T^ is even quasinilpotent, 
meaning that rsp{T^) = for all e > 0. For such an Anosov map, the spectral radius 
of Tg is therefore "unstable" or "discontinuous" in the limit e — > 0, while in the same 
limit the (radius of its) pseudospectrum as(T^) (for 6 > fixed) is "stable". 
We end this section by determining the coarse-grained dissipation time for non weakly- 
mixing maps. We have 

Proposition 2.15 Let F be a measure-preserving map. If F is not weakly-mixing 
then fc = oo for small enough e > 0. 

Proof. Let h G Lq(T'^) be a normalized eigenfunction of Up, then 

||Ti")/i|| = WG^w^ih+iG^-iMiyWGM-WGeUUGe-iM 

> l~2\\iG,-I)h\\. 

Since the RHS above is independent of n, we see that ||Te || is close to 1 for all times 
and sufficiently small e > 0. ■ Thus as opposed to the noisy case (see Prop. I2.10|) . 
the coarse-grained evolution through a non-weakly-mixing map does not dissipate. 

2.5 Local expansion rate and general lower bound 

We saw in the previous section that there exists no general upper bound for coarse- 
grained dynamics fc- On the opposite, we will prove below a general lower hound for 
both coarse-grained and noisy evolutions, valid for any measure-preserving map F of 
regularity G^. We note that Propositions 12.81 and I2.l5t ) (which are valid indepen- 
dently of any regularity assumption) do not provide an explicit lower bound. 
First we introduce some notation. For any map F E G^, DF{x) is the tangent map 
of F at the point x G T'^, mapping a tangent vector at a; to a tangent vector at 
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F{x). Selecting the canonical (i.e. Cartesian) basis and metrics on T(T'^), this map 
can be represented as a d x d matrix. The metrics naturally yields a norm v G 
Tx{T'^) 1-^ \v\ on the tangent space, and therefore a norm on this matrix: \DF{x)\ = 
max|.„l=i \DF{x) ■ v\. We are now in position to define the maximal expansion rate 
of F: 

HF = limsup pF"||^", where ||/^F"||oo = sup |(DF")(a;)|. 

Since F preserves the Lebesgue measure, the Jacobian Jf{x) satisfies | Ji7'(a;)| > 1 at 
all points. In the Cartesian basis, Jf{x) = det{DF{x)), so that we have ||D-F"(a;)|| > 
1 for all a? G T'^, n > 0. One can actually prove the following: 

Remark 2.16 Although \{DF'^){x)\ and \\DF\\oo tnay depend on the choice of the 
metrics, fiF does not, and satisfies 1 < fJ^F ^ ||-DF||oo- 

From the definition of fiF, for any fi > fiF there exists a constant A > 1 such that 

Vn G N, ||I^i^"||oo < ^At". (2.32) 

In some cases one may take n = fiF ^t^ the RHS. In case /Xi? = 1, HDF^Hoo can 
sometimes grow as a power-law: 

WDF^^KAn^, nen (2.33) 

for some j3 > 0, or even be uniformly bounded by a constant {[3 = 0). 

The relationship between, on one side, the local expansion of the map F and on the 
other side, the dissipation time, can be intuitively understood as follows. A lack of 
expansion (||DF||oo = 1) results in the transformation of "soft" or "long-wavelength" 
oscillations into "soft oscillations", both being little affected by the noise operator 
Ge- On the opposite, a locally strictly expansive map (||DF||oo > 1) will quickly 
transform soft oscillations into "hard" or "short- wavelength" , the latter being much 
more damped by the noise. 
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The following theorem precisely measures this relationship, in terms of lower bounds 
for the dissipation times. 

Theorem 2.17 Let F be a measure-preserving C^ map on T'^, and assume that the 
noise generating density g satisfies i2.1^) or Ii2.14\ ) for some a G (0, 2]. 
i) If \\DF\\ao > 1, resp. fip > 1, then there exist a constant c, resp. constants fi > fir 
and c, such that for small enough e, 

rc> , .||„p|| , ln6-i +c, resp. r, > In 6"^ + c 2.34 

ln(||L)F||oo) m/i 



If F is a C^ diffeomorphism, then \2. 34 ^ holds with \\DF\\ao replaced by \\DF\\ao A 



||D(F~-'^)||oo, resp. with some fi > fip A fip-i . 

a) If \\DF\\oo = 1 then T^ has slow dissipation time, Tc > e^("^i). If the noise 

kernel satisfies the condition \2.14\ ) for a G (0, 1], then the dissipation time is simple, 

Hi) If fip = 1 and \\DF^\\oo grows as a power-law as in Eq. ()2.33|) with /9 > 0, then 
fc > e^("^i)//3. If ||_DF"||oo is uniformly bounded above by a constant, then fc = oo 
for small enough e. 

Remark 2.18 This theorem shows that classical systems on T'^ (i.e. C'^ diffeomor- 
phisms) cannot have a dissipation time growing slower than Cln(e^"'^). In view of the 
results for toral automorphisms (cf. Proposition 4), this lower bound on the dissipa- 
tion time is sharp and consistent with Kouchnirenko's upper bound on the entropy 
of the classical systems, namely all classical systems have a finite (possibly zero) 
Kolmogorov- Sinai entropy (Theorem 12.35. in \ll^ - see also lllf - YHJ )- 

Proof of the Theorem 12. 17L The following trivial lemma (similar to Lemma f2.1H) 
will be crucial in the proof. 



46 Chapter 2. Dissipation time - definition, properties and general results 

Lemma 2.19 Assume that there exists some a > and a strictly increasing function 
'j{n), 7(0) = such that 

\/n>l, IITJ^II > l-e"7(n). (2.35) 

Then the dissipation time is bounded from below as: 

r. > 1^-'^ {^f^) ' (2.36) 

where 7^"-^) is the inverse function of 'j. 

The same statement holds for the coarse-grained version. 

Our task is therefore to bound \\T^\\ (resp. \\Te ||) from below. A simple computation 
shows that for any / G C^(T'^), \\G^f\\ 00 < ll/lloo- Since convolution commutes with 
differentiation, for / G C^ we also have ||V(Ge/)||oo < ||V/||oo- We use this fact to 
estimate the gradient of T^f: 

||V(TJ)|U = ||V(G,f/^/)lU 

< ||V(/oF)|U = ||(V/)oF-DF|U 

< ||(V/) o F\U\DF\\^ = ||V/||oopF||oo. 

Repeating the above procedure m times, we get 

||V(Tr/)|U < \\Vf\UDF\\Z \\V {UpT^f )\\oo < ||V/|UpF|ir'- (2-37) 

We now choose some arbitrary / G Cq{T'^), with ||/|| = 1. We first apply the triangle 
inequality: 

||T;/|| = WG^UpTr'fW > WUpTr^fW - \\{G,-I)UFTr'f\\. 

To estimate the second term on the RHS we use the bound (|2.19|) and the estimate 
(jTTTfl to obtain 

iit:/ii > iiTrvii -c^e"^'iiv/iioopFii^. 



Section 2.5. Local expansion rate and general lower bound 47 

Applying the same procedure iteratively to the first term on the RHS, we finally get 
(remember ||/|| = 1): 

n 

iit:ii > iit;/ii > 1 - ce"^i v/iu Yl w^nz- (2.38) 

m=l 

The computations in the case of the coarse-grained operator are even simpler: 

||rw/|| = \\G,U^GJ\\ 

> 1 _ Ce-^i||V/|U - Ce"^iV(GJ)|UpFl|oo 

> l-2Ce"^iV/||oo||/^i^"||oo. (2.39) 

Notice that from the assumptions on /, || V/||oo cannot be made arbitrary small, but 
is necessarily larger than some positive constant. We choose some arbitrary function, 
say / = Gfc with k = (1, 0) which satisfies || V/||oo = 27r. 

The estimate ()2.38|) has the form given in Lemma 12.191 The growth of the function 
7(n) depends on whether ||DF||oo is equal to or larger than 1, which explains why 
the lower bounds are qualitatively different in the two cases. 

In case ||Z)F||oo is strictly larger than 1, then the function 7(n) grows like an exponen- 
tial, therefore the lower bound is of the type ()2.34|) . For the coarse-grained version, 
a growth of ||DF||oo of the type ()2.32|) yields the lower bound for fc in ()2.34|) . 

In the case ||i5F||oo = 1, 7(^) is a linear function, so that Tc > J'li^fii e~^°'^^\ 

In the coarse-grained version, ii ^p = 1 and HDF^Hoo grows like in ()2.33|) with f3 > 0, 
the dissipation is slow: fc > Ce^*-"^^)/^. In the case where ||DF"||oo is uniformly 
bounded by some constant, the norm of the coarse-grained propagator stays larger 
than some positive constant for all times, so that for small enough noise fc is infinite. 
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2.6 Decay of correlations and general upper bound 

For any two functions f,g & Lq(T'^) the dynamical correlation function for the map 
F is defined as the following function of n G N (see e.g. jT^): 

CfA^) = Cl^{n) = mUU-pg) = {f,U^9) = {P^f,9). 

The same quantity may be defined for the noisy evolution: 

Cl^{n)=m{fT:g). 

We recall that a map F is mixing iff for any /, g ^ L^, 

Cf^giji) — * 0, as n — > oo. 

The correlation function can easily be measured in (numerical or real-life) experi- 
ments, so it is often used to characterize the dynamics of a system. 
To focus the attention, we will only be concerned with maps for which correlations 
decay in a precise way. We assume that there exist Holder exponents ■s*,s G IR+, 

< s^. < s together with some decreasing function T{n) = Ts^s{n) with T{n) >0, 

such that for any observables / G Cq*(T'^), g G Cq(T'^) and for sufficiently small e > 
(sometimes only for e = 0), 

Vn G N, \Cl,{n)\ < ||/||c.. ||(7||c.r(n). (2.40) 

In general, such a bound can be proved only if the map F has regularity C*"^^. The 
reason why we do not necessarily take the same norm for the functions / and g will 
be clear below. 
We will be mainly interested in the following two types of decay 

i) Power-law decay: there exists C > 0, f3 > such that, 

r{n) = Cn-I^. (2.41) 
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This behavior is characteristic of intermittent maps, e.g. maps possessing one or 
several neutral orbits [T^ . 

ii) Exponential decay: there exists C>0,0<a<l such that, 

r{n) = Ca"". (2.42) 

Such a behavior was proved in the case of uniformly expanding or hyperbolic 
maps on the torus (see Section IT^ . as well as many other cases jT3]. 

The central result of this section is a relationship between, on one side, the decay of 
noisy (resp. noiseless) correlations and on the other side, the small-noise behavior of 
the noisy (resp. coarse-graining) dissipation time. The intuitive picture is similar to 
the one linking the local expansion rate to the dissipation: namely, a fast decay of 
correlations is generally due to the transition of "soft" into "hard" fluctuations of the 
observable through the evolution, which is itself induced by large expansion rates of 
the map. Still, as opposed to what we obtained in last Section, the following theorem 
and its corollary yields upper bounds for the dissipation time. 

Theorem 2.20 Let F be a volume preserving map on T'^ with correlations decaying 
as in Eq. \2.40^ for some indices s, s* and decreasing function T{n), at least in 
the noiseless limit e = 0. Assume that the noise generating function g is {[s] + 1)- 
differentiable, and that all its derivatives of order |q:|i < [s] + 1 satisfy 

\D-g{x)\<^, |c.|»l, 

with a power M > d. 

Then there exist constants C > 0, Co > such that the coarse-grained propagator 

satisfies 

Ve<e„ Vn>0, \\fi-)\\<C^^^. (2.43) 
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// the decay of correlations (j2.4(jp also holds for sufficiently small e > (and assuming 
the Perron- Frobenius operator Pp is hounded in C'^{T'^)), then the noisy operator 
satisfies (for some constants C > 0, e^ > Oj; 

Tin) 
Ve<e., Vn>0, \\Tn\<C^^^. (2.44) 

From these estimates, we straightforwardly obtain the following bounds on both dis- 
sipation times (the assumptions on F and the noise generating function g are the 
same as in the Theorem): 

Corollary 2.21 I) If the correlation function satisfies the hound ()2.40|) for e = 0, 
then the coarse-grained dissipation time is well defined (fc < oo^. Moreover, 

i) ifV{n) ~ n~^ then there exists a constant C > such that 

~ d-\-S-\-Sif: 

fc < Ce — 

ii) ifT{n) ~ a" then there exists a constant c such that 

Tc < —Y, r- In e~ + c, 

I ^'^'^1 

II) If Eq. fl2.40|) holds for sufficiently small e > 0, then 
i) ifT{n) ~ n~^ , there exists a constant C > such that 

_d+s+a» 

rc<Ce 8 

ii) ifV{n) ~ 0"", there exists a constant c such that 

d + s + s^ 1 

Tc < — . — ln(e ^) + c. 

I ^^"^l 
Proof of Theorem ?I7M 

P* step: We represent the action of T" (resp. Te ) on an observable / G Lq(T'^) 
in terms of the correlation functions C^ln) (resp. C{n)). To do this we Fourier 



Section 2.6. Decay of correlations and general upper bound 51 

decompose both T^+'^f and /i = Upf, and use Eq. (j2.6j) : 

= E E fiik)9eiJmk){PFej,T:e^)ej. 

(remember that ^ is a real function). A similar computation for the coarse-grained 
propagator yields: 

^i"V= E E fikmJ)9e{k){ej,U-e^)ej. 
Taking the norms on both sides, we get in the noisy case: 

^ E ( E lA(fc)i') E \{PFej,T:e,)\r9.{m{k)\' 
^rr^vf < ii/if E iCp,e_„e.Hri^(ej)^(efc)r, (2.45) 

and in the coarse-graining case 

l|2;(™VlP < ll/f E |Ce_„e,(n)P|^(6J)^(6fc)p. (2.46) 

o^j.fceZ'* 

These two expressions explicitly relate the dissipation with the correlation functions. 



2"'^ step: We now apply the estimates ()2.4()|1 on correlations for the observables e^, 
e_j, PpB^j. In the coarse-grained case, it yields (using simple bounds of the type of 
Eq. dZSZI)): 

Vj, fc G Z'^ \ {0}, |a_,,e,(r2)| < C \JY\kY'V{n). 
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In the noisy case, we need to assume that the Perron- Frobenius operator Pp is 
bounded in the space C'^(T'^). This property is in general a prerequisite in the proof 
of estimates of the type (j2.4(Jj) . so this assumption is quite natural here. 

Vj, fcGZ^\{0}, |C|,^e_„e,WI < C \\PFe^i\\cs\\eu\\c^,T{n) 

< C\\PF\\cAJ\'\k\'*T{n). (2.47) 

We insert these bounds on the decay of correlations in the expressions (J23323EI), for 
instance in the coarse-grained case we get: 

Vn>0, ||Tj")||2<Cr(r2)2/'e-(^+^*) J] \ek\'+''g{ek)A . (2.48) 

3'"'^ step: We finally estimate the e-dependence of the RHS of the above inequal- 
ity. Up to a factor e^'^, the sum in the brackets is a Riemann sum for the integral 
/ \^\'^^'^''9{^)^d^ < oo. A precise connection is given in the following lemma, proved 
in Section! 



Lemma 2.22 Let f G C'^(M"') be symmetric w.r.t. the origin and decaying at infinity 
as |/(a;)| < |a3|^^^ with M > d. Then the following small-e estimate holds in the limit 



e' 



Y, hekf = f fi^fd^ + O(e^). (2.49) 

Let 771 G N satisfy 2m < s + s* < 2m + 2 (notice that m < [s] since we assumed 
-s* < s). From the obvious inequality 

we may replace in the RHS of ()2.48|) the factor |efc|*+** by lefep^ + lefep™"^^. Applying 
Lemma r2.22l to the derivatives of g of order m and m + 1, we end up with the following 
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upper bound, which proves the first part of the theorem: 



i^i"^r < ^r(")'(^l(i^ 



M^ 



< c 



Tin) 



\g\\fJ7n + l- 



The computations follow identically for the case of the noisy operator, yielding the 
second part of the theorem. ■ 



2.7 Dissipation time and optimization problems 

In general the problem of computing the dissipation time is rather complicated. In 
some cases it can be reformulated as an asymptotic optimization problem. To see it, 
one can represent the action of a given unitary operator U in the Fourier basis 

(2.50) 






where for each k 



Oy^k'eZ'i 



(2.51) 



Next we introduce the notation 

C>n[k 

Then for any / G Lq(T'^) we have 



07^fci,...,fc„_ieZ'* 1=1 

{fcoGZMo}:W„(fco,fc„)^0}. 



J2 />o)Trefc„ 

o^fcoezd 



E 

o^k„eZ'' 

E 



X^ /(^o) Y W„(fco,fcn)efc„ 
2 



y^ fiko)l(niko,kn) 

^ f{ko)Uniko,kr 

koeSn(kn) 



(2.52) 



(2.53) 
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The following general upper bound for ||T"/|| holds. 
Lemma 2.23 For any f G Lg(T'^), 

\KJf < E E l/>0)l' E \^n{ko,kr.)f. (2.54) 



For the proof we refer to Section | 

In order to see how this lemma can work in practice let us consider a concrete example. 

To this end we focus on a case when u^./,' is a Kronecker's delta function 

Uk,k' = SAk,k', (2.55) 

where A : Z'^ i-^ Z'' is a linear surjective map. 

Under this assumption the upper bound ()2.54p can be used to obtain an identity for 

||T"||. Indeed, first observe that 

n 

Un{kQ,kn) = Ylge{A^ko)5A"ko,k„ 
1=1 

and hence ()2.54p becomes 

n n 

rr/ir < e i/>o)pn^^(^'^o)<ii/ir max^n^^(^'^) 

On the other hand for any nonzero k E'L'^, one can take in (J2.52J) f = e^ and get 

n 

\\T:fr = l[9M'k) 
1=1 

and therefore 

n 

||T"|| = max YlgJA^k). (2.56) 

Oy^keZ'' \^ 

Let us now determine the class of maps F such that the corresponding Koopman 
operator Up satisfies (|2.55p . The relation (|2.55p implies 



Upek = BAk = e 



2m{Ak,x) 
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On the other hand 

Thus 

(fc, Fx) = {Ak, x) mod 1, Mx eW^, fee Z^, 

that is, A is hnear and A'^ equals the hfting of F from T"' onto M°'. Moreover, the 
matrix A has integer entries and determinant equal to ±1, i.e., A (and F) is a toral 
automorphism. In the next Chapter we will use formula I2.5(il to derive an exact 
asymptotics of the dissipation time for virtually all toral automorphisms. 

2.8 Technical proofs 

Proof of Lemma 12.11 

We use the following upper bound: for any a G (0,2], there is a constant Ca such 
that 

VxeM, 0< l-cos(27rx) <C„|x|". (2.57) 

Besides, one has the asymptotics 1 — cos(x) ~ ^ ^^ small x. We simply apply these 
estimates to the following integral: 

1-^(0 = I {l-cos{2-KX-i))g{x)dx 

< [ C^\x-^rg{x)dx 

< Cal^l" [ \x\''g{x)dx = CaM^\^\''. 
In the case g admits a second moment, we have in the limit ^ — > 0: 

/ {1 - cos{2nx ■ $,))g{x)dx ^ / 2n'^{x ■ $,Yg{x)dx 

« 27r^|^P / {x.^^fg{x)dx, 
where we have used the notation ^ = |l| ioi any ^ 7^ 0. ■ 
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Proof of Proposition 12.21 

The statement i) is standard in tlie context of distributions |125l p. 157]. In our case, 
assume that / G L^ is normalized to unity and consider an arbitrary small 6 > 0. 
Since / G L^(T'^), there exists i^ > s.t. J2\k\>K l/(^)P < ^- Since g is continuous 
and ^(0) = 1, there exists r] such that (1 — 5'(^))^ < 6 ii |^| < rj. Thus using spectral 
decomposition ()2.(j|l of G^, we obtain for all e < -^ 

\\GJ-fr=J2i^-g{ek)r\m\'<6j2 l/>)l'+ E l/>)l'<25- (2-58) 

keZ^ \k\<K \k\>K 

To prove the next statement, first notice that if g satisfies the estimate ()2.13|) for 
the exponent a, it also satisfies it for the exponent 7 A a. Using once again spectral 
decomposition of G^, and applying the estimate ()2.13p with the latter exponent we 
get 

WGJ-ff < J2{C,^^M,^.\ek\'r^^f\f\k)\' 
keZ'i 
< (C,^„M,^„)2e2(7Aa) J2 |fcp(^^")|/(fc)p (2.59) 

^ (^7Ao^7Aa) ^ II / ll/f^Aa. 

To obtain the last statement, we notice that any / G C^(T'^) is automatically in 
H^(T'^), and that its gradient satisfies 

iiv/iiL > iiv/ir = 47r^ E ifcri/»r > ^^' E i^r^'^-v^r. 

The inequality ()2.59j) with 7 = 1 then yields the desired result. ■ 

Proof of Proposition 12.51 

The proof will be based on the Riesz convexity theorem (see |131j . pp. 93-100) which 
states that for any operator T defined on Lp{T'^), 1 < p < 00, ln||T||p is a convex 
function of p^^. On the space Lp(T'^) we consider the operator T := T^^a — (■) and we 
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have the relation T"/ = T^aif ~{f))'^f ^ LP(T^),n > 1 because T^^a is conservative. 
Now since ||/ - {f)\\p < 2||/||p, it follows that 

\\f"% < 2||t;;j|p,o < 2 (2.60) 

\\t:j,,o < \\n\, (2.61) 

for 1 < p < 00, n > 1. The Riesz convexity theorem implies that if p < g < oo 

In llf "IL < ^ In ||r"|L +(1-^1 In llf "IL (2.62) 

q \ qJ 

while if 1 < g < p 

ln|rL< ([^)ln||fX+(l-i^)lnl|rii. (2.63) 

From (|2.62|) - (j2.63p we have the interpolation relations 

||T"||g < ||T"||P/^||f"||;^P/'?, p<q <oo (2.64) 

IIT"II < ||T"l|(l~9-i)/(l-p-')||'fn||l-(l-9"')/(l-p-') T ^ n <^ n (9 a^'] 

Ik Ik — Ik lip IK 111 ) J^ ^ y ^ F y^.vo) 

which, along with ()2.60p - ()2.6H) . imply 

iit::ji,,o < 2\\t:Xi p<q<oo 

WT-JUo < 2\\T-X,'^-'^^^'-'-'\ l<q<p 

This proves that the order of divergence of Udissip) are the same for 1 < p < oo. 
Estimates ()2.64|) - ()2.65p also show that the order of divergence oiudissi^) and ridissico) 
is at least as high as ridissip), I < p < oo. ■ 

Proof of Proposition 12.71 

We prove the limit d^{l) >0 by contradiction. Assume that there is some constant 

a G (0,1) such that for all e > 0, the distance d^{l) > a. We will show that the 
following triangle inequality holds: 

Ve>0, d,{l-a/2)> a/2. (2.66) 
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First of all, notice that the assumption d^{l) > a means that for any X E S^, 
\\R^{X)\\ < a"^. We apply the following identity |125j : 

i?,(A') = R,{\){1 + J](A - XTReiXT} 

n>l 

with A' = rX, for 1 — a < r < 1. Taking norm of both sides yields the bound 
||-Re(A')|| < r-(i-a) ' Uniformly w.r.t e. Since this upper bound holds for any |A'| = r, it 
shows that the spectral radius rsp{T^) < 1 — a, and proves ()2.(i(j|l by taking r = 1 — a/2. 
We can now use ()2.fi(i|l in the upper bound ()2.29j) of Theorem 12 .91 this e-independent 
upper bound shows that r* remains finite in the limit e ^ 0, which contradicts 
Proposition 12 .81 ■ 

Proof of Lemma 12.221 

Considering its decay at infinity, the function / is automatically in L^(R'^). The 
function /^ is the Fourier transform of the self-convolution /*/. Therefore, using the 
parity of / and applying the Poisson summation formula to the LHS of ()2.49|1 yields 

6^5^/>fc)^=/f(^K+ Y: (/*/)(7)- (2-67) 

A simple computation shows that (/ * f){x) also decays as fast as |a;|~^''^. This piece 
of information is now sufficient to control the RHS of ()2.67|1 . yielding the result, 
Eq. TIJ^ . m 

Proof of Lemma IZ^ 

Using the notation introduced in Section I^TI one has 
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We note that for any n and fc„ G U^, the sequence lAn{kQ,kn) (indexed by feg G U^) 
belongs to P{IJ^). Indeed, using the Cauchy-Schwarz inequahty and identity (|2.51|) 
one gets for ra = 2, 

2 

-e(|fclP" + |fc2|2«) 



OT^fcieZ'* 



5^ |W2(fco,fc2)r= 5^ 

where K denotes a constant which depends only on e and a. Similar estimates hold 

for n > 2. 

Now applying the Cauchy-Schwarz inequality in ()2.52p we get 

\KJ\\" < E E l/>o)l' E \^n{ko,K)\'. (2.68) 



OjLk„eZ'i ko&Snikn) 



fcog5„(fc„) 
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Chapter 3 

Dissipation time of classically 
chaotic systems 



3.1 Dissipation time of toral automorphisms 

3.1.1 Preliminaries 

It is well known (see pQ) that (the lifting map corresponding to) any toral homeo- 
morphism H : T'^ \-^ T'^ can be decomposed into three parts H = L + P + c, where L, 
the linear part, is an element of SL{d, Z) - the set of all matrices with integer entries 
and determinant equal to ±1, P is periodic i.e. P{x + v) = P{x) for any v G Z*^, 
and c is a constant shift vector. 

Every algebraic and measurable automorphism of the torus is continuous. Each con- 
tinuous toral automorphism is a homeomorphism with zero periodic and constant 
parts and hence can be identified with an element of SL{d, Z). And vice versa, each 
element of SL{d, Z) uniquely determines a measurable, algebraic toral automorphism. 
Thus from now on the term toral automorphism will simply be reserved for elements 
of SL{d, Z). We recall here that all Anosov diffeomorphisms on T*^ are topologically 
conjugate to the toral automorphisms ([2^1, |HS1)- 
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Below we summarize a few definitions from ergodic theory along with some well known 
ergodic properties of toral automorphisms (cf. jHSl P- 160, [HZj and [TU]). 

Definition 3.1 ^4 dynamical system (T'^,/i, F) is called a K-system, (possesses K- 
property) if there exists subalgebra A of the algebra M. of all ^-measurable sets such 
that 

• z ) V ^ C F^'A 

ngZ+ 

• ii) V ^"^ = M 

n>0 

• ii) A F-''A = C 

n>0 

where C denotes the algebra of sets of measure O or 1. We note that the definition 
easily extends to abstract noncomutative version (in this case C stands for c\). 

Proposition 3.2 Let F be a toral automorphism. 
The following statements are equivalent 

a) no root of unity is an eigenvalue of F. 

b) F is ergodic. 

c) F is mixing. 

d) F is a K-system. 

In the sequel we will use the following result (cf. [126_ ). 

Proposition 3.3 The entropy h{F) of any toral endomorphism F is computed by the 
formula 

h{F)= 5^1n|A,|, (3.1) 

|A,I>1 

where Xj denote the eigenvalues of A. 

From the formula p.lj) one immediately sees that a toral automorphism has zero 
entropy iff all its eigenvalues are of modulus 1. In fact much stronger result holds (for 
proof see Section IT3j) . 
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Proposition 3.4 A toral automorphism has zero entropy iff all its eigenvalues are 
roots of unity. In particular all ergodic toral automorphisms have positive entropy. 

Given any toral automorphism F we denote by P its ciiaracteristic polynomial and by 
{Pi, ..., Ps} the complete set of its distinct irreducible (over Q) factors. Let dj denote 
the degree of polynomial Pj and hj the KS-entropy of any toral automorphism with 
the characteristic polynomial Pj. For each Pj we define its dimensionally averaged 
KS-entropy as 

hi 



_ J2 
^^ d, 



hj = -f. (3.2) 



Definition 3.5 Assuming the above notation, we define minimal dimensionally av- 
eraged entropy of F (denoted h{F)) as 

h{F) = min hj 

j=l,...,s 

3.1.2 Main theorems 

In this section we state two main theorems of Part I of this work. Both theorems 
concern the asymptotics of the dissipation time for toral automorphisms in arbitrary 
dimension. Our main task is to derive not only the logarithmic order of the asymp- 
totics in the case of chaotic maps but also to find the exact value of the constant 
and relate it to dynamical properties of the map via the connection with the minimal 
dimensionally averaged KS-entropy of its irreducible blocks. 

In order to be able to perform exact calculations we need to specify more concretely 
the family of noise kernels we are going to work with. Namely, we assume here that 
the noise kernel is a-stable, cf. fl2.15|) . that is 

^,,„(a.):=^e-l^'=l)"efc(x), (3.3) 

fcezd 

Under this assumption we have the following results regarding the dissipation time of 
||T"^|| (in order to emphasize that our calculations depend in an essential way on the 
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choice of the noise kernel we denote noisy operators associated with a-stable kennels 
by of T,_„) 

Theorem 3.6 Let F be any toral automorphism, Up the Koopman operator associ- 
ated with F , G^^a a-stahle noise operator and T^^ = G^^aUp. Then 
i) T^^a has simple dissipation time iff F is not ergodic. 
a) T^^a has logarithmic dissipation time iff F is ergodic. 

Hi) IfT^^a has logarithmic dissipation time then the dissipation rate constant satisfies 
the following constraint 

1 ^ 1 



h{F)~ ~ h{F) 
where h{F) is a positive constant satisfying h{F) < h{F). 

The natural question arises, whether the lower bound for the dissipation rate constant 
given in the above theorem is best possible. The next theorem and its corollary 
provides a strong argument in favor of this conjecture. 

Theorem 3.7 If F is ergodic and diagonalizable then 

Tc^J—ln{l/e). 
h{F) ^ ' ' 

That is, the dissipation rate constant ofT^^a is given by 

1 



lir 



HF) 

In this case ergodicity of toral automorphisms is equivalent to hyperbolicity (in higher 
dimension it is not true, since there exists so called quasihyperbolic automorphisms, 
which although ergodic, possesses eigenvalues of modulus one P^llUSj ). Two dimen- 
sional hyperbolic toral automorphisms are usually referred to as the cat maps. 
Using Corollary 13.171 and applying Theorem 13.71 in two and three dimensions one gets 
the following 
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Corollary 3.8 (Cat maps) Let F be any ergodic, two or three dimensional toral 
automorphism. Then 

Tc ~ ln(l/e), 

h{F) ^ ' ' 

We end this section with a remark that toral automorphisms provide a good example 
on which robustness of the dissipation time (i.e. its independence of unimportant de- 
tails of the underlying conservative dynamics) can be tested. To this end we compare 
the dissipation time with another characteristics of chaoticity - the decay of correla- 
tions. As we will show in Section KlHI the dissipation time has the same logarithmic 
asymptotics among general class of Anosov diffeomorphisms while the decorrelation 
may be exponential (generic Anosov case) or super-exponential depending on partic- 
ular map. 

Indeed, we illustrate this fact below by the following result (obtained in the simi- 
lar way as the above theorems) on the decay of correlations for c?- dimensional toral 
automorphisms (for a proof see Section IT^ . 

Proposition 3.9 Let F he a diagonalizable ergodic toral automorphism, Up its Koop- 
man operator and A any constant such that < A < h{F). Then for any f,h E 
Lq(T^'^) the correlation function for noisy dynamics generated by F and and any 
a-stable noise decays superexponentially i.e. 



C}M = {f,T:,^h)<\\f\\\\h\\e 



-e"A' 



a \cy.n 



Moreover, let f,hE G ^{L'^i^'^'^)) be smooth observables, where G^ denotes Gaussian 
noise operator. Then the decay of correlations of unperturbed operator Up i-G- of toral 
automorphism itself is still superexponential 



Gf^u{n) = {f^Uph) <\\G:'f\\\\G:'h\\e 



e2A2" 
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3.1.3 Asymptotic arithmetic minimization problem 

In this section we prepare the ground for the proof of theorems stated in previous 
section. To this end we need to derive a concrete version of the general formula 
(|2.56|) obtained in Section IT7I for any toral automorphism and arbitrary noise kernel. 
Using the fact that here we consider only a-stable kernels the formula (j2.56|) can be 
rewritten as follows 

n 

||T"|| = max n^eU'fc) = e-^"°^'Vfc«dEr.i|A'fcr_ .g^^N 

Hence, for toral automorphisms and a-stable kernels, the calculation of the dissipation 
time reduces to the following nonlinear, asymptotic (large n) arithmetic minimization 
problem 

n 

min VM'fcl", (3.5) 

where A G SL{d,7j). When A is not ergodic the asymptotics of (j3.5|) is clearly of 
the order 0{n). Thus we will be only concerned with the ergodic case. For d = 2 
the problem (j3.5|) can be solved easily as follows. Consider first the case that A is 
symmetric and a = 2. From det(yl) = 1 we see that eigenvalues are A, A~^ with 
|A| > 1. We have 

2n+l n 



mm 

1=1 l=—n 



mm 

0^k=ki+k2&'i 



mm 

Oj^keZ'i 



y \A^k\' = min y \A^k\' 

l=—n 
n n 

fc|2 + Y, iAr'ifci|2 + \\\~^'\k2? + Y. i^r'^i^ii' + i^i'^i^^ 
1=1 

n 

^|A|2'|fc|2= ^|Ap'. 



2 



1=1 1=1 



l=—n l=—n 

Hence there exist constants Ci and C2 such that 



Cie'^(^)"< min V U'fcP < C2e'^(^)". 
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where h{A) denotes the KS-entropy of A. The estimates for the general case of 
non-symmetric A and a ^ 2 are similar. 

In higher dimensions, the solution to ()3.5|1 is much more involved because of the 
presence of different eigenvalues with absolute values bigger than one. We have the 
following general estimate 

Theorem 3.10 Let A G SL{d,7j) be ergodic. There exist constants Ci and C2 such 
that for any < 5 < 1 and sufficiently large n 

n 

fj^^{i-5)ah{A)n ^ miu V I A'fcl" < Cane"'^^^)" (3.6) 

where as before h{A) denotes minimal dimensionally averaged entropy of A and h{A) 
denotes a constant satisfying < h{A) < h{A), with equality achieved for all diago- 
nalizable matrices A. 

The question whether the equality h{A) = h{A) holds for all ergodic matrices remains 

open. 

The proof of the theorem relies on nontrivial use of three number-theoretical results 

stated below. 

I. Minkowski's Theorem on linear forms 

Let Li,...,Lii be linearly independent linear forms on M'^ which are real or occur in 

conjugate complex pairs. Suppose ai,a2,...,ad are real positive numbers satisfying 

aia2...a(i = 1 and ai = Qj, whenever Lj = Lj. Then there exists a nonzero integer 

vector k &l/ such that for every j = 1, ...,d, 

\Ljk\ < Daj, (3.7) 

where D = | det[Li, ..., Lj]^/'^. 

Minkowski's Theorem on linear forms will be used to obtain a sharp upper bound on 

the asymptotic solution of the arithmetic minimization problem. The proof of the 
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above theorem and its generalization to arbitrary lattices can be found in P^ (Chap. 

VI). 

II. Schmidt's Subspace Theorem 

Let Li, ..., Lrf be linearly independent linear forms on W^ with real or complex algebraic 
coefficients. Given 6 > 0, there are finitely many proper rational subspaces ofM."^ such 
that every nonzero integer vector k with 

d 

Y[\Ljk\ < \k\-^ (3.8) 

i=i 
lies in one of these subspaces. 

Schmidt's Subspace Theorem will be used in conjunction with Van der Waerden's 

Theorem on arithmetic progressions (see below) to obtain a sharp lower bound for the 

asymptotic solution of the arithmetic minimization problem. The proof of Schmidt's 

Subspace Theorem can be found in pTTTj (Theorem IF, p. 153). 

Definition 3.11 For a given set of linear forms and for fixed 6 > 0, the smallest 
collection of proper rational subspaces of M.'^ which contain all nonzero integer vectors 
satisfying i3.8^} . is called the exceptional set and denoted by Es. 

A main difficulty to be resolved in using Schmidt's Subspace Theorem is to show that 
the minimizer of either the original problem ()3.5p or an equivalent problem does not 
lie in the respective exceptional set which is in general unknown. We will pursue the 
latter route by using Van der Waerden's Theorem on arithmetic progressions to show 
that one can always construct an equivalent minimization problem whose minimizer is 
guaranteed to lie outside the corresponding exceptional set. To this end we note that 
Schmidt's Subspace Theorem is true when the standard lattice Z*^ is replaced by any 
other rational lattice, that is any lattice of the form A = Q{Z'^) where Q G GL{d, Q). 
Schmidt's subspace theorem can be generalized to this situation by considering the 
set of new forms Lj = LjQ. The fact that Q G GL{d, Q) implies immediately that Lj 
are still linearly independent forms on R'^ with real or complex algebraic coefficients. 
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III. Van der Waerden's Theorem on arithmetic progressions 

Let k and d he two arbitrary natural numbers. Then there exists a natural number 

n^,{k, d) such that, if an arbitrary segment of length n > n^ of the sequence of natural 

numbers is divided in any manner into k (finite) subsequences, then an arithmetic 

progression of length d appears in at least one of these subsequences. 

The original proof was published in |124j : Lukomskaya's simplification can be found 

in IH3I. 

Before presenting the proof of our main results we state a number of technical facts 

concerning the structure of toral automorphisms. 

3.1.4 Algebraic structure of toral automorphisms 

In this section we denote by GL{d, Q) the group of nonsingular d x d matrices with 
rational entries or the group of linear operators on Euclidean space M'^, which are 
represented in standard basis by such matrices. We generally use the same symbol 
to denote both operator and its matrix. 

In the sequel a vector a; G M'^ will be called an integer (or integral) vector if all its 
components are integers, and similarly a rational, an algebraic vector if all its com- 
ponents are rational or respectively algebraic numbers. The term rational subspace 
of W^ will then refer to a linear subspace of W^ spanned by rational vectors (cf. |117j 
p. 113). 

Definition 3.12 A G GL{d,Q) is called irreducible (over Q) if its characteristic 
polynomial is irreducible in Q[x]. 

Lemma 3.13 The following statements about a matrix A G GL{d,Q) are equivalent. 

a) A is irreducible. 

b) A does not possess any proper rational A-invariant subspaces ofW^. 
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c) No rational proper subspace ofM.'^ is contained in any proper A-invariant subspace 
ofR'^. 

d) For any nonzero q G Q"' and any arithmetic progression of integer numbers 
ni,...,nd, the set {A'^'^q,A^^q,...,A'^''q} forms a basis ofM.^. 

e) A'^ is irreducible. 

f) No nonzero q & Q^ is orthogonal to any proper A-invariant subspace o/M"^. 

g) No proper A-invariant subspace ofW^ is contained in any proper rational subspace 
ofR'^. 

Definition 3.14 We say that operator A G GL{d, Q) is completely decomposable over 
Q if there exists a rational basis ofM.^ in which A admits the following block diagonal 
form 

Ai ... 
A2 ... 



(3.9) 



... Ar 

where for each j = 1, ..., r < d, Aj G GL{dj, Q) is 



irreducible and X]f=i ^j = d- 



In general, any matrix A G GL{d, Q) admits a rational block diagonal representation 
[^j]j=i,...,r- The smallest rational blocks to which A can be decomposed are called 
elementary divisor blocks. The characteristic polynomial corresponding to any ele- 
mentary divisor block is of the form p"*, where p is an irreducible (over Q) polynomial 
(see, e.g., [IE])- Although elementary divisor blocks cannot be decomposed over Q 
into smaller invariant blocks, some elementary divisor blocks may not be irreducible. 
This happens iff ?7i > 1 iff A is not completely decomposable over Q. One has the 
following elementary fact (see Section IX^ for a proof). 



Proposition 3.15 A G GL{d, 
nalizable. 



is completely decomposable over Q iff A is diago- 



(3.10) 
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However, even ii A E GL{d, Q) is not completely decomposable, each elementary 
divisor block of A can be uniquely represented (in a rational basis) in the following 
block upper triangular form 

^B C 
D 

where B is the unique rational irreducible sub-block associated with A-invariant ra- 
tional subspace of that elementary divisor and C, D denote some rational matrices. 

Proposition 3.16 All the eigenvalues of an irreducible matrix A G GL{d,Q) are 
distinct (complex) algebraic numbers. In particular all irreducible matrices are diag- 
onalizable. 

The proofs of the above propositions can be found in Appendix B. 
Finally we note that since the leading coefficient and constant term of a characteristic 
polynomial of any toral automorphism are equal to 1, the only possible rational 
eigenvalues of such map are ±1 or ±i. The latter fact implies that ergodic toral 
automorphisms do not possesses rational eigenvalues. Thus we have the following 

Corollary 3.17 Let F be an ergodic, two or three dimensional toral automorphism. 
Then F is irreducible (and hence diagonalizable) . 

3.1.5 The Proof of Theorems EIH EI3 and lOTll 

This section is entirely devoted to the proofs of the main theorems of this chapter. 
We start with the proof of Theorem 13.101 which constitutes the main ingredient in 
the proofs of Theorems 13.61 and 13.71 

Let [Aj]j=i....,r be a rational block- diagonal decomposition of A into elementary divisor 
blocks. Since A G SL{d.,'L)., there exist a transition matrix Q G SL{d,Q) such that 
for every Z G Z, 

A^ = Q-\[A,]yQ 
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and moreover each elementary divisor block [A]j is represented in its block upper 
triangular form (j3.1(jp . 

The matrix Q defines a new lattice A = Q(Z'^) and acts bijectively between this 
lattice and the standard lattice Z'^. Hence 

n n n 

min VU'fc|"= min V |g-i(M,-l)'gfcr = min V |g-VM.l)'g|". 
Moreover 

\\Q\r\i[A,]yQr < \Q-\[A,]yqr < \\Q-Tm,]yqr^ vz,j,a. 

Now we decompose A into the direct sum of lower dimensional sublattices Aj corre- 
sponding to invariant blocks [Aj]. So that 

n n 

min V|([A,])'qr = . min min 5^|(A,)'qr- (3.11) 

1=1 J \ , , S ^ J i^^ 

Thus, without loss of generality, we may specialize to the case that A is already 
indecomposable over Q i.e. A does not possesses any proper elementary divisor blocks. 
To simplify the notation we will work with the standard lattice A = Z'^. According 
to the remarks following the statements of Minkowski's and Schmidt's Theorems the 
proof can be easily adapted for any rational lattice A = Q{Jf-\ 
Since the technique of the proof differs depending on diagonalizability of A we consider 
two cases: 

Diagonalizable case 

Here we concentrate on the case when A is diagonalizable and hence due to its in- 
decomposability irreducible (cf. Proposition 13. 15() . 

We denote by Aj (j = 1, ..., d) the eigenvalues oi A. Following Proposition l3. 16l we note 
that \j are distinct (possibly complex) algebraic numbers and hence there exists a 
basis (of C^) {fj}j=i,...,d composed of normalized algebraic eigenvectors corresponding 
to eigenvalues Aj. 
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We denote by [-Pj]^=i the projections on [Vj], and by [Lj] the corresponding hnear 
forms. It is easy to check that [Lj] are given, in the Riesz identification, by the 
eigenvectors [uj] of the matrix A'' which are co-orthogonal to [Vj], i.e., {ui,Vj) = 
for i ^ j. [uj] and [vj] are real or occur in complex conjugate pairs. We have 

d d d 

x = Y^ PjX = J2iLjx)vj = ^(a;, Uj)vj, Wx G M^. 

The equivalence between any two norms in a finite dimensional vector space, implies 
the existence of absolute constants Ci, C2 such that 

d d 

Using the above inequalities, the monotonicity of a map x h^ x'^ and an obvious 
inequality (a + 6)" < a" + 6", which holds for all positive a,b and a G (0, 1] one 
obtains the following estimates 

n n / d \ n / d 

1=1 1=1 \ j=i J 1=1 \j=i 

n d d / n \ 

and on the other hand 

n / n \ '-' / n d 

j2\A^kr > ij:\A^k\A > ij2c,j2\p^A^f^\ 
1=1 \i=i J \i=i j=i 

(n d \" / d / n \ N" 

j:j:iA,i-ip,fcio =cAy^[y^\hA\pm' 
i=\ i=i / \i=i \i=\ ) 

Now we introduce some notation 

Aj := max{l,|Aj|}, (3.12) 

/ d \ i/<i 
Vo := fn^^j • (3.13) 
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One can easily observe that there exists a constant C such that 

n 

1=1 
In the sequel we do not distinguish between particular values of constants appearing 

in computations. The symbols Ci,C2,.. are used to denote any generic constants 

independent of n. 

The normalization condition \vj\ = 1 implies the following relation 

\Pjx\ = \Ljx\. (3.14) 

Combining the above estimates one gets the following general bounds 

(d \ " n d 

^Af|L,fcn <^|A'fc|"<C2n^A^""|L,fc|". (3.15) 

j=i J 1=1 j=i 

Therefore in order to estimate ()3.5|) it suffices, essentially, to estimate 

d 
min VAriL.fcl". (3.16) 

We denote by Zn the sequence of minimizers i.e. nonzero integral vectors solving 

duni). 

Upper bound. 

For the upper bound we assign to the set of linear forms Lj the set A composed of 
all real vectors a = (ai, ...,ad) satisfying the conditions aj > 0, for j = 1, ...,d and 
tti = ttj whenever Lj = Lj and 

d 

l[aj = l (3.17) 

From Minkowski's theorem on linear forms, we know that for any a E A, there 
exists nonzero integral vector ka satisfying \Ljka\ < Daj, j = l,...,d, where D = 

|det[Li,...,Lrf]|V'^. 
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Thus 

d d 

Y, >^T\^jkX <dY, \""«i • (3-18) 

The minimizing property of 2„ imphes that for any a E A, 

d d 

J2 Ari^.-^^r < Yl ^^ti^m''- (3-19) 

Thus combining (jH.lSp and ()H.19|1 . and applying the Lagrange multiphers minimiza- 
tion with the constraint (jH.17|) (and using the fact that Aj = \j whenever Li = Lj), 
we get 

d d / d \^/'^ 

J2 >^T\^jzX <DminY \""«J = ^^ n ^"" = ^^^?eo"- (3-20) 

Thus the following upper bound holds 

n 

min V U'fcr < C2n\lZ- (3.21) 

Lower bound. 

Let m denote an arbitrary natural number. Using the fact that A acts bijectively on 
Z'^ we can restate the minimization problem ()3.16jl in the following form 

d d 

min VAriL.fcr = min V AriL^^-'^A'^fcr (3.22) 

d 

= min VAr|A,r°"'|L,A™fc|° (3.23) 

That is 

d d 

^AJ"|LjZ„|° = ^AJ"|Aj|""'"|LjA™z„|". (3.24) 

i=i i=i 

We choose arbitrary 6 > and consider the exceptional set Es (see Definition 13.111) 
associated with the system of linear forms [Lj]. Since [Lj] correspond to the eigen- 
pairs [Xj, Uj] of A'^ they are linearly independent linear forms with (real or complex) 
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algebraic coefficients. Thus the subspace theorem asserts that E^ is a finite collection 
of proper rational subspaces of M*^. We denote by k^ the number of subspaces forming 

Now we want to show that for all sufficiently large n there exist an integer m < n 
such that A'^Zn does not lie in any element of Es- To this end we assume to the 
contrary that all A^Zn lie in the subspaces forming Es and we divide the sequence 
of natural numbers 1, ...,n into ks classes in such a way that two numbers mi and 
7712 are in the same class if A^^Zn and A^^Zn lie in the same element of E(,. Now 
let 7i^,{ks,d) be the number given in the van der Waerden theorem and let n > n,,. 
Then there exists an arithmetic progression 7ni, ...,7nd in one of these subsequences. 
By Lemma f3. 131 d) the set of vectors {A'^^Zn, A^'^z^, ..., A"^'^Zn} forms a basis of the 
whole space M*^, which contradicts the fact that they lie in one fixed rational proper 
subspace. Hence for any 6 > and n>n^, there exists m* < n such that A"^* Zn does 
not lie in any element of E^. 
Now, introducing the notation 

Zn = A^'Zn (3.25) 

one concludes from ()3.24|) that for any 5 > and all n > n^, the following equality 
and estimate hold 

d d 

J]A-"|L,;2j" = 5^Af1A,r""*|L,£„r (3.26) 

(3.27) 



with some / : 



;=1 


i=i 


d 

_\\LjZn\ > 


1 


i=i 




■ be rewritten as 




d 


1 


/(^nl)^ 


R+ such that /(r) 


<r,'ir> 0. 



(3.28) 
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Using (J3.25J) and (j3.2Up we obtain the existence of a constant A > 1 such that 

J\\^n\) ^ \Zn\ — 1^ Zn\ ^ ■^maxl^n] 
d 

< Ka. Yl \"I^J-^"I ^ dDiXmaAeoT < A"- (3.29) 

Note that Hj A^ = 1. So, by (jT^ the quantities 5j-„ = (|Aj|-'^*/(|2;„|)'^/'^|Lj2;„|)", 
j = 1, ...,d satisfy the constraint 

d 

Y[Bj^n = l, \/n>n^. (3.30) 

i=i 

Thus applying ()3.29|) and the Lagrange multiphers minimization with the constraint 
()3.30|) one gets 

d d 

E\ara|\ l—arritlT z. |ck fflZ. \\^ctS/d\ "" \an ry ^ \—2anS/d\an . \an(l—5) 

\j \\j\ \LjZn\ =/(|z„|j 2_^X,j i^j,„ > A Agg„ =. Agg„^ >. 

i=i i=i 

This and equahty ()3.26|) yields the following lower bound for ()3.16|) 

min VAriL.fcr > A^""^!-''). (3.31) 

Non-diagonalizable case 

We move on to the general case where A is not irreducible (but, as assumed at 
the beginning of the proof, indecomposable over Q). We denote by B the invariant 
irreducible sub-block of A given by its block upper triangular decomposition (j3.1(Jj) 
and by S the rational invariant subspace associated with this block. We note that B 
as an irreducible matrix is diagonalizable. 

Upper bound. 

Note that 

n n 

minVU^fe|°< min V |5'fc|°. (3.32) 
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The corresponding upper bound (j3.2ip for B is then also an upper bound for the 
whole matrix A. We note that geometric average of Xj over S is equal to the geometric 
average of all Xj associated with matrix A (i.e. over the whole space W^). 



Lower bound. 

According to our assumption A is indecomposable and thus the characteristic poly- 
nomial of A is of the form p™ for some irreducible p. All Jordan blocks of A have the 
same size m and different Jordan blocks correspond to distinct eigenvalues. Denote 
by b the number of the Jordan blocks in A and by Xj, where j = l,..,l all these 
distinct eigenvalues . Since each Xj has algebraic multiplicity m, we get d = mb. Let 
{vjh}j=i^i,.h=o,...,m-i be a basis (of C^) in which A admits the Jordan canonical form. 
As usually Lj ^ will denote the corresponding linear forms. Each Vj^^ can be regarded 
as a generalized eigenvector of A associated with an eigenvalue Xj. We assume that 
these generalized eigenvectors are ordered according to their degree i.e. Vj^h satisfies 
the equation {A — XjI^^^Vj^h = 0. Reordering the eigenvalues, if necessary, we can 
also assume that Ai has the largest modulus among all eigenvalues of A and hence 
Ai = |Ai|. Let Zn be the sequence of minimizers solving ()3.5|) . We first note that 
for each n there exists < h < m — 1 such that Li^h^n 7^ 0. Indeed, otherwise for 
all h = 0,...,m — 1, Li^hZn = and consequently for any n and h Li^hA"'Zn = 0. 
The latter implies that the set of consecutive iterations {zn,A^Zn,A'^Zn, ■■■} spans a 
proper rational A-invariant subspace of W^ which does not have any intersection with 
the subspace spanned by the generalized eigenvectors of A associated with eigenvalue 
Ai. This clearly contradicts the irreducibility of p. Now, for given n we denote by 
h{n) the biggest index h for which the condition Li^h^n 7^ holds. 
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We have the following estimate 

I b m—1 m~l~h / -. \ 

I ^ — , .^ — , .^ — , . I n \ \ 

(3.33) 



b m—1 
7-=l h=0 



ri—i — ii / \ 



J-'j,h+i^n 



n 



< Ci|A"2„r<Ci5^|A'z„r<C2nA^:„, (3.34) 

1=1 

where the last inequality follows from previously established upper bound. 

From the Diophantine approximation and the assumption that \Li^h{n)Zn\ 7^ 0, there 

exists P > such that (see |117j p. 164) 

\Ll^h(n)Zn\ > I 1^- (3.35) 

Thus combining ()3.33|) with ()3.35|) one gets 

\ cm I _ I —a/S ^ \an\ T ^ l"^/^^\"" 

After rearrangements one obtains the following lower bound estimate for ()3.5|) 
C 



A"" < C\znV < V lA^zX = min V U'fcr, (3.36) 

where 



n z_. o^fc,^. .^^^ 



A ^'^' 

We note that ergodicity of A implies Ai > A^eo > 1 (see ()3.13p . ()3.1|) and Proposition 
13. 4|) which ensures non-triviality of this lower bound. 

Now in order to finish the proof is suffices to combine the estimates ()3.2H) . ()3.31|) and 
()3.36|) . and note that 

\an a ^ ' n ah(A)n 

^geo ~ ^ ~ ^ 

which yields ()3.6p . ■ 
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Proofs of Theorems |M1 and UTTI 

We start with part i) of Theorem 13 . 61 which follows as a simple consequence of Theorem 
12.121 Indeed, it suffices to construct an eigenfunction of Up which belongs to L^iT'^) fl 
H°(T'^). Directly from Proposition 13.21 one concludes that F, and hence also A, 
possesses a root of unity in its spectrum. This means that A'^ko = ko, for some m 
and certain nonzero vector fco, which can be chosen to be an integer. Now we define 

/ = /fco ~^ JAko + ••• + /A™-ifco 

Obviously / G -^oC^*^) ni7"(T'^), for any a. To complete the proof it suffices to notice 
that 

Upf = f Ako + f A^ko + ••■ + f A^ko = /fco + f Ako + ••• + f A"^-'^ko = /• * 

Now we apply Theorem I3.1UI to prove part ii) and iii) of Theorem 13.61 and Theorem 

In order to determine the dissipation time of T^^q, one has to determine the asymptotics 
of IITJ'q,!! when n goes to infinity. According to formula (|3.4j) this problem reduces to 
problem (j3.5|) which has been solved by Theorem 13. lUp . This in view of this theorem 
there exist constants Ci and C2 such that for any 5, 5' > and sufficiently large n 

n 

o^fceZ'* -^ 
Using formula (j3.4j) 

^ — II-' e, all — ^ 

Now when n = Tr, we have 



C-l 



and 



fj^^{l-S)h(A)Tc < 1 < (j^^{l+5')h{A)T, 

~ e ~ 



ln(l/e) - In C2 ) < r, < J-^— - ( ln(l/e) - In Ci , , 



l + 6')h{A)\ J ~ ~ il-S)h{A) 
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which proves part ii) of Theorem 13.61 i.e. the logarithmic growth of dissipation time 

as a function of e~^. 

Moreover, using the definition of dissipation rate constant 

Rr = hm 



.-0 ln(l/e) 
we obtain 



<Rr< 



il + S')h{A) ~ ~ {l-6)h{A) 
Finally letting S ^ and 6' —>■ we arrive at the following results: 



The general case - Theorem 13.61 iii 



d „ 1 

<Rr<^ 



h{F) ~ ~ h{F) 
• The diagonalizable case - Theorem 13.71 

Rr = . 

h{F) 
This completes the proof. ■ 

3.2 Generalizations and Applications 

In this section we generalize or apply our main results i.e. Theorems 13.61 and 13.71 in 
the following situations 

1. Coarse grained dynamics (Section 13.2.11) . 

2. Affine toral maps (Section 13. 2. 2j) . 

3. Degenerate noises f Section l3. 2. 3|) . 

4. Kinematic dynamo problem f Section 13. 2. 4|) . 
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3.2.1 Dissipation time of coarse-grained dynamics 

The uncertainties in the initial preparation and the final measurement of the noiseless 
system give rise to non-cumulative random perturbations to the system. The dynam- 
ics of such systems can be modeled by the coarse-grained family of noisy operators 
Te^ introduced in Section 0.1.4[ In this section we compute corresponding dissipation 
time fc- We show that, remarkably, for ergodic toral automorphisms both dissipation 
times are asymptotically the same, despite considerable difference in the structure of 
the corresponding noisy dynamical systems. 

To prove this result one can represent the action of Up or more generally Up in the 
Fourier series along the lines introduced in Section refOpt 

Tin ST^ (n) 

where m^J,, coincides with Uj^^' defined previously (cf. (j2.5(Jj) ) and 



''kM ~ 

0^fci,...,fc„_i( 



Uii, — 2_^ Uk,ki'^kiM---'^kn-i,k' 



which satisfies 



Y, l<lf = l, Vn, k. (3.37) 



Then 



7^(n)„ _ n TJ^n p, — n Tjn ~e\ko\'^ _ „-e|fcoP(^ Y^ „ (") „ 

Now we define 

5n(fc„) = {ko e zMo} : «£,,^ ^ 0}. 

Similar computations to these performed in Section 12.71 give the following general 
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upper bound for ||Ti,a || 

ll^?/f< E E l/>o)P E \utluf- (3.38) 

For a toral automorphism one easily sees that 

<!.„ = e-^(l'=°l"-^l-^"'»°l")4„.A^.o (3.39) 

and hence 

||T(")|| = p-^™i'lo5^feGZd(l'»o|" + |A'fcor) 
II e,a II 

The arithmetic minimization problem (|3.5p corresponding to the dissipation time of 



Te a now becomes 



min (|fc|° + U"fcr). (3.40) 

Oj^kezd 



The key observation is that, by the same arguments as before, similar estimates to 
these given in ()3.15p hold 

/ d \ " d 

^1 E \'"i^i'^i' < I'^r + i^"fcr < ^^2 E Aj"i^jfcr. (3.41) 

The remaining computations are the same verbatim so the dissipation time of T^^a 
and the family T}^ are equal asymptotically. 

3.2.2 AfRne toral maps 

In this section we present a slight generalization of the main results concerning the 
asymptotics of the dissipation time of toral automorphisms. Namely, we consider 
here general affine transformations of the torus. The term affine transformations 
will be used here to refer to homeomorphisms of the torus with zero periodic but 
not necessary zero constant part (cf. Section I3.1.1|) i.e. transformations of the form 
F = F + c, where F is a toral automorphism and c is a constant shift vector. 
We begin with a short discussion of the ergodicity of affine transforms. 
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The relation between ergodicity of a given affine transform F and associated with it 
toral automorphism F is summarized in the following proposition (for the proof we 
refer to appendix B) 

Proposition 3.18 Let F he any toral automorphism. Then 

i) If F is ergodic then F is also ergodic. 

a) If F is not ergodic then F is ergodic iff 1 is the only root of unity in the spectrum 

of F and c- k ^7J^ for any integer eigenvector k of F"^ . 

Proof, i) Assume F is ergodic and for some c, F = F + cis not ergodic. Then there 
exists non-constant / G Lq(T'^) satisfying f = f o F or in the Fourier representation 

Y, f{k)e^ = E e'^^^"^-7(A-ifc)efc (3.42) 

where A = F^ . Comparing the absolute values of the coefficients we get 

\f{k)\ = \hA--k)\ (3.43) 

for any integer n and any k. However, ergodicity of F implies that A~"fc ^ k for all 

fc 7^ 0, which contradicts our assumption that / G Lq(T'^). 

ii) We will use the following fact, which can be proved by simple application of 

rational canonical decomposition. For any A G SL{d,Zi) the following conditions are 

equivalent 

a) A possesses in its spectrum a root of unity not equal to one. 

b) There exists nonzero fc G Z'^ and a positive integer n such that k+Ak+.-.+A^-^^k = 
0. 

Now assume that 1 is the only root of unity in spectrum of F (and hence of A) and 
c • fc ^ Z"' for any integer eigenvector k of A, and that both F and F are not ergodic. 
The latter assumption implies the existence of a non-constant / G Lq(T'^) satisfying 
equations ()3.42|) and ()3.43p . Relation ()3.43|) clearly implies that if f{k) ^ then 
A"-k = k for some n. Moreover, since 1 is the only root of unity in spectrum of A, 
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we have, in view of b) that Ak = k. Thus the only possible non-constant invariant 
functions of F are single Fourier modes e^ corresponding to integer eigenvectors of 
A. But if such a Fourier mode is invariant under F then directly from (j3.42p one 
concludes that e^'^*'"'^ = 1 or equivalently k ■ c E Z'*, for some integer eigenvector of 
A. To prove the converse we assume that F is not ergodic and consider two cases: 
Case 1. A possesses in its spectrum a root of unity not equal to one. In this case 
according to condition b) there exists nonzero k E U^ and a positive integer n such 
that k + Ak + . .. + A^'^^k = 0, which implies in particular that A"-k = k and Ak ^ k. 
Now we define the function 

/ = efc + e e^fc + ... + e v^j-o ; e^n-i^ 

which clearly satisfies the condition f = f o F. This proves that F is not ergodic. 
Case 2. There exists integer eigenvector of A such that k ■ c E Z^. Then clearly for 
such k, f = Ck is -F-invariant and hence again F is not ergodic. ■ 

We recall that c = {ci,..,Cd) generates ergodic shift on the torus iff l,Ci,..,Cd are 
linearly independent over rationals. Thus as a direct consequence of the above propo- 
sition we get 

Corollary 3.19 If F is not ergodic and 1 is the only root of unity in the spectrum of 
F then F is ergodic for all vectors c generating ergodic shifts on the torus. 

Now we are in a position to state and prove the generalization of Theorem 13.61 from 
Section 13.1.21 to the case of affine transforms (the corresponding generalizations of 
Theorem 13.71 and Corollarv 13.81 are straightforward) 

Theorem 3.20 Let F be any affine transformation on the torus T'^, F associated 
with F toral automorphism and T^^t = G^^aUp. Then 
i) T^^a has simple dissipation time iff F is not ergodic. 
ii) T^^a has logarithmic dissipation time iff F is ergodic. 
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Hi) IfT^a has logarithmic dissipation time then the dissipation rate constant satisfies 
the following constraint 

<R,<^ 



h{F) ~ ~ h{F) 
where h{F) < h{F) is certain positive constant. 

Remark 3.21 The dissipation time of an affine transformation F is determined by 
ergodic properties of its linear part F and hence not by ergodic properties of F it- 
self. In particular all ergodic affine transformations associated with nonergodic toral 
automorphisms (cf. Proposition \3.1^} have simple dissipation time. 

Proof of Theorem 13.201 Specializing the general calculations of dissipation time 

presented in Section I^TTI to the case of affine transformations F = F + c, with nonzero 
c, one easily finds the following counterparts of formulas ()2.55|) and ()2.56|) 

Uk,k' - e OAk,k', 

Now, in order to determine the dissipation time of T^^a = G^aUp one has to determine 
the asymptotics of ||T"^|| as n goes to infinity. According to the above formulas and 
formulas ()2.54|) and ()3.4|) from Sections 12.71 and 13.1.31 the value of ||T"„|| does not 
depend on c, which reduces the proof to the case we already considered i.e. c = 0. 



3.2.3 Degenerate noise 

In this section we compute the dissipation time for non-strictly contracting gener- 
alizations of a-stable transition operators. Instead of considering standard a-stable 
kernels of the form ()2.15|) one can allow for some degree of degeneracy of noise in 
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chosen directions by introducing the following family of noise kernels 

g.,^,B{x):=Y.e-\^^^\''et.{x), (3.44) 

Where B denotes any dx d matrix with det 5 = 0. 

We denote by Ge,a,B the noise operator associated with ge,a,B- The degeneracy of 

B immediately implies that ||G'e,a,_B|| = 1 and hence the general considerations of 

sections 1 and 2 do not apply here. The answer to the question whether or not the 

dissipation time is finite depends on the choice of matrix B. 

For simplicity we concentrate on the case when B is diagonalizable. 

We call the eigenvector of B nondegenerate if it corresponds to nonzero eigenvalue. 

Theorem 3.22 Let F be any toral automorphism and T^,a,B = Ge,a,BUF- Assume 

that B is diagonalizable. Then 

i) If all nondegenerate eigenvectors of B* lie in one proper invariant subspace of F 

then dissipation does not take place i.e. Tc = oo. 

a) Otherwise the following statements hold. 

a) T^a h(is simple dissipation time iff F is not ergodic. 

b) T^^a has logarithmic dissipation time iff F is ergodic. 

c) IfT^^a,B has logarithmic dissipation time then the dissipation rate constant satisfies 
the following bounds 



h{F)~ h{Fy 

with some constant h{F) < h{F). The equality is achieved for all diagonalizable 
automorphisms F. 

Proof. 

We continue to use the convention A = F'^ . The general formula derived previously 
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for ||T"q,|| (see (13 .41) ). will now take the form 

n 

||T" .|| = sup e-^"^"=il^^''=l" =exp-e" inf V iBA'fcl". (3.45) 
Thus we need to estimate 

n 

inf V|5A'fc|". 

To this end we denote by /Xj (j = 1, ..., d) the eigenvalues of B and we construct a basis 
(of C^) {f j}j=i,...,d composed of normalized eigenvectors corresponding to eigenvalues 
\ij. We denote by Pja^^ ^ the set of eigen-projections on Vj, and by Lj the set of 
corresponding linear forms, given by the eigenvectors Uj of B\ which are of course 
co-orthogonal to i;^, i.e. (itj, Vj) = for i 7^ j. We have 

d d d 

x = J2 Pj^ = ^iLjx)vj = ^{x, Uj)vj, \/x e R'^. 

In subsequent computations the symbols Ci, C2 denote some absolute constants values 

of which are subject to change during calculations. 

We consider two cases. 

i) All nondegenerate eigenvectors of 5^ lie in one proper subspace of F. We have the 

following estimates 

d d 

i=i i=i 

d d 

2 



and 






d d 

d d 
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Since at least one of fij is zero and all nondegenerate vectors Uj lie in a proper 
invariant subspace of F, one easily sees that for each fixed n 

n n d 

inf V|5A'fcr= inf VV|/i,r|(fc,F''ii,)r = 0. 
i=\ i=\ jf=i 

ii) In this case we have the following upper bound 

n n n 

inf y|5A'fcr<||fir inf y U'fcr = ll^r min y lA'fcr. (3.46) 

In order to provide an appropriate lower bound we note that the set of vectors {F'^Uj}, 
where 1 < h < d and j runs through the indices of all nondegenerate eigenvectors 
of B, spans the whole space (otherwise all nondegenerate Uj would lie in one proper 
invariant subspace of F). We denote by {F'^'itj.} {1 < i < d) a. basis extracted from 
the above set. We can define now a new norm | ■ |u on M'^ by 

d 

and compute 

dn n—l d n~l d d 

1=1 1=0 h=l 1=0 h=l j=l 

n—l d d n—l d 

1=0 h=l j=l 1=0 i=l 

n—l d n—l 

1=0 i=l 1=0 

Using the equivalence between norms | ■ | and \ ■ \u and combining (|3.46|) with the 
above estimate we get 

n—l dn dn 

Ci min yU'^'fcr< inf V lEA'fcr < ||5r min V U'fcl". 
;=o 1=1 1=1 

This together with the obvious fact that /i(A'^) = dh{A) and the general estimate 
fl3.6|) reduces the proof back to the nondegenerate case considered in the previous 
section. ■ 
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3.2.4 Time scales in kinematic dynamo 

In this section we briefly discuss tlie connection between tlie dissipation time and 
some cliaracteristic time scales associated with kinematic dynamo, which concerns 
the generation of electromagnetic flelds by mechanical motion. For a general setup 
and discussion we refer the reader to IHU] and [Z21 and references therein. Here we 
restrict ourselves only to necessary deflnitions. 

Let B G Lq(T^, M.'^) denote periodic, zero mean and divergence free magnetic fleld and 
let F be the time-1 map associated with the fluid velocity. We deflne the push-forward 
map 

F,B{x) = dF{F-\x))B{F-\x)). 

The noisy push-forward map P^a on Lq(T°',]R'^) is then given by 

P,,«B := G^^^F^B, (3.47) 

where the convolution (the action of G^.a) is applied component-wise. 

It is said that the kinematic dynamo action (positive dynamo effect) occurs if the 

dynamo growth rate is positive i.e. 

1 



Moreover if 



Rdyn= lim -ln||P,"J| > 0. 



lim lim -In \\P"a\\ > 0, 

e— >0 n— ►oo fl ' 



then the dynamo action is said to be fast; otherwise it is slow. The anti-dynamo 
action takes place if 

lim -ln||P;„|| <0. 

Now we introduce the threshold time scale as 

nth = max{n : ||P,';^|| > e such that ||P"^i or ||P"+i < e}. 
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The threshold time nth{^) is of order 0(1) as e ^ for all fast kinematic dynamo 
systems. In the case of anti-dynamo action, nth{e) captures the longest time scale 
on which the generation of the magnetic field still takes place. Finally nth{^) is not 
defined for systems with do not exhibit any growth of magnetic field throughout 
the evolution. In the case of anti-dynamo we consider the time scale on which the 
generation of the magnetic field achieves its maximal value 



Up = min{n : ||-PJ^ctll — ^up ||-P( 



m I 
t,a\ 



which is called the peak time of the anti-dynamo action. 

Our next theorem establishes the relation between np, nth and Tc for toral automor- 
phisms. Thus dF = F and 



-l^ 



P,,aB = (7,,„ * F{B o F 



Theorem 3.23 Let F be any toral automorphism. Then 

i) If F is nonergodic and has positive entropy then for all < e < Rclnpp the fast 

dynamo action takes place with dynamo growth rate satisfying 

Rdyn = lnpF- eR~^ ^^ Inpp > 0, 

where pp denotes the spectral radius of F. The threshold time Uth is of order 0(1) 
and if F is diagonalizahle then nth ~ [Rdyn] ■ 

a) If F is nonergodic and has zero entropy then anti-dynamo action occurs and for 
nondiagonalizahle F , 

nth 1 

JXr — . 



In(nift) e 

Moreover there exists a constant < 7 < rf such that ||-Pe"a|| ~ (1/e)'^. If F is 
diagonalizahle then \\PlW\ is strictly decreasing (in n) and, hence, n^, = and nth is 
not defined. 
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Hi) If F is ergodic then anti-dynamo action occurs and Up ^ t^. In particular if F is 
diagonalizable then 

rip ^ nth- RJn{nth) ^ RJn{l/e) 

= -^ln(l/e) 
ah{F) 

and 



r;:sii~(iA)"'^(-). (3.48) 

We see that even in the case of anti-dynamo action the magnetic field can still grow 
to relatively large magnitude when the noise is small (power-law in 1/e). 
Proof. Representing the initial magnetic field B = (bi, ..., bd) in Fourier basis 

B= Y^ B{k)ek 

one obtains 

P,^B= Yl FB{k)e-^''\^^\''eAu. 



where we set A = (F ^^ . After n iterations 



Thus 



P:aB= Yl F-B{k)e-^"^'-^\^'''\"eA^^. 



KaBf < Y ii^"^(fc)re-'^"^"-^'^''^'" 

< max e-2^"Sr^il-^''=l" V \F"B{k)\' 
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where A;„ denotes a solution of the minimization problem 

n 

min VlA'fcr. 
The above calculation provides the following upper bound 

On the other hand let f„ denote a unit vector satisfying ||-F"|| = |-F"v|. One imme- 
diately sees that the above upper bound for ||-P"(^|| is achieved for magnetic field of 
the form B = v„efc^. Thus 

IIP" \\ _ ^-e"E?=i\A'k„\"\\pnu (OAO) 

ll-* €,a\\ ^ ll-' II- I^<J.'±^; 

Now we consider the cases mentioned in the statement of the theorem 
i) Nonergodic, nonzero entropy case. 
For any nonergodic map we have 



1=1 
This implies the following asymptotics 

II pn II ^ ~e"Rc^n\\ rpnu ^ (-e"Rc^+hipF)n+ci\nn+C2 (n ^(\\ 

\\^e,a\\ ^ ^ 11-^ II ~ "^ ; l^O.OUJ 

where Ci,C2 > are constants (both equal iff F is diagonalizable) . Thus for e < 
i?clnpir we have 

Rdyn = \npF- e°'R~^ ^^ Inp^ > 0. 

The threshold time is clearly of order 0(1) and in diagonalizable case can be written 

as 

1 e-.0 1 

In pf — e°'R^^ Inpp 

ii) Nonergodic, zero entropy case. 
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In this case Inpp = 0. Thus if F is nondiagonahzable then (j3.5(jp reads 

II p" II ~ p-<^°'Rc^n\i rpnu ^ -e°'Rc^n+ci\nn+C2 
ll-' e,all '^ C ll-T II ^^ e , 

with < Ci < d. This immediately yields 

D Ci nth 1 

e ln(nt/,) e 

And moreover ||-Pe"a|| ~ {1/eY'^. 

If F is diagonalizable then ||F"|| = 1 and in this case ||-P"a|| ~ e~'^°'^'' " which implies 

Up = 0. 

iii) If F is diagonalizable, then from ()3.6|) we know that for any < S < 1 and 

sufficiently large n 

n 

Thus for large n we have 

maxe +^pp < max||P"^|| < maxe "Vf- 

n n ' n 

We obtain the following constraints for Up 

1 ,JhP!L]^^,Jl)<„.^,Jh£^)^^,J' 



lnA+5 \lnA+5/ InA-^ \e/ lnA_A \lnA_A/ InA^^ Ve 

This gives 



lnA+5 e^o ln(l/e) InA.^' 

Now since \±s — > e"'*^^-' for 5^0 the above estimation yields the following asymp- 
totics 

^p ~ T7^ lii( l/e) ~ Tc, nth - Re In(nift) ^ Tc. 

h{F) 

Similar asymptotic estimates (except for the constant) hold for nondiagonahzable F. 
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3.3 Dissipation time of general C^ Anosov maps 

We recall that a diffeomorphism F : T'^ i-^ T^ is called Anosov if it is uniformly 
hyperbolic: there exist constants A > Q and < A^ < 1 < A„ such that at each 
£C G T'' the tangent space T^T'^ admits the direct sum decomposition TxT''' = E'^Q)E!^ 
into stable and unstable subspaces such that for every n &N, 






(d.f)(e:) = e^^, \\{DxF-)\es\\ < ax: 

{DxF){El) = El^, II (Z^.F-") 1^,11 < AK\ 

These inequalities have obvious consequences on the expansion rates of F and F~^, for 
instance they imply ||i5F||J^ > ||DF"||oo > A'^X^. As a consequence, the quantities 
of interest in Theorem 12.171 i) satisfy 

||DF|UA||DF-ioo > a„aa;\ 

lip A /iF-1 > A„ A A^"*^- 

All these expansion rates are > 1, so both noisy and coarse-grained dissipation times 
admit logarithmic lower bounds as in Eq. ()2.34|) . 

To obtain upper bounds, we use the mixing properties of this dynamics. Exponential 
mixing has been proved for Anosov diffeomorphisms of regularity C^^"^ (0 < r^ < 1) by 
Bowen ^Hl; using symbolic dynamics; the exponential decay is then valid for Holder 
observables in C^ for some < rj' < r]. 

Because we are interested in the noisy dynamics as well, we also refer to a more recent 
work [21j concerning C^ Anosov maps on T"', which bypasses symbolic dynamics. The 
authors construct an ad hoc invariant Banach space B of generalized functions on the 
phase space, such that the Perron-Frobenius operator is quasicompact on this space. 
The difficulty compared to the case of expanding maps, is that the space B explicitly 
depends on the (un)stable foliations generated by the map F on T^. Vaguely speaking, 
the elements of B are "smooth" along this unstable direction F^, but can be singular 
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("dual of smooth") along the stable foliation E^. The space B is the completion 
of C^lT'^) with respect to a weaker norm || ■ ||g adapted to these foliations. The 
definition of that norm is given in terms of a parameter < (3 < 1, the choice of 
which depends on the regularity of the unstable foliation. In general, the latter is 
only Holder continuous on the torus, with some exponent < r < 2 (note that the 
regularity of the foliations has little to do with the regularity of the map itself: a real- 
analytic map may have foliations being only Holder). Then, one must take (3 < t Al, 
and the authors prove that the essential spectrum of Pp on B has a radius smaller 
than ri = max(A~^, Af). This upper bound is sharper if /5 can be taken close to 1, that 
is, if the foliation is at least C^. This is the case for smooth area-preserving Anosov 
maps on T^, for which the foliations have regularity C^'^ for any 6 > jSH]- The 
operator Pp may have isolated eigenvalues 1 > |Aj| > ri, corresponding to eigenstates 
in B which are genuine distributions, outside Lq. As opposed to the radius ri, there is 
(to our knowledge) no simple general upper bound for largest resonance |Ai| in terms 
of (A„, Xs). 

The space C^(T'^) can be embedded continuously in both B and its dual B*, so that 
one can take s = s* = 1 in Eq. ()3.54|) . Therefore, for any 1 > a > max(|Ai|,ri), there 
is some constant C > such that for any f,gE Cq(T'^), 

Wn > 0, \Cj,,{n)\ < C \\g\\c^ \\f\\c^ a\ (3.52) 

In the proof of Theorem 12. 2(JI f Step 3), for the case s = s* = 1 we only need to assume 
that the noise generating function g is C^ with fast-decaying first derivatives. The 
fast decay implies that the first moment of g is finite (that is, one can take a > 1). 
The results of J2I] also concern the noisy dynamics. If the unstable foliation has 
regularity C^+^ with r^ > (for instance for any C^ Anosov diffeomorphism on the 2- 
dimensional torus), and assuming the noise generating function g G C^(]R^) with com- 
pact support, the authors prove the strong spectral stability of the Perron-Frobenius 
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operator Pp on the space B defined with the parameter (3 < rj. Therefore, the esti- 
mate (j3.52p also apphes to the noisy correlation function Cj^gi^) as long as e is small 
enough. 

Using these estimates and applying Theorem 12.171 (Section 12. 5|) and Corollary 12.211 
(Section 12. 6|) we obtain the following result regarding the dissipation time of C^ 
Anosov maps on the torus: 

Theorem 3.24 Let F be a volume preserving C^ Anosov diffeomorphisms on T'^ and 
let g be a C^ noise generating function with fast decay at infinity. 
I) Then there exist yU > A„ A A^"*^, < a < 1 and C > such that the dissipation time 
of the coarse-grained dynamics satisfies 

^ ln(e"i) -C <f,< ^-^ ln(e-^) + C, 



In // I 111 cr I 

II) If in addition F has C^'^'^ -regular foliations and g G C^(R'^) is compactly supported, 
then the dissipation time of the noisy evolution satisfies for some C > and small 
enough e: 



ln(e-^) -C <T,< —-— \n{e-') + C 



InllDFlU - I In 



CXD 



O" 



3.4 Examples and general comments on other chaotic sys- 
tems 

3.4.1 Noiseless correlations 

Let us first consider the correlation functions Cf^g{n) for the noiseless dynamics. A 
common route to prove that a map F is mixing consists in studying the Perron- 
Frobenius operator Pp on some cleverly selected space B of densities (this may be a 
Banach or Frechet space). The objective is to prove that the spectrum of Pp on B is of 
quasicompact type. More precisely, one shows that Pp admits 1 as simple eigenvalue 
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(with the constant eigenfunction), and that the radius of its essential spectrum is 
smaller than some < ri < 1. The spectrum outside the disk of radius ri is made of 
a finite number of isolated, finitely-degenerate eigenvalues {Aj} of moduli ri < |Aj| < 1 
(usually called Ruelle-Pollicott resonances). We order these eigenvalues according to 
their decreasing moduli, so that the largest modulus is |Ai|. This spectral structure 
of Pp on B implies that for any 1 > cr > max(|Ai|,ri) (if the largest resonance Ai is 
semisimple, one may take a = |Ai|), there is some constant C = C(cr) such that for 
any f e Bq, g e B*: 

< II^IU*l|i^;iUoll/IU 

< C\\g\\s^\\fhcr": (3.53) 

The choice of the space B is crucial here: in general the spectrum of Pp on Lq 
intersects the unit circle, so there is no chance to prove exponential decay within the 
L^ framework. B can be a Banach space of bounded variation. Holder or C^ functions 
which embeds continuously in L^ (this is the choice made for F uniformly expanding, 
see Section I3.4.3|) : it may also be a space of generalized functions lying outside L^ 
(case of Anosov maps, see Section IX5|1 . 

In general, there exist Holder exponents < s* < s such that C^ (resp. C^*) embeds 
continuously in B (resp. its dual B*). As a result, the upper bound (|3.53|) induces 
that for any f & Cq, g E Cq* and any n > 0, 

\Cf,gin)\<C\\g\\cs,\\f\\csa\ (3.54) 

This is the form of upper bound we used in Theorem 12.201 When the radius a is 
given by a resonance Ai, it can often not be decreased when one takes observables 
of higher regularity. In this case, it is preferable to take in the above estimate the 
weakest norms possible, that is take s and s* as small as possible (to obtain better 
upper bounds in Corollarv l2.2H) . 
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This strategy of proof has been apphed to several types of maps [121 [13] , including 
the (noninvertible) expanding maps and the Anosov or Axiom-A diffeomorphisms 
on a compact manifold. We will give some details and examples of these two types 
of maps in the next subsections. Exponential decay of correlations has also been 
proved (using various methods) for piecewise expanding maps on the interval, some 
nonuniformly hyperbolic/expanding maps, like "good" logistic maps on the interval, 
or "good" Henon maps; some expanding or hyperbolic maps with singularities. In 
those cases, the decay rate may have no (obvious) spectral interpretation, as opposed 
to the formalism described above fH^ . 

Different types of decay, like the stretched exponential Cf^g{n) < Kr^ with < r < 1 
and < .^ < 1 have been proved for some Poincare maps of hyperbolic flows and some 
random nonuniformly hyperbolic systems; yet, it seems that the bound may not be 
optimal in some of these cases, but rather due to the method of the proof. On the 
opposite, the polynomial decay Cf^giji) < n^^ was shown to be optimal for some 
"intermittent" systems, like a one-dimensional map expanding everywhere except 
at a fixed "neutral" point (such maps are sometimes called "almost expanding" or 
"almost hyperbolic" ) . 

Remark: In general, an expanding or Anosov map does not preserve the Lebesgue 
measure, so the first step before dealing with correlations is to precise the invariant 
measure with respect to which one wants to study the ergodic properties. In the 
"nice" cases, one can prove the existence and uniqueness of a "physical measure", 
which is ergodic for the map F, and then study the correlation functions with respect 
to this measure. The formalism of Perron-Frobenius operators applies also to this 
general case, the physical measure being related to the simple eigenstate of Pp (on 
the space B) associated with the eigenvalue 1. As stated in the Introduction, in this 
work we only consider maps for which the physical measure is the Lebesgue measure. 
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3.4.2 Noisy correlations 

There exist fewer results on the decay of correlations for stochastic perturbations of 
deterministic maps, like our noisy evolution T^. In general, one wants to prove some 
sort of strong stochastic stability, that is stability of the invariant measure, and of the 
rate of decay of the correlations in the limit when the noise parameter e vanishes. 
In the case of exponential mixing, this strong stochastic stability implies that for 
small enough e, the upper bound (j3.54|) is "stable" when switching on the noise: for 
small enough e > there exists a radius cTe-^-^o- such that for any f E C^, g E Cq* , 

Wn > 0, \Clg{n)\ < C \\g\\cs, \\f\\c^ a,". (3.55) 

This stability has been proved for smooth uniformly expanding maps |14j, some 
nonuniformly expanding or piecewise expanding maps (see the review in [^ or the 
book P2])- It has been shown also for uniformly hyperbolic (Anosov) maps on the 
2-dimensional torus [21]. The proof uses perturbation theory: one shows that the 
isolated eigenvalues of G^ o Pp on B of moduli |Aje| > ri (and the associated eigen- 
states) vary continuously w.r.to e when e -^ 0. Therefore, one can choose a rate 
cr^ > max(|Ai^e|)''"i) which is a continuous function of e. 

In the next two sections, we describe in some detail the exponential decay of cor- 
relations for smooth uniformly expanding maps and Anosov diffeomorphisms on the 
torus. 

3.4.3 Smooth uniformly expanding maps 

Let F be a C'^"^^ map on T'^ (with s > 0). Assume that there exists A > 1 such that 
for any a; G T^ and any v in the tangent space T^T"^, \\DF{x)v\\ > X\\v\\ (we assume 
that A is the largest such constant). Such a map is called uniformly expanding. In 
general, it admits a unique absolutely continuous invariant probability measure; here 
we will restrict ourselves to maps for which this measure is the Lebesgue measure. 
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Ruelle |112j proved that the Perron-Frobenius operator Pp of such a map is quasi- 
compact on the space C"'(T'^), and that its essential spectrum is contained inside the 
disk of radius ri = A^'^. In general, one has little information on the possible discrete 
spectrum outside this disk (upper bounds on the decay rate have been obtained in 
the case of an expanding map of regularity C^+'' [12]). Strong stochastic stability for 
such maps was proved in ^^, with a more general definition of the noise than the 
one we gave. 

For all these cases, one can take s^, = 0, since the continuous functions are continu- 
ously embedded in any space (C"*)*. 

Case of a linear expanding map 

We describe the simplest example possible for such a map, namely the angle- doubling 
map on T^ defined as F{x) = 2x mod 1. This map is real analytic, with uniform 
expansion rate A = 2. Due to its linearity, the dynamics of this map (as well as its 
noisy version) is simple to express in the basis of Fourier modes ek{x) = e^*'^'^^: 

Wk G Z, Upek = e2k 

=^ T.efc = g{ek)e2k 

n 

^T^Bk = [I[g{e2^k)]e2nk 
i=i 

The computation is even simpler for the coarse-grained propagator: 

Ti")efe = g{k)g{2^k)e2r^k. 

To fix the ideas, we take for the noise generating function g{^) = e~'^'° for some 
< a < 2. One easily checks that for any n > 1, 

nna i 

'Till IIT^"^ II r^^r-r^ ) ^Q: 



T" = TVJ =exp<^ -e 

II e II II E ill i- I ^ _2-a 

|f;(")|| = exp{-e"(2"" + l)}. 
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For any e > 0, these decays are super-exponential: the spectrum of T^ on Lg is reduced 
to {0} for any e > (the spectrum of Up is the full unit disk). From there, we get 
explicit expressions for both dissipation times: 

rc = :^ He-') + 0{1), f. = ji^ ln(e-^) + 0{1). (3.56) 

For this linear map, 2 = ||Z}F||oo = fJ'F, so this estimate is in agreement with the 
lower bounds ()2.34|) . the latter being sharp if a G [1,2]. On the other hand. In 2 is 
also equal with the Kolmogorov-Sinai (K-S) entropy h{F) of F. Therefore, for this 
linear map the dissipation rate constant exists, and is equal to ^7^. 
To compare these exact asymptotics with the upper bounds of Corollary 12.211 we 
estimate the correlation functions Cf^g{n) on the spaces C"'(T^). We give below a 
short proof in the case s > | . We will use the following Fourier estimates |131j : 



3C>0, ,V/eCo^(T^), Vfc^O, \f{k)\<C 



c 



Therefore, writing the correlation function as a Fourier series, we get: 

^WPm < C'^^. (3.57) 

This estimate yields a decay of the correlation function as in Eq. ()3.54|) . with a rate 
a = 2~* and s* = 0. One can check that this rate is sharp for functions in C^: 
indeed, any z E C, \z\ < 2^* is an eigenvalue of Pp on that space. Applying the 
Corollarv 12.211 / ii), we get that for any s > 0, there exists a constant c such that 

rc<^ln(e-i) + c (3.58) 

sm2 

for sufficiently small e. The exact dissipation rate constant l/ln2 is recovered only 
for large s. 

A straightforward computation shows that the estimate ()3.57|) also holds if one re- 
places Pp by PpoGe', hence the noisy correlation function dynamics satisfies the same 
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uniform upper bound as the noiseless one, with the decay rate o^ = 2^*. As a result, 
the upper bound on Tc given by Corollary 12. 211 /7 ii) is the same as in Eq. (|3.58|) . 

3.5 Technical proofs 

Proof of Proposition I3.13L 

For the purposes of the proof we use the following abbreviation 

• PRS{M.'^) - proper rational subspace of M'^. 

• PIS{A, M.'^) - proper A-invariant subspace of M*^. 

• PRIS{A, M'^) - proper rational A-invariant subspace of M'^. 

a) ^ b) . Suppose there exists PRIS{A,W^) Si. Let Ai be a matrix representing 

invariant rational block associated with 5*1. Then Ai is rational matrix and its 
characteristic polynomial Pi belongs to Q[x]. Let P denote the characteristic 
polynomial of A. Then P = P1P2 and since both P, Pi G Q[x] then also P2 G 
Q[x], which means P and hence A is not irreducible. 

b) ^ c) . Assume there exists PRSiR"^) S contained in PIS{A,R'^) V. Take any 

rational vector q G S and let do = dimV then the set {q,Aq, ..,A'^°^^q} spans 
PRIS{A,R^). 

c) ^ d) . Assume that for given q and an arithmetic sequence ni,...,nd, the set 

S = {A"'^q,A"''^q, ...,A'^'^q} does not form a basis. Since for some fixed integer 
r, ni = rii + {I — l)r, we have A'^'-q = {A^'Y'^ A'^'^ q = {A^Y'^q, where q = A^^q. 
Now consider the biggest subset 5*0 = {q,A^q,{A^Yq,...,{A^Y^^^q} such that 
do < d and 5*0 is linearly independent. Obviously 5*0 spans a PRISMA') which 
is also a PRIS{A). 



Section 3.5. Technical proofs 103 

d) ^ a) . Suppose that characteristic polynomial P of A is not irreducible in Q[x]. 

Then P = P1P2, with Pi, P2 ^ Q[x]. From the Cayley-Hamilton theorem we get 
that = P{A) = Pi{A)P2{A). Hence for any nonzero rational vector q, either 1) 
P2{A)q = or 2) qr := P2(A)qr ^ and Pi(A)qr = 0. Since max{rfe^(Pi, P2)} < 
d, there exists a nonzero rational vector q (namely q or q) such that the set of 
iterations {q, Aq, A?q, ..., A'^~^q} does not form a basis of M'^. 

e) ^ f ) . Assume there exist nonzero q E Q'^ orthogonal to certain PIS{A, W^) V. 

Then for any n and any f E V, {{A^Yq.f) = {q,A'^f) = and hence S = 
{q, A'^q, (y4"'')^q, ..., {A'^Y^^q}, cannot form a basis, which in view of equivalence 
a) -v^ d) implies reducibility of A^. 

f ) ^ g) . Suppose there exists PIS {A, M'^) V contained in certain PRS(M.'^) S. 

Since 5* is rational, S-^ is also rational. Consider any rational vector q G S*-*", 
then (q, /) = for any f EV. 

g) ^ b) . If there exists PRIS{M.'^), then this subspace is A-invariant and contained 

in PRS(R'^) i.e in itself. 

Now since b) is equivalent to a) it is enough to establish the equivalence between a) 
and e) to complete the proof. But the latter equivalence is obvious in view of the fact 
that A and A'' have the same characteristic polynomial. ■ 

Proof of Proposition 13. 4L 

Suppose A is a toral automorphism of zero entropy. The latter property is equivalent 
to the fact that all the eigenvalues of A are of modulus 1. Let Pa be a characteristic 
polynomial of A. Consider any irreducible over Z factor P of polynomial Pa and 
construct a toral automorphism B such that its characteristic polynomial is equal 
to P . Obviously all the eigenvalues of B are also the eigenvalues of A, and each 
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eigenvalue of A can be found among eigenvalues of some matrix B of this type. 

Irreducibility of P implies irreducibility and hence diagonalizability of B. 

Thus for any nonzero vector fc G Z*^ and any positive integer n the following estimate 

holds |i?"k| < |k|, which implies the existence (for each k) of some integer r such 

that B^'k = k. 

The latter shows that all the eigenvalues of B (and hence also of A) are roots of unity. 



Proof of Proposition I3.15L 

We first show that irreducible polynomials P G Q[x] do not have repeated roots. 
Indeed suppose A is a root of P of multiplicity greater that 1, then A is also a root of a 
derivative polynomial P' G Q[x]. Since the minimal polynomial of A must divide both 
P and P' and deg{P') < deg{P) one immediately concludes that P is not irreducible. 
Now, suppose A G GL{d, Q) is completely decomposable over Q and let ()3.9|) be its 
block diagonal decomposition into irreducible blocks. Each Pa , as a characteristic 
polynomial of Aj, is irreducible over Q and hence does not possesses repeated roots, 
which implies diagonalizability of each Aj and hence of A. On the other hand if 
A is diagonalizable then its minimal polynomial does not possesses repeated roots, 
which implies that all characteristic polynomials associated with elementary divisors 
are (first powers of) irreducible polynomials. This implies irreducibility of each block 
in representation ()3.9|) . ■ 

Proof of Proposition I3.16L 

Let Pa be the characteristic polynomial of an irreducible matrix A G GL{d, Q). Since 
Pa is an irreducible element of Q[x] it does not possesses repeated roots (see the proof 
of Proposition I3.15|) . ■ 
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Proof of Proposition 13.91 

Combining formula 13.41 and Tlieorem 13.101 in case of a-stable noises we get tliat for 
any 6 > 

\\T:J < e-^"''"-"'"''"'" . (3.59) 

Using this estimate one immediately gets 

CUn) = {f.T-^h) < 11/11 ||/.i|||T;;j| < 11/11 ||/^||e-^"^"". 

Now let / = G*e/o and h = G^Qq. Since the estimate (j3.59|) holds also in coarse-grained 
version, we have for a = 2 

Cf,hin) = (/,f/;/i) = (Gjo,t/]^GeM = (/o,Ti")M 

< ||/o||||/^o||||Ti")||<||/o||||/io||e-^'"". ■ 
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Chapter 4 



Interludium 



This chapter - different in its form from the rest of the work - is meant as an intro- 
duction and historical overview of the development of the idea of quantum mechanics 
on the torus. We put an emphasis on the semiclassical analysis and we select only 
most important contributions to the subject directly pertaining to the problems con- 
sidered in this work. Hence, although a lot of other important topics are going to 
be omitted, this introduction should be helpful for someone who wants to enter the 
field or at least understand some of its results without prior familiarity with it and 
without an intention to study all subsidiary contributions and developments which 
have been introduced during the now almost 25-years long history of the subject. 

In the late seventies and early eighties both mathematicians and physicists re- 
alized the importance of and started to discuss seriously the problem of chaoticity in 
quantum dynamical systems jSH 1127] . 

By that time the analysis of chaotic behavior of classical discrete-time dynam- 
ical systems on compact phase spaces (represented usually by a torus) constituted 
one of the most developed and best understood branches of the general theory of 
dynamical systems fOl IHS EHl- In fact, even today, most of the well-known and 
deeply understood examples of fully chaotic systems are described in terms of toral 
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maps (ergodic toral automorphisms, Baker and sawtooth maps, angle doubhng maps, 
general Bernoulli shifts etc.). 

Despite some controversy concerning physical interpretation (periodicity in mo- 
mentum), from a mathematical point of view it became clear that chaotic toral maps 
should provide a convenient testing ground for a newly born field of study, now usu- 
ally referred to as 'quantum chaos'. The quantization on the torus became at that 
point a necessity. As in the general case of quantization of any dynamical system 
the procedure is highly non-unique and from the very beginning different approaches 
had been taken and developed. Here we consider the two most popular ones, de- 
scribed respectively in finite and infinite dimensional settings. We start with the 
finite dimensional case. 

In 1980 M.V. Berry and J.H. Hannay [U^ made a first breakthrough in the area 
by introducing the notion of the finite dimensional quantum Hilbert space associated 
with a classical symplectic 2-torus. The dimension of the space in their model de- 
pends on the value of the Planck constant which itself is restricted to reciprocals of 
integers. The later condition results from the assumption, made by the authors, that 
quantum mechanical wave functions of any system on the torus should be periodic 
in both position and momentum representations. This assumption restricts the set 
of admissible pure states of the systems to periodic Dirac delta combs with equally 
spaced spikes. The distance between two neighboring delta functions in the comb is 
equal to the Planck constant (the smallest possible quantum resolution of the phase 
space) and since the whole comb is 1-periodic ('lives' on the torus) the constant can 
assume only inverse integer values h = 1/N. Each delta comb can be identified with 
a set of its A^ complex amplitudes, i.e., the strengths associated to its 'delta spikes'. 
The corresponding quantum Hilbert space of all pure states is then A^-dimensional 
and can be identified with C^. 

This way the kinematic step of the quantization (i.e., description of the quantum 
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phase space) had been completed. The next step was to define the quantum dynamics 
on that space. Hannay and Berry chose to work with the well-known from the classical 
setting ^U] discrete-time cat map dynamics. The assumption of the strict periodicity 
of their quantum phase space restricted in a considerable way the class of quantizable 
maps. The authors came to the conclusion that in order to be quantizable a cat map 
had to satisfy a so-called 'checkerboard condition', i.e., its matrix had to assume one 
of the following forms [1°] or [eo]? where 'e' and 'o' denote respectively even and 
odd integer entries. 

Under this assumption the authors constructed the quantum propagator (A^ x A^ 
unitary matrix) associated with a given map. One of the key points of choosing the 
cat map dynamics for quantization lay in the fact that the Hamiltonian of the clas- 
sical map is in this case quadratic (see Section lA.ll in Appendix Ej) and hence the 
semiclassical approximation to the Green function of the corresponding Schrodinger 
equation on M^ is exact [89j. This allowed the authors to define the quantum prop- 
agator on the torus as a 'projection' of the standard propagator on the plane. The 
projection, which in this case amounts to the discretization (applying the propagator 
to a delta comb) and periodization (wrapping the result around the torus) was per- 
formed formally by averaging of the kernel of the planar quantum propagator over 
all integer winding numbers associated with a classical path on the tours. This way 
the authors constructed a discrete and finite kernel (A^ x A^ matrix) representing the 
unitary quantum propagator of a toral map. The entries of the matrix were expressed 
in terms of oscillatory Gauss sums and the method introduced considerable number- 
theoretical complications as well as some difficult to resolve normalization issues. It 
thus became clear that except for some trivial cases the explicit formula for the ma- 
trix of the propagator was unattainable. Nevertheless, the authors were able to prove 
periodicity of such a quantum propagator, i.e., the existence of an integer P{N) such 
that f/-^(^) = e**^*-^-*!. This result yielded the following important information re- 
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garding the spectrum of U: each eigenvalue e*"-' of U is constrained to take one of 
the possible P{N) values specified in terms of their eigenangles 

^ P{N) ^ ' 

Not all of the values have to be assumed however, and on the other hand, some eigen- 
values may be (and in fact usually are) highly degenerate. Thus in order to determine 
the spectrum one needs to find the period P(N) and the multiplicities dj associated 
with each Uj. Hannay and Berry conjectured (and supported the conjecture with 
numerical evidence) that for regular maps P{N) would be either bounded (elliptic 
case) or grow linearly with A^ (parabolic case) and that it would grow 'in average' as 
A^ but erratically in the chaotic (i.e., hyperbolic) case. 

This line of research had been continued by Tabor |12()j and further developed 
in 1991 by J. P. Keating |HH1 IH^ . Using probabilistic number-theoretical approach 
supported by numerical computations Keating argued jHE] that for chaotic cat maps 
P{N) grows in fact slightly sublinearly (see below) and indeed highly erratically - the 
fluctuations themselves being of the order of A^. More precisely, Keating studied the 
asymptotic behavior of the average order of P{N) denoted by {P{N)) and defined as 

1 ^ 

n=l 

The reason to study {P{N)) instead of P{N) is that as opposed to P{N) the cumulant 
function C{N) := Yln=i -P(^) behaves very regularly and grows in A^ 'almost' like A^^. 
The normalized cumulant, i.e., the averaged order {P{N)) is then particularly suitable 
for asymptotic analysis. The result obtained by Keating states that for arbitrary small 
6 > and arbitrary big p > 0, 

N'-'<{PiN))<i\nNy^N. (4.2) 

From this Keating concluded that in average P{N)/N tends in fact to zero, although 
at a very slow rate. 
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As to the precise behavior of P{N) itself (for large A^) not much can be said 
because of its chaotic dependence on A^. The size of the fluctuations of P{N) had 
been estimated analytically by Kurlberg and Rudnick in [77| [TH] . The result states 
that there exist constants c and C such that for every A^ 

cln(A^) < P{N) < CA^ln(ln(A^)). (4.3) 

It is worth pointing out that the bounds are sharp. In particular it has been proved 
recently in [H3] that for every hyperbolic cat map there exists a sequence {N^} such 
that P{Nk) ~ ln(A"fc). As will be described later in more details the existence of such 
sequences plays a crucial role in the semiclassical analysis of ergodic properties of 
quantum cat maps. 

Once some knowledge about the properties of the quantum period function had 
been accumulated the next step was to study the degeneracies dj of each eigenvalue 
to obtain the information about the spectrum of U. 

It is easy to show (see ^^) that dj can be represented in terms of the trace of 
the powers of the propagator U: 

P{N) 

d^ = plN)J2^'^U'^^~"''- (4.4) 

In jnni Keating estimated the value of dj for cat maps by deriving an exact version 
of the Gutzwiller's formula 



Tr(f/'=) = -L V exp ('-^{c'^m^ - h^n" + 2{d^ - l)mn)\ , (4.5) 



m,n£'P 
' a b] D. _ ^of l' T?k _ T\ _ o ^ "Ty T?k 



where F = ["^], Rk = det(F — 1) = 2 — TrF and V denotes the fundamental 
parallelogram formed by the action of the planar linear map F^ — I on the unit 
square. This parallelogram encompasses all the representatives of the periodic orbits 
of F of period k (all nonequivalent fixed points of F''). 

Thus just like in the case of the general Gutzwiller's formula, ()4.5p relates the 
trace of U^ and hence in view of ()4.4|) and ()4.1|) the distribution of the eigenvalues 
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(the spectrum) of the quantum propagator U to its average over periodic orbits of 
the classical map (i.e., over the fixed points of F^). It has however two important 
features which distinguish it from the general Gutzwiller's formula. First of all, the 
summation does not constitute merely a semiclassical approximation but gives the 
exact value of the trace for all (allowed on the torus) values of Planck constant 
h = 1/N. Secondly, the formula is explicit. Thus in this particular case the most 
important feature of Gutzwiller's formula, i.e., the fact that it translates 'hard-to- 
compute' iterations of the quantum propagator into 'easy-to-compute' powers of the 
classical map, can be used to compute exactly the spectrum of U. For a fixed and 
small value of A^ this can be done numerically. For large N (which is important for 
semiclassical analysis) the number of periodic orbits proliferates exponentially and 
the numerical approach becomes impractical. But the formula remains useful for 
asymptotic analysis especially since the behavior of the periodic orbits of classical cat 
maps is now well understood (for further details in this direction we refer to |l(J5j and 

m)- 

Having described the first stages of the development of the finite dimension 
quantization on the torus we now compare it with another approach - the infinite 
dimensional one, developed almost simultaneously but by different authors and for 
different reasons. 

Just like the spectral approach to quantum chaoticity provided the main mo- 
tivation for the introduction of the finite dimensional quantization on the torus, the 
infinite dimensional approach was initiated during the attempts to extend the notion 
of the entropy and K-property from classical ergodic theory to the quantum setting or 
more generally to the case of dynamical systems defined on noncommutative algebras. 

Below we briefly discuss the role which the study of quantum entropy played in 
the development of the infinite dimensional noncommutative mechanics on the torus. 
We want to stress that the discussion is not meant to summarize the whole (much 
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longer) history of the notion of quantum entropy but to mention the developments 
relevant to quantization on the torus. In particular we will not have space here to 
discuss (but instead take for granted) von Neumann's definition of the entropy of 
the density matrix ^T7\ IHH] , its extension to normal states on von Neumann algebras 
(obtained after the work by Araki 9j), Lieb's results on strong subadditivity and 
WYD conjecture [501 IM| etc. For a comprehensive overview of these notions and 
results we refer to |lU2j . 



In 1975 A. Connes and E. St0rmer ^^ introduced the notion of the entropy 
for automorphisms on certain noncommutative von Neumann algebras, namely on a 
hyperfinite type J/i factor (see Definition 15. 6|) . These algebras occur in the context 
of quantization on the torus. As was suggested by the authors, the main difficulty in 
generalizing the notion of the classical KS-entropy to non abelian setting seemed to 
lay in the fact that two finite dimensional subalgebras of the nonabelian von Neumann 
algebra need not generate a finite dimensional algebra (for a simple example see j^lj ) , 
and hence there was no immediate analogue of the classical operation of the refinement 
Py Q oi two finite partitions P, Q of a phase space - a key notion used in the classical 
setting. We will see later in the discussion that an appropriate analogue in fact can be 
constructed by means of operational partition of unity but it leads to a different notion 
of quantum entropy [so-called ALF-entropy] . The solution the authors provided at 
that time was the following: instead of replacing the notion of the refinement of 
two partitions by its quantum analogue, the notion of its entropy hxsiP ^ Q) was 
replaced by a nonabelian counterpart H{P,Q). The way the replacement can be 
made is, as one can expect, non- unique and different choices led to nonequivalent 
notions of quantum entropy. An interesting feature is that different but classically 
equivalent approaches can still yield nonequivalent quantum counterparts. Indeed, 
Connes and St0rmer constructed a replacement for hKs{P V Q) based on its original 
definition, while independently but at almost the same time Emch jlTj introduced 
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another replacement, based on a (classically equivalent) approach via conditional 
entropy. The resulting quantum notions did not coincide. In the original work by 
Connes and St0rmer the noncommutativity of the algebra was 'overshadowed' to 
some extent by the requirement that the invariant state of the dynamics be tracial. 
In this case the state itself does not see the noncommutativity of the underlying 
algebraic structure. The generalization to nontracial states and to arbitrary C*- and 
W* algebras was achieved by Connes, Narnhofer and Thirring in 1987 [3^ and resulted 
in the now well-known CNT entropy (in 1992 Sauvageot and Thouvenot introduced 
an alternative definition |116j and showed that for hyperfinite algebras all three, i.e., 
CS, CNT and 'ST' entropies coincide). 

The main idea of the CNT construction is to introduce an abelian model of 
the noncommutative algebra and to transfer via a completely positive unital (CPU) 
map the invariant state and the dynamics from the nonabelian to the abelian setting. 
After appropriate corrections (taking into account the 'entropy defect' introduced by 
the CPU map) the entropy is computed in a classical manner. 

After the definitions had been established the authors of the newly introduced 
notions started to look for examples of dynamical systems on which the notions could 
be tested and compared. As in the previous (finite dimensional) case the well known 
examples of classical chaotic toral maps once again came to the attention. At that 
time, however, it seemed obvious that an infinite dimensional model is absolutely 
necessary [that it is not exactly the case will become clear later in the discussion]. 
The main argument was that, just like in the classical case, in principle, the system 
can have nonzero entropy only if the generator (Koopman operator) of the dynamics 
possesses a continuous spectrum (cf. [3^]). The latter requirement obviously calls for 
an infinite dimensional algebra of observables [this was the main reason why CS and 
CNT entropies were introduced in infinite dimensional settings]. 

Thus in 1991 Benatti, Narnhofer and Sewell ^H] introduced a non-commutative 
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version of the Arnold cat map. Unlike in the case of the finite dimensional quantiza- 
tion, the authors didn't pay attention to the underlying physical space of pure states 
of the system and instead began the construction by introducing the appropriate 
noncomutative algebra of (quantum) observables. The most natural choice was the 
discrete (countably generated) Weyl-Heisenberg algebra spanned by the elements W^ 
(Weyl translations) indexed over the integral lattice Z^'^ and satisfying the standard 
canonical commutation relations (CCR) 

WkWrr, = e-^^'»^"^iyfc+^, (4.6) 

WkWm = e^-^'^'^^^^iy^l^fc. (4.7) 

[The detailed technical description including the explanation of the seemingly non- 
standard placement for the Planck constant will be given in the next Chapter]. 

After introducing the cat map dynamics and an appropriate tracial state corre- 
sponding to the Lebesgue measure of a classical system the authors constructed the 
GNS representation, and hence in particular the Hilbert space ((/^(Z^))) of states for 
the system. It turned out that the unitary map implementing the quantum auto- 
morphism on this space coincided with the Koopman operator of the classical map. 
Since the result turned out to be independent of the value of the Planck constant it 
became clear that the model would not be suitable for semiclassical considerations. 

Nevertheless, the model still seemed adequate for the computations of CS-CNT 
entropies. Indeed, the above mentioned algebra is hyperfinite and either reduces to a 
tensor product of L°°{T^) and a finite matrix algebra (if h is rational) or is irreducible 
and of type J/i in irrational case (cf. Definition 15. 6|) . In either case both CS and 
CNT entropies are well defined and coincide. In the rational case the CS entropy of 
a quantum cat map was shown jT^I to coincide with a classical one. But the result 
did not contribute any new information as far as the quantum setting is concerned. 
Indeed, the cat map dynamics factorizes into a purely classical part acting on L°°(T^) 
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with its classical entropy and a finite dimensional quantum remainder with trivially 
vanishing quantum entropy. 

In hope to obtain some nontrivial information one needed to look at 'irrational- /i' 
case. It soon became clear, however, that in this case exceptionally strong properties 
of clustering and asymptotic abehanness jHSl HH IHE] were necessary to secure the 
possibility of non-vanishing CNT entropy. Without entering into the details of the 
definitions, let us only recall that the notion of strong asymptotic abelianness was 
introduced in [THj and that it was basically shown in [8^ that its lack implies zero CNT 
entropy. It was then found [^21 that indeed, except for a possibly measure zero set of 
the values of the Planck constants, the CNT-entropy of the cat map is in fact zero. 
On top of that, in the course of the investigations, it turned out [95j that using the 
above mentioned property one can show that CNT-entropy fails to be additive w.r.t. 
the tensor products. Thus although CS-CNT entropies had been and are successfully 
applied to different dynamical systems another notion of noncommutative entropy 
was still needed. 

In 1994 R. Alicki and M. Fannes [H] introduced a new approach to quantum 
dynamical entropy. Finite partitions used in the definition of the classical KS entropy 
had been replaced by the authors by (also finite) operational partitions of unity in 
corresponding quantum algebra. Such partitions could be evolved with a dynamics 
and composed among themselves yielding at each step finer but always finite new 
operational partitions. Using the standard notion of the von Neumann entropy of 
a quantum state and applying the procedure similar to the one known from the 
classical setting the authors constructed a new quantum dynamical entropy, now 
usually referred to as ALF entropy (in recognition of the connection of their idea to 
an earlier work by G. Lindblad |82j). 

The definition was clearly compatible with a classical one in abelian case. More- 
over it was shown |3] that, e.g., in the case of toral automorphisms in arbitrary di- 
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mensions the computations of the classical KS entropy using the operational approach 
could be significantly simplified in comparison to the traditional approach \V2b\ . The 
entropy has also been successfully computed for quantized cat maps in the infinite 
dimensional setting ,3; with the result identical to the classical one for all values of 
Planck constant (the result is consistent with the fact that the generator of the dy- 
namics is independent of its value). The latter example indicated that the notion 
differs significantly from CS-CNT entropies. In fact there even exist examples of 
systems for which Half = oo while Hqnt = (see 0)- The two entropies differ 
also in the cases where the latter one is nonzero (e.g., for the shifts on quantum spin 
chains JHIIH])- Nevertheless, similarly as all previously mentioned quantum dynam- 
ical entropies, ALF entropy also vanishes for any fixed finite-dimensional quantum 
dynamical system - the phenomena known as saturation. The word saturation is a 
key notion here and it is worth to pause for a while and explain the source of the 
phenomenon in more detail. 

In the classical setting KS entropy is defined as an infinite-time limit of the 
appropriately rescaled entropy production which occurs during the evolution of the 
system. Since the classical phase space can be refined into arbitrary small partitions, 
the rescaled entropy production converges to a well-defined limit and that limit is 
understood as the entropy of the system (in some, but not most interesting, cases the 
limit may be infinite). 

In the finite dimensional quantum setting, with the phase-space resolution con- 
strained by the value of the Planck constant, it is clear that after an initial stage of 
positive entropy production the growth of the entropy of the evolved partition comes 
to an end as further iterations do not introduce any new information. When this 
stage is achieved (i.e., the saturation takes place), the rescaling factor dominates the 
asymptotics and the infinite-time limit is zero. 

The above analysis suggests the following solution. In order to get nontrivial 
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entropy results in the finite dimensional setting, instead of taking an infinite-time 
limit with a fixed Planck constant, one should consider simultaneous semi classical 
and long time limits computing the entropy production for each value of positive 
but eventually vanishing Planck constant only up to certain 'saturation' time. If 
the limits are taken in this particular way and order, the infinite time limit should 
be nonzero and should provide truly valuable information about chaotic behavior 
of the underlying system. In particular, the longest time scale on which the above 
procedure recovers classical entropy can be thought of as a kind of "breaking" time, 
after which classical and quantum evolutions no longer agree. Some recent preprints 
(THinillTj confirm the usefulness of the idea. In particular in jTHI the authors compute 
semiclassical limit of both CNT and ALF entropy productions (on logarithmic in h 
time scales) for cat maps and prove that both notions lead to the same classical result. 
In particular the results remain in full agreement with our results (see Theorem 15. 17j) 
providing the base for the interpretation of the quantum dissipation rate constant as 
quantum dynamical entropy of the system (at least in 2-dimensional setting). It is 
also worthwhile to note at this point that the agreement holds exactly in the situation 
when the dissipation time does not exceeds the "breaking" time of the system. 

The problem of "breaking" time is in some sense universal and emerges in all 
aspects of semiclassical analysis of quantum chaos. In order to understand it better we 
need to come back for a while to the description of the finite dimensional quantization 
procedure. 

After initial works and promising results by Hannay, Berry, Tabor and Keat- 
ing a need arose to clarify technical issues of finite dimensional quantization and in 
particular to generalize the scheme to allow quantization of at least all toral auto- 
morphisms (not only cat maps satisfying the 'checkerboard' condition) but eventually 
also Anosov maps and in general all canonical maps (symplectic diffeomorphisms) on 
the torus. 
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The first major breakthrough in this direction came in 1993 with a paper by 
Mirko Degli Esposti j331 in which the author suggested a representation-theoretical 
approach to the finite dimensional quantization (for more modern treatment see 
IS1E3)- Just hke in the infinite dimensional case the key ingredient of the quan- 
tization was the Weyl representation of the discrete Heisenberg group, with the dis- 
tinction that the representation was considered here over finite dimensional Hilbert 
space L^(T, /i), with purely atomic measure fi. The resulting *-algebras are indexed 
by rational values of h. Unlike the standard Weyl representation over the usual in- 
finite dimensional quantum Hilbert space L^(M), where by the Stone- von Neumann 
theorem there is essentially unique irreducible infinite dimensional representation, for 
any fixed rational h there are infinitely many irreducible Weyl representations over 
finite dimensional space (in particular L^(T, /i)). The representations are indexed by 
a parameter G T^ and all share the same dimension N {ii h = p/q then A^ is the 
smallest integer such that Nh is an integer jHI])- Using this "extra space" provided 
by the non-uniqueness of the infinite dimensional representations the author was able 
to remove the old checkerboard quantization condition and hence proved that for 
every cat map (i.e., 2-dim hyperbolic automorphism of the torus) there exists a fi- 
nite dimensional Weyl representation on which the corresponding quantum Koopman 
operator acts as an (inner) *-automorphism (see Section 15.2.11) . It was soon realized 
that the abstract Hilbert space L^(T, /i) can be considered as a generalization of the 
original Hannay and Berry space in a sense that strict periodicity of a wave function 
is replaced by a weaker condition of 6'-quasiperiodicity defined as follows 

^(g + m) = e2"'^^'"V(g), 

where m is an integer and 6p G T. An analogous condition is required for the momen- 
tum representation (with a corresponding Bloch angle ^g G T). For a given cat map 
F = [^ ^] the appropriate quantization condition on the joint Bloch angle 6 = {6q, Op) 
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has been established (see e.g. 1^31 170]) 

N I a ■ b 



2 ^c-d, 



Fe = mod 1, 



and after the work of Keating, Mezzadri and Robbins [70] it is usually referred to as 
'quantum boundary condition'. In the same paper the authors generalized the quan- 
tization scheme to arbitrary canonical maps on the torus (some additional technical 
issues of the multiplicativity of the unitary propagators implementing the dynamics 
had been addressed later by Mezzadri, Kurlberg and Rudnick jUOl [77] , but since they 
do not pertain directly to our study, the results will not be described here). 

In the meantime and afterward, different authors introduced finite dimensional 
quantization procedures for almost all other (e.g., discontinuous) well-known classical 
systems on the torus (see, e.g., [T31 11141 I115J for quantization of the Baker map and 
IH] for quantization of the Sawtooth map). 

Thus, after almost twenty years of development, the finite dimensional quanti- 
zation on the torus had finally been established as a general and rigorous procedure 
and successfully applied to a large variety of toral maps. 

In 1996 A. Bouzouina and S. De Bievre ^24. introduced an elegant and highly 
efficient mathematical formalism and notational setting in which the procedure could 
be described in a natural and transparent way. The authors started with a standard 
Weyl quantization on the usual Hilbert space L'^{W), extended unitary Weyl transla- 
tion operators (generators of the algebra) to the space of all tempered distributions 
S'{R) and then restricted their action to the space of quasiperiodic wave functions. 
The construction seems more natural then previous ones since instead of starting from 
an abstract representation theory, the appropriate algebra of quantum observables is 
derived directly from the traditional one. Thanks to this, any quantum observable 
corresponds in an natural way to its classical counterpart through its decomposition 
in the quantum Fourier basis formed by Weyl translation operators. 
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The method also provided a transparent hnk between the finite and the infinite 
dimensional settings (although the authors did not address this issue in their work 
directly). Indeed, if we denote by A the infinite dimensional algebra of the model 
introduced by Benatti, Narnhofer and Sewell JHl and restrict the values of the Planck 
constant to rationals then in the formalism introduced by Bouzouina and De Bievre 
we simply have (cf. Section I5.1.5|) 

A= / An{0), 

where A^id) denotes finite dimensional algebras of observables over the spaces of 
0-quasiperiodic wave functions (and in Degli Esposti's formalism AN{d) run through 
all the irreducible representations of the discrete Heisenberg group as varies in T^). 
Moreover as will be shown in Sections 15 . II and l5 .21 both finite and infinite dimensional 
quantizations can be constructed within the framework similar to the one introduced 
by Bouzouina and De Bievre, provided one defines in an appropriate way the Hilbert 
spaces of pure states to which the action of the Weyl operators is to be restricted. 

After clarifying quantization issues we can now return to the discussion of the 
breaking time. We did not define this notion precisely yet and in fact there is still 
no agreement in the literature on how it should be defined and even named. In an 
attempt to formalize this notion let us recall that according to Bohr's correspondence 
principle, the quantum evolution of any observable should approximate its classical 
evolution better and better as the Planck constant tends to zero. In mathemati- 
cal literature the principle is expressed in terms of the so-called Egorov property 
f |42[l64[llU9j ). which states that for any smooth classical observable / and the clas- 
sical evolution operator U, and for their corresponding quantizations Op{f), U the 
following estimate holds 

\\U'OpU)-Ov{U\m\<Ctjh. (4.8) 

In the above estimate the constant in general depends in an essential way on both 
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t and /. The key point is that it is ^-independent and hence in the classical limit 
the quantum and the classical evolutions "commute" . We want to stress that even in 
case of very reasonable quantization procedures the Egorov property cannot be taken 
for granted and has to be formally proved (see remarks after Theorem 4.2 in j21] for 
an appropriate counterexample, cf. also jHHl; for the most up to date results on the 
dependence of Ctj on t and /, see P5|). 

Using Egorov property we can introduce one possible definition of the breaking 
time by asking for the largest t = t{h) for which RHS of ()4.8p remains bounded. The 
corresponding time scale te is usually referred to as the Ehrenfest time which had 
been introduced many years ago pT|ll27j when physicists conjectured that for regular 
systems te should diverge as a function of /i at a power-law rate, while for chaotic 
ones, at a logarithmic rate. It turns out that the conjecture, although intuitively 
appealing, is extremely hard to verify on rigorous grounds, even for one of the simplest 
possible models of quantized fully chaotic systems; i.e., toral cat maps. Most of the 
results confirm the agreement between the quantum and the classical evolutions up 
to the logarithmic times. There is however still no general result stating that after 
the logarithmic time is passed, the quantum system diverges in its behavior from its 
classical counterpart. 

Below we briefly review the results obtained in this matter up to the present 
time. The first interesting observation is that the problem of estimating the Ehrenfest 
time is closely related to the problem of determining the semiclassical behavior of 
the eigenstates of the quantum system, on the one hand, and with the problem of 
finding the quantum period function on the other. The first problem (semiclassical 
behavior of eigenstates) is usually approached from the following two, slightly different 
perspectives. In the first approach one hopes to prove so called Quantum Unique 
Ergodicity (QUE), which states that a// eigenstates of a chaotic quantum propagator 
equidistribute (i.e., converge weakly to the Lebesgue measure) in the classical limit. 
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In the second approach one relaxes a httle bit the above requirement and aims at 
proving the so-called Schnirelman property f |118| IH3]). which states that if a quantum 
system has ergodic classical limit, then most of its eigenstates (with given energy 
level) equidistribute in the classical phase-space when h —>■ 0. QUE results are very 
difficult to prove. The first example of QUE for a quantum dynamical system with 
ergodic classical limit has been presented by Marklof and Rudnick in their paper 
|88j . published in 2000. The example is however very specific. The authors consider 
irrational skew translations of the two-torus. The key point which simplifies the 
matter considerably here is that the classical limit is not only ergodic but also uniquely 
ergodic, meaning that there is no other than Lebesgue measure on which quantum 
eigenstates can possibly concentrate on in the classical limit. Thus the problem 
of proving QUE in this particular case was essentially reduced to the problem of 
proving the Egorov property. Moreover irrational translations, although ergodic, are 
not weakly mixing and hence do not represent a typical example of a chaotic system. 
In fact QUE has not been proved yet for any quantized fully chaotic system !2B1 p. 4], 
and to the contrary in many cases, including as we will see below all cat maps, it 
simply does not hold. 

As far as cat maps are concerned the "closest" results to QUE had been proved 
either under the very strong additional restriction that the classical limit {h = 1/N -^ 
0) is taken only over special sequences of reciprocals of primes A^ (and under the 
assumption that generalized Riemann Hypothesis holds) |13], or under weaker condi- 
tions on h (density one sequence in N) but with quite restrictive number-theoretical 
conditions on admissible cat maps [TH], cf. also [77j. In both cases restrictions aim 
at exactly the same point - to ensure the absence or sufficiently slow growth of the 
degeneracies of the corresponding eigenvalues, which corresponds to long (of order A^) 
periods of the quantum propagator. Rapidly growing degeneracies of the eigenspaces 
of the propagator caused by short (e.g., of order InA^) quantum periods existing 
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for some particular values of the Planck constants allowed for construction of the se- 
quences of eigenstates which concentrate in the semiclassical limit on measures with a 
nontrivial pure-point component. The existence of such sequences - the phenomenon 
referred to as strong scarring - has recently been proved by Faure, Nonnenmacher 
and De Bievre in j33] for an arbitrary hyperbolic toral automorphism in 2-dim (for 
some results in higher dimension see |2S])- Thus QUE does not hold in these cases. 
The scarred eigenstates cannot concentrate however totally on pure-point measures. 
In jHSI Faure and Nonnenmacher show that the Lebesgue measure component of the 
support of any such state must account for at least 1/2 of the total weight (and the 
bound is sharp). Moreover, the sequences of scarred eigenstates are very exceptional 
(in the sense that to construct them one needs to choose these exceptional values of 
h for which the quantum propagator has minimal, i.e., logarithmic in h period). The 
conclusion is that although QUE in general fails, some version of the Schnirelman 
property usually holds for general chaotic quantum systems. For cat maps the prop- 
erty had been proved already in [23 • It has been also confirmed in a wide variety 
of other contexts (e.g. for ergodic geodesic flows on compact Riemannian manifolds 
|1181 11281 I3B] . for ergodic billiards |SHl I13UJ . and recently for maps with mixed dy- 
namics [H7j). 

As far as the Ehrenfest time is concerned the discovery of the scarred states 
turned out to be directly connected with the construction of an example in which the 
breakdown of the classical-quantum correspondence happens sharply on a logarithmic 
scale. The appropriate example was constructed within the framework of cat maps 
by Bonechi and De Bievre in |^|. The main idea was to consider the action of the 
quantum propagator on coherent states supported on the region in the phase space 
of the diameter not exceeding \/h. For times tending to infinity but no faster than 
te = (27)"^ ln(^~^) (7 is the Lyapunov exponent of the cat map) the whole support 
of the state must shrink and in the classical limit the state's evolution concentrates 
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on the classical orbit of its center - the behavior shared by both classical and quan- 
tum evolutions (the Wigner function of the coherent state converges weakly to the 
delta function). On time scales slightly longer than te strong mixing properties of 
the classical map prevail and cause the support of the coherent state to stretch and 
the resulting semiclassical limit corresponds to a constant function over the whole 
phase-space (the Wigner function of the coherent state converges to 1). For classical 
dynamics this scenario will continue regardless the length of the time scale. To the 
contrary, the quantum system is periodic and the initial concentration phenomena 
must repeat itself once the period of the map is completed. At that point quantum 
and classical evolutions depart from each other. It is then enough to find a sequence 
of Planck constants for which the quantum period function is of order In(^). This 
is exactly the minimal possible period and as was mentioned above the appropri- 
ate sequence can be constructed for any cat map yielding the desired logarithmic 
asymptotics of the breaking time. 

A few words of caution are necessary here. First of all, the above described phe- 
nomena need not reflect general properties of quantum systems including cat maps. 
Indeed, according to earlier described results by Keating, Kurlberg and Rudnick (see 
14.21 and 14. 3|) for a great majority of the values of h, quantum period is much longer 
than In(^) and hence for 'most' of classical limits the phenomenon may not be visible 
on this time scale at all. As the time progresses more and more trajectories will start 
to complete their periods and the breaking of the correspondence, understood as a 
statistical phenomenon, may happen on longer time scales. 

It is also important to distinguish here between two different semiclassical ap- 
proaches to the question of determining when the breakdown between classical and 
quantum mechanics occurs. The above- described results belong to the category of 
direct 'quantum-classical' limits. 

Another approach is to study the validity of semiclassical approximations to 
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quantum propagators. Here the situation is different and the correspondence on 
times scales much longer than logarithmic have been observed in several cases. For 
numerical results in these directions see IV21\ \V2']l llUlj . These results cannot shed 
any light on the question of the breaking time for cat maps, however, since in this 
particular case the semiclassical approximation to the quantum propagator is exact 
and hence no breaking time can be observed through this analysis. 

A general conclusion is that depending on the observed property and the level of 
the intensity of the phenomenon on which the correspondence is tested one can expect 
that a whole spectrum of different breaking times may exist ranging from logarithmic 
to power-law scales in h. In Section 1^^ Corollarv 15 . 1 81 f point IV) we give an example 
of the situation when the quantum and classical noisy toral automorphisms behave 
in the same way from their dissipative properties point of view, while the evolutions 
have already exceeded the Ehrenfest time. It would be of particular interest now 
to see whether 'CNT-entropic' breaking time coincides with 'ALF-entropic' breaking 
time and how these two are related to Ehrenfest and dissipation time scales. The 
results which we are going to derive in the next chapter are closely related to this 
question. The complete solution is however still elusive. 
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Chapter 5 



Quantum dissipation time 



The main goal of this chapter is to introduce the notion of quantum dissipation time 
and to study its semiclassical and noisy asymptotics for quantized toral maps. The 
first step is to construct appropriate quantum model of the phase space fSection lS.lj) 
and of the dynamics (Section 15. 2|) so that it can be considered as a quantization of 
the classical systems studied in Part I. In the second step (Section 15. 4j) we will then 
focus on definitions and semiclassical analysis of the corresponding quantum notions. 

5.1 Quantization on the torus. The kinematics 

In this section we present in a systematic and rigorous way kinematic step of the 
quantization of classical systems on the torus. That is, we construct appropriate 
Hilbert spaces of pure states and algebras of observables. In the following section we 
will concentrate on the dynamical step i.e. on quantization of canonical toral maps. 
The presentation is based on the approach by Bouzouina and De Bievre |24j and 
generalizes it in two directions: 

1. We quantize the systems with arbitrary phase-space dimension (for equivalent 
approaches to quantization in multidimensional setting see e.g. |lU8t l61p. 
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2. The quantization scheme is extended in such a way that both finite and infinite 
dimensional settings are included as particular cases. 

Throughout this chapter we will use the terms Model I and Model II to refer respec- 
tively to finite and infinite dimensional quantization schemes. 

5.1.1 Weyl quantization on W^ 

We briefiy recall the standard Weyl quantization of systems with d degrees of freedom 
and with M^'^ as a classical phase space (for systematic presentation see 0). 
As before, h will denote the Planck constant. Whenever convenient we will also use 
the notation h = ^. 

ZTV 

Consider usual quantized position and momentum operators 

Q = iQi,...,Q,), P = (Pi,...,P,) 

on the Hilbert space L^(M'^) of square integrable wave functions 

dip 
Qjip{x) = Xjip{x), Pjip{x) = -i^-^ — (x). 

OXj 

The domains of these operators contain the Schwartz space iS(M'^) and hence are 
dense in L^(R'^). Although all are essentially selfadjoint, they do not possess a com- 
mon domain of selfadjointness. In Weyl quantization one considers exponentiated 
versions of the above operators, which helps to avoid the domain problems and, more 
importantly, introduces in a natural way the notion of quantum translation operators 

The explicit action of Uq and Vp on L'^{M.'^) is given by 

{Uqij){x) = ^{x - q), {Vp'il:){x) = e^P-^VW- 

The commutation relations for Uq and Vp follow easily 

UqVpiPix) = e^P-^'-^Vl^ - q), VpUqiPix) = e^P'Xx - q). 
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and hence 

UgVp = e-i'^WpUq. (5.1) 

Using BCH formula and the fact that [p ■ Q^q ■ P] = ihq ■ pi we also get 

UgVp = e-49Pet{pQ-9P) 

Let V = {q,p) G M^*^. The Weyl translation operators T^ are defined as the following 
symmetrized products of Uq and Vp 

Using the notation convention X = {Q, P) and v A X = p ■ Q — q ■ P the Weyl 
operators can be written in a compact form 

T, = e^"^^. (5.2) 

T^ can naturally be extended to the space of tempered distributions 5'(M'^). 
The explicit action of T^ on a (generalized) wave function ip G iS'(]R'^) is given by 

{T^ij){x) = e^''-PesP-(^-9)^(a; - q) = esP-(^-'^/2)^(a; - q) (5.3) 

and yields the following property 

T^T^, = e^n-^^'T^+^,, (5.4) 

which in turn implies Weyl-type Canonical Commutation Relations 

T^T^, = e>^"'T^,T^. (5.5) 

5.1.2 The space of pure states on T'^. 

In this section we determine a quantum analog of the notion of periodic phase-space 
of a classical system by introducing the space of generalized quasi-periodic wave 
functions (distributions). 
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The idea of quantization on the torus will be reflected here in the requirement that 
the wave function be quasiperiodic in its position or momentum representation. 
A distribution ip G iS'(M"') is called quasiperiodic in position representation if there 
exists a constant 6p G T'^ such that for any mi G Z'^ 

The set of all such distributions will be denoted by Sj^ {Op). 

Similarly a distribution ip G iS'(]R^) is called quasiperiodic in momentum representa- 
tion if there exists a constant Oq G T'^ such that for any 7712 G 7/ 

where J^h denotes the quantum Fourier transform 

(•^^^)(P) = -f^j^me-'^'"^dq. (5.6) 

The corresponding set will be denoted by Sf^ (Og). 

In our approach the space of all admissible pure states denoted by HhiO) of the 

quantum system on the torus will always be understood as a linear space generated 

by (not necessary all) elements oi Sf^ {dp)USf^ (Og) and equipped, according to usual 

requirements of quantum mechanics, with some Hilbert structure. The choice of the 

space will depend on the model one intends to work with. 

We note that if ■?/' G S)^ {Op) then its momentum representation J^hi^ has a discrete 

uniformly /i-spaced and ^p-shifted support and hence is represented by a Dirac delta 

comb (or more generally brush) of the form 

where c = {cs}sei'^ denotes a sequence of complex numbers. Thus any distribution 
■4) G Sy^ {Op) is uniquely determined by a triple {c,h,Op). We will sometimes write 
ip = (c, /i. Op). Similar remark applies to any element of Sj^ {Oq). 
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One easily notices that for any tp G S^ [Op] and any mi E U^ 

Tim.M = e-2-^--"^^^ (5.7) 

Similarly for any ^ G Sj^ (Og) and any m2 G Z"* 

T(0,rr^,)V' = C^-^^'^^^V'- (5-8) 

That is, the spaces Sfl' (dp) and Sf^ {dq) consists of eigenstates of integral translations. 

5.1.3 The algebra of observables on T"^. 

In this and subsequent sections we will frequently use the following terminology. 

Definition 5.1 An algebra A is called 

i) * -algebra if it is equipped with an involution * : ^ i-^ ^, A** = A, A E A. 

a) B* -algebra if it has a Banach space structure and \\AB\\ < ||A||||i?||. 

Hi) C* -algebra if it is a B*-algebra that satisfies \\A*A\\ = ||v4p, for all A E A. 

iv) H*-algebra if it is a B* -algebra with an inner product satisfying 

{A,BC*) = {AC,B) = {C,A*B) 

(e.g. {A,B) = t{A*B) if t is a faithful tracial state). 

v) W* -algebra (von Neumann algebra) if A!' = A, where A" is a bicommutant. 

The Weyl quantization on the torus will consist in the restriction of the action of 
Weyl translations from the whole 5'(M°') to 7ih{0). 

This restriction is well defined only if translations Ty preserve S^ {Op) and S^ {6q) 
spaces. In view of (15. 7j) and (j5.8|) the condition is equivalent to the commutativity of 
T^ with integer translations. The latter condition can be stated as follows 

Proposition 5.2 

T^ : S''^\ep) ^ Sl^\0p) ifJvER'^x hZ'^. 
T, : SJ:['\eq) ^ SJ^\eq) iffvE hZ'^ X M^. 
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Remark 5.3 If 1-Lh{0) has nontrivial intersection with Sf^ [Op] and S^ [Oq] then 

T^:nh{9)^nH{e) zjfvehz^''. 

Motivated by the above considerations we define microscopic quantum phase-space 
translations acting on the space Hh{0) 

Wk := T,fc = 62-*=^^. (5.9) 

The operators Wk are indexed by the elements of the integral lattice Z^*^ and can be 
thought of as quantum counterparts of classical Fourier modes. The corresponding 
commutation relations are now given by 

WkWrr, = e-^^'»^^l^fc+^, (5.10) 

WkWrr^ = e^^''""'"'W^Wk. (5.11) 

Note that the placement of the Planck constant has changed after this rescaling (cf. 

dEl and (IH3|) V 

The formal algebra of observables of our quantum system is generated by the set of 

operators {Wk}kGZ2'' acting on TihiO) and will be denoted by Ah{0). We note that 

Ah{0) is a *-algebra with the involution defined by W^^ = W-^. Equipped with the 

standard operator norm Ah{0) would become a C*-algebra and its weak closure - a 

iy*-algebra. 

The reference state defined on the generators of Ah{0) by 

T{Wk) = 4,0 

corresponds to the standard Lebesgue measure in the classical system and can be 
uniquely extended to a tracial state on the whole algebra Ah{0)- 
In our approach, instead of studying Ah{0) as a C* or iy*-algebras, we prefer to take 
advantage of its natural if *-algebraic structure, associated with the state r. If r is 
faithful the inner product on Ah{0) is introduced as follows 

{A,B) = t{A*B). 
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In particular 



(W^fc, Wrr.) = 4,r 



and hence automatically {Wk}f,^^2d becomes an orthonormal basis for Ah{0)- 
This completes the discussion of the kinematic step of our general quantization 
scheme. Below we derive as particular cases two models described in Chapter 0] 
We call them here Model I and Model II. 



5.1.4 Model I 

In this section we derive the original Hannay-Berry model including all its later de- 
velopments and generalizations described in Chapter EJ To this end we set 

Thus a distribution ip G iS'(]R'') belongs to TiKO) iff for any m eT, 



2d 



{J^hip){p + rn2) 



2TTi0p-m,i 



^(g), 



-2TTi0q-m/2 



(^h^)(p) 



The parameters Op, 6q are sometimes called Floquet or Bloch angels. It is well known 
(see e.g. pij) that the set of all distributions satisfying the above mutual quasiperi- 
odicity conditions contains nonzero elements iff 






(5.12) 



where N G Z+. In literature ()5.12|) is sometimes referred to as the Bohr-Sommerfield 
condition. Thus whenever we discuss the Model I we assume that ()5.12|) holds. 
We pause for a while to introduce some notation. Let 



J_ : = 



-/ 
/ 



J^:-- 



I 

1 
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where / denotes the identity matrix on M'^. 
For any pair of vectors v,v' E M^^ we set 

V Av' := vJ^v', V y v' := vJ+v'. 

Of course i? A t;' is a standard symplectic product on M^'^. 

Using the above notation quasiperiodicity conditions can be naturally encoded in 

terms of the action of translation operators 

Proposition 5.4 Let ip G S'{R'^), then ip G ni{e) iff for all m G Z^'^, 

Trr^V' = e2"<^"'^"'+"'^^)^. (5.13) 

Proof. Obviously if ()5.13|) holds then ip is 0-quasiperiodic (consider m = (mi,0) 
and m = (0, 1712)). On the other hand ii ip E Ti-l^O) then 

_ 27r^(Af/4)T7^VT7^ -27^^0p■T7^l 27rJ0,-TTi2„;. _ 2ni{N/4,m\/m+mA0) 1 h 

The general form of an element of 'Hj^{6) can easily be determined. One finds (for 
detailed derivation see Section lA.21 of Appendix IX|) that ip G Ti-l^O) is necessary a 
quasiperiodic Dirac delta comb (brush) of the form 

^(9) = ]^ Yl (^^^^+e,/N{q), (5.14) 

where c^ is a quasiperiodic sequence of arbitrary numbers supported on Z'^/N lattice 
and satisfying Cs+n = e^'"''^^"^Cs. Thus, although S^^'{Op) and S^^'{Op) are not Hilbert 
spaces, Ti.l{0) can naturally be identified with the Hilbert space C^ , by introducing 
the following L^-norm 

2 •- jyd Z^ l^^l • 
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We note that, the crucial difference between the present model and the one considered 
in the next section lies in the fact that full quasiperiodicity of the state space Ti.{{d) 
implies also quasiperiodicity of quantum Fourier modes. Indeed we have 

Proposition 5.5 For any m G Z^'^, 

W^fc+iVrr. = e^^C^'^^'^^M^fc, (5.15) 

where 

1 A^ 
a{k, m, 0) = -k Am -\ m \/ m + m A 6. 

Proof. Using Proposition 15.41 we have for any ip G 7Y((0), 

_ ^TTikAm^2ni{N/4mVm+mAe)Tk/M ^ _ g27ria(fc,m,0)|y / g 

The algebra of observables of Model I will be denoted by Aj^f{6). Due to the quasiperi- 
odicity of the set of its generators Aj^f{6) is finite dimensional and as a linear space 
can be identified with C{Hi{e)) = M^^y^N^i = C^'" . 

We note that Ajy{6) is still a *-algebra with the involution defined by W^ = W^^. 
The fact that operation * is consistent with the quasiperiodic structure of the set of 
generators follows from the property «(— fc, — m, 6) = —a{k, m, 6) mod 1. 
The tracial state on A]^{6) can be defined now explicitly 

And as above {Wk, k G Z^} (Z^ := Z^"^ mod A^) becomes an orthonormal, but this 
time finite basis for an if*-structure on Aj^{6). 

The H*-noTm. on A^j^{6) will be denoted by || • \\hs- One needs to keep in mind 
that II ■ Whs does not coincide with the standard operator norm, hence Alf{0) is not 
(considered here as) a C*-algebra. 



Section 5.1. Quantization on the torus. The kinematics 137 

Now we choose the fundamental domain of periodicity Z^ for our quantum Fourier 
lattice. The choice centered around the origin seems to be the most natural one (cf. 
|lUUj ). That is, we assume that if fc = (fci, ..., k2d) & Z^ then for every j G {1, ..., 2d} 



kje 



{-N/2 + 1, ..., A^/2}, for N even 

{-{N - l)/2 + 1, ..., (A^ - l)/2}, for A^ odd. 



The set {Wk, k G Z^} forms an orthonormal basis for An{0)- 

To any classical observable / G C°°(T^'*), or more generally / G L^(T^'^), with 
Sfc 1/(^)1 < °° there corresponds an element of A]^{6) i.e. its Weyl quantization, 
denoted by OpN{f), and defined in terms of its Fourier expansion 



OpNif) 






fceZ2d k&jf \TngZ2'i 



The map Op^ is not invertible. One can define however the isometry W^ : Aj^f{6) \-^ 
L^(T^'^) which associates with each observable A G Aj^{0) its polynomial Weyl symbol 



w^{A) = y: 



akWk, 



feezed 



where Wk '■= e^'^*'"^^ denote classical Fourier modes and a^ = {Wk,A). The operators 
OpN and W^ are inverse of each other when the domain of Op^ and the codomain 
of W^ are restricted to In = W^o Opn{L^{T^'^)). 

5.1.5 Model II 

In their original paper [1^ the authors did not specify any particular Hilbert space 
of pure states for their model and instead performed the quantization starting on the 
algebraic level. 

In order to emphasize a close link between the two models within the framework 
considered here, we first construct the physical Hilbert space for this model and then 
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show that the natural restriction of Weyl quantization on 5'(]R'^) to this space yields 
the algebra of observables considered in jTHj . 
First we define the following two Hilbert spaces 

nt\e,) := {{cKO,) G s^\e,) ■. c e /^(z^)}, 
H^:\e,) := {(c, h, e,) G st\e,) ■. c e I'iz')}. 

The values of the parameters Op and Og do not play any significant role in this model 
and hence we will be working with the following spaces 

As as we proved in the previous section, if (c, /i, 0) G 5^'' (0) fl Sj^ {Qi) then either 
c = or c ^ l'^{l/)}, since any such c must be periodic. Thus Ti}^ fl nf = {0}. 
We then define 

In opposition to 7i(, the space Tij^ is infinite dimensional. 

Now similarly as in Model I the quantization consists in the restriction of the standard 

Weyl quantization on iS'(M°') to the space 7i^^. The construction of the space 7i\^ 

insure that the set of all admissible quantum translations is discrete and coincides 

with the family of operators Wkki^T?^ considered in the previous section (see Remark 

Q. 

As it was mentioned above (cf. ()5.10|) ) Wk satisfy the relations 

The Planck constant h plays the role of deformation parameter Q considered in [T^ 
(unrelated, of course, to our &). Putting h := 26 one recovers exactly the relations 
assumed in |I2j. The algebra of observables of model II Aj,^ can now be defined as 
a *-algebra generated by elements Wk- Taking the weak closure (the bicommutant) 
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in BiTil^) yields a iy*-algebra isomorphic to the algebra of observables introduced in 

m- 

For all h > 0, A\^ is always infinite dimensional and hyperfinite i.e. is generated by 
an ascending sequence of finite dimensional algebras. However the internal structure 
of A\^ depends in an essential way on h is rational or not. Before we discuss the 
classification of A]^ we need to introduce some definitions jTTj . 

Definition 5.6 

i) A W*-algebra is called a factor if its center is trivial (the algebra is irreducible). 

a) A factor is of type I if it is isomorphic to B{7i) for some Hilbert space Ti. 

Hi) If dim{7i) = N then the factor is called to be of type In- 

iv) Two projections Pi, P2 in a W* -algebra are said to be equivalent if there exists an 

element A in the algebra such that Pi = A* A and P2 = AA* . 

v) Projection is said to be finite if it is not equivalent to any proper subprojection of 

itself. 

vi) A factor is called type IIi if the identity operator 1 is finite. 

If h is rational then A^' factorizes into its nontrivial commutative center and a type 
/jyd factor (matrix algebra) 

Ai' = L\T^'')®Mn,^n, 

where A^ is the smallest integer such that hN G Z_|_ (see j73] and [HI] for transparent 
proofs). Moreover it can be shown (cf. [23) that in this case The algebras of both 
models are related by the following decomposition formula 



JT2 



If h is irrational Aj^^ is a factor of type III (cf. 

In both cases Ah admits a faithful tracial state, which in case of rational h factorizes 

into standard Lebesgue measure on T^*^ and the normalized trace on Ai^d-^Nd. 
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Moreover it can be shown that 

We end this section with a few remarks regarding the choice of the Hilbert space for 

the Model II. 

Remark 5.7 The choice of the physical space for Model II is not unique. One can 
consider e.g. the space Hj/ := Tijf ® 7i)f and the algebra generated by elements 

Remark 5.8 In contrast to the case of Model I, for all irrational h, the algebra of 
Model II is of type III is never isomorphic to any full (i.e. type I) algebra BiTih) 
regardless the choice of the Hilbert space Tih. 

5.2 Quantization on the torus. The dynamics 

In this section we perform the dynamical step in the process of quantization on the 

torus. We consider arbitrary canonical toral maps i.e. area and orientation preserving 

homeomorphisms on T^'^. 

Let $ denote such a map. As was already mentioned in Section 13.1.11 $ can be 

decomposed into the product of three maps $ = F o t^ o $i, where F G SL{2d,Z) 

is a symplectomorphism, t^ denotes a classical translation by vector v and $i(a;) = 

X +p{x), where p is an arbitrary zero-mean, continuous and periodic function. 

We assume that $i represents a time-1 flow map associated with a periodic Hamil- 

tonian. In 2-dimensional case this assumption is equivalent (cf. (TO]) to the above 

requirement that p be of zero mean (w.r.t. the Lebesgue measure on T^'^). This 

can always be achieved by replacing P with p' = p — (p) and adjusting accordingly 

translational component of $. 

To quantize $ one first quantizes F, t^ and $1 separately. The quantization of $ is 

then defined as a composition of corresponding quantum *-automorphisms 
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The quantization procedure will be prescribed in such a way that the correspondence 
principle will hold. In the case of Model I this will be expressed in terms of so called 
Egorov property, which states that for every / G C°°{T^'^) there exists C > such 
that 

\\U{0pM))-0pN{Uf)\\<^, (5.16) 

where Uf = / o $ is a classical Koopman operator of $. 

5.2.1 Quantization of symplectomorphisms 

Here we describe the quantization of toral symplectomorphism i.e. symplectic au- 
tomorphism on T^"'. In classical setting the action of a toral automorphism on the 
algebra of classical observables was defined by means of the Koopman operator given 
by Upf = f o F, where / was an arbitrary observable, usually assumed to be an 
element of L°°(T^°') or C°°(T^^). In our approach however it was more convenient 
to consider the action of Up on a slightly bigger space L'^{T'^'^). The natural Hilbert 
structure of the later allowed us to determine the dynamics by specifying it on the 
basis of classical Fourier modes Wk{x) = e^'^**'^^ 

where F'^ denotes the transposed map. Setting F' = J_F'^J_ we get 

The most natural way to define the quantum counterpart of this dynamics is to 
consider the formal superoperator version Up of the classical Koopman operator 

KpWk = Wf^u- (5.17) 

Such dynamics will be well defined i.e. will define a *-automorphism of Ah{0) only 
if the action of Up is consistent with its algebraic structure (Weyl commutation rela- 
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tions). We note that 

^,iTihkAm^—iTihF'kAF'mTjr ti/ 

— e e yVp'kyvF'm 

Thus in order for Uf to be a *-autoniorphisni the following condition has to be satisfied 

2 
F'kAF'm = kAm mod-. (5.18) 

ft 

Since the condition (j5.18|) has to be valid for arbitrary small h we have 

Proposition 5.9 If a map F G SL±{2d,Z) is quantizable then it is symplectic. 

Symplecticity is then necessary for quantization regardless the model i.e. the space 

7ih{0) one choses to work with. We note that for symplectic maps F' = F~^. 

Depending on the choice of Tih^O) the condition may be also sufficient. Indeed, 

this is exactly the case in Model II. In view of the lack of additional conditions on 

the generators of the algebra A\^ there are no quantization restrictions other than 

symplecticity of the map. 

In some cases, however, the structure of Ti.h{6) impose additional relations on the 

generators of the algebra Ah{0) and then additional conditions are needed. This is 

the case in Model I, where the assumption of qasiperiodicity in both position and 

momentum resulted in qasiperiodicity of algebra. 

Thus in this case we have to ensure the compatibility of the action of Uf with the 

quasiperiodic structure of A^^{0). To this end we note that for a given quantum 

Fourier mode Wk+Nm, one can compute the value of UFWk+Nm in the following two, 

in general different, ways: 

- on one hand using linearity of F' we have 

UfWu+Ntti = WF'k+NF'm = & "" ' ^' ^F'k 
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- on the other hand, by hnearity of Up 

That is, for every A; and m in Z^^ we must have 

a{F'k,F'm,e) = a{k,m,e) modi. (5.19) 

The map F will be called quantizable in Model I if for every A^ G Z+ there exists 
6 G T^*^ (possibly depending on N) such that (j5.19|) holds for all k and m in Z^*^. 
Below we summarize the quantization condition in Model I for general toral auto- 
morphisms. 

Proposition 5.10 A toral automorphism represented hy F E SL±{2d,'Z) is quanti- 
zable iff it is symplectic. For any given h = N"^, the corresponding has to satisfy 
the following condition: 

N ( A-b\ 

— \ +Fe = e mod 1, (5.20) 

2 \C-DJ 

where A, B, C, D denote block-matrix elements of F, that is 

A B 
C D 

and A ■ B denotes a contraction of two matrices into a (column) vector, defined as 
follows 



{A-B\ = Y,Ai3B^r 



The existence of solutions of equations of the type ()5.20|) is easy to establish. We 
note that if N is even then one can simply choose = 0. The same solution can be 
chosen whenever the vector 

' A-B 
C -D 
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has all even components (and this condition reduces to Hannay and Berry's 'checker- 
board' condition stated in ^^ in d = 1 case). Otherwise one considers two cases. If 
F — / is invertible then for any k G Z^'^ 

•="-"- (ft::)" 

There are exactly | det(F — /)| distinct solutions. In particular the solution is unique 

if F — / G SL±{2d, Z) (in c? = 1 case det(F — /) = 2 — Tr F, hence uniqueness holds 

iff TrF = 1 or 3). If F — / is singular one can decompose the matrix F into an 

identity and nonsingular block and construct an appropriate 6 by applying the above 

considerations to each block separately. 

We want remark that in view of the defining condition (|5.17|) . the Egorov property 

()5.1fi|l is automatically satisfied (with no error term). 

We end this section with a few comments about the above quantization conditions. 

The conditions were imposed to ensure that Up is a *-automorphism of the algebra 

In Section 1331 we will prove the following simple 

Proposition 5.11 Any *- automorphism on finite matrix algebra is inner. 

Thus in view of this proposition our conditions ensure the existence of a physical 
quantum propagator Up implementing the dynamics on the underlying physical space 

UpA = ad{UF)A = WpAUp. 

It can be shown that tJp (which is only determined up to the phase factor) coincides 
with the propagator introduced by Hannay and Berry in their original paper [63]. For 
further details regarding this construction we refer the reader to [HUj . 
We also note that there exists a geometric interpretation of these quantization con- 
ditions. Indeed, consider the evolution of the Wigner transform of a given wave 
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function under the cat map dynamics. It is Well known that the Wigner function 
evolves according to a classical map (see Appendix ^ and forms a {2NY x [2NY 
periodic delta brush supported on the half-integer lattice (if N is length of the side 
of the fundamental domain of its periodicity). Symplecticity insures that the evolved 
delta brush represents once again a Wigner function. 

In case of odd A^ or when the wave function is 0-quasiperiodic {0 ^ 0) which means 
that the supporting lattice of Wigner function is shifted by on the coordinate 
plane, one wants to ensure that this supporting lattice remains on the same place 
throughout the evolution. The latter property is equivalent to condition ()5.20|) . and 
can be thought of as the conservation of the initial 'quantum boundary conditions' 
(see ^). 

5.2.2 Quantization of translations 

As explained in Section 15.1.11 a translation t^ is quantized on L^(M'^) through a 
Weyl operator T^. We have noticed that such quantum translations act inside the 
algebra A]^{0) only if i; G hZ'^'^. If this condition is not satisfied then quantization 
depends on the model. In case of Model I there are several possibilities to quantize 
such translations [21]. We will choose the prescription given in jSH]: we take the 
vector v^^^ E N~^Zi'^^ closest to v (in Eudlidean distance), which can be obtained 



/) 



by taking, for each j = l,...,2d, the component v, = ■'-^, where [x] denotes 
the closest integer to x. One then quantizes t^ on Hn{G) through the restriction of 
T^{N) on that space (this is the same operator as W[nv]{N,0)). The corresponding 
*-automorphism on AN{d) is given by 



Ut, = ad{T^ 



(N) , 



It was proved in jHH] that the Egorov property ()5.16p holds for this quantization. 
In case of Model II the situation is simpler. Even though T^ itself does not act as 
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an inner *-autoniorphisni on Aj^^ it introduces an external *-autoniorphisni on this 
algebra 

UuWk := T;WuT^ = e-'™^'=W^fc (5.21) 

and this automorphism can be taken as a quantization of t^- Note that definition 
given by ()5.21|) cannot be applied in Model I since the action of Ut^ is not consistent 
with quasiperiodic structure of ^^(^). 

5.2.3 Quantization of time-1 flow maps of periodic Hamiltonians 

We present here the quantization in case of Model I. Let $i denote the time-1 flow 
map associated with a periodic Hamiltonian H{z,t), meaning that $( : T^*^ —>■ T^'^ 
satisfies the Hamilton equations: 



dt 



V^H{^t{z),t), $0 = /. 



To quantize $i, one applies the Weyl quantization to the Hamiltonian H{t), obtaining 
a time- dependent Hermitian operator OpN^e{H{t)). From there, one constructs the 
time-1 quantum propagator on 1-Ln{6) associated with the Schrodinger equation for 
OpN,e{H{t)): 

UN,e{'^i) := Te-^'^^^^o OPiv,e(H(t))dt 

(T represents the time ordering). As above, the corresponding *-automorphism on 
A^j^{d) is defined by 

UiA = adiUN^^i)) A = U*^^e{^i)AUN,e{^i). 

For such a propagator, the Egorov property is holds, yet with a nonzero error (see 
|lU9j . Theorem IV. 10 for the case of a time-independent Hamiltonian on R^*^). 
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5.3 Quantum noise 

In this section we construct quantum noise operator for our system. We first recall 
some standard facts regarding quantum noise in general (cf. |5].|32j). 
The influence of a noise on a quantum system is described through the interaction 
between the system and its environment. If the system is in a state given by a density 
matrix p and the state of the environment is penv then the state of the open system 
(universe) is given by their tensor product p ® Penv As a principle the quantum 
evolution of the whole universe is assumed to be unitary. Thus there exists a unitary 
operator U such that the evolved jointed density is given by Up ® penvU*. The noisy 
quantum evolution F of the small system is recovered by tracing out the environment 

T{p) = Tr,^,{Up(S> PenvU*). (5.22) 

In general, F is not unitary. It is however always (I) trace preserving and (II) com- 
pletely positive. Complete positivity implies positivity which together with trace 
preserving property ensures that F maps densities into densities. 
Usually the environment is not specified explicitly and the noisy evolution of the 
system is described entirely in terms of some abstract operator F acting only on 
the small system. Although positivity and trace preserving property would suffice to 
insure that such operator preserves densities, it would not insure its representability in 
the form (j5.22|) . The sufficient condition for such representation to exists is complete 
positivity (which provides robustness of positivity of F w.r.t. tensor products). Any 
quantum noise operator is thus characterized by and should satisfy properties (I) and 
(II). 

By Kraus theorem condition (I) is equivalent to the fact that F admits the following 
operation-sum representation 

T{p) = J2GkpGl (5.23) 
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where Gk-, called operation elements, are arbitrary bounded operators. The above 
representation is not unique and in some cases it is useful to use continuous parameter 
(integral) representation 

r(p)= [dkGkpGl, (5.24) 

Jk 

The corresponding dynamics on the observables of the system is given by 

A^J2 GlAGk A^ f dkGlAGk. 
k -'^ 

Condition (II) is then equivalent to the requirement that the family Gk constitutes 

operational partition of unity 

Y,GlGk = Id. (5.25) 

k 

We now proceed with the construction of the quantum analog of the classical noise 
operator introduced in Section 12.1.21 We start with Model I and follow the quan- 
tization method used in |1U(JJ . For given noise generating density g we will use the 
following notation 

^^(*) •= ^d^ (~) ' ge{x):= ^ g,{x + n). 

Also the following notation regarding Fourier transform will be utilized 

lik) := / ~gix)e-'-''^"dx, 
It is easy to check that ge{k) = ge{k) = g{^k) for k G Z^'^ and hence 
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Let US note that the classical noise operator G^ on L^(T^'^) can be represented as 
follows 



GJ= I U'"){tJ){x)dv, 

JT2d 



where, as before, t^ denotes the Frobenius- Perron operator of the classical phase space 
translation. A natural way of quantizing the noise operator is to formally replace t^ 
with quantum phase space translations T^. However due to the discrete character of 
quantum translations one has to discretize the classical operator before quantization 






k&?^ 



and then quantize it by introducing the following superoperator 



rZ^£d j^j K.^£d j^j 



where, as always, ad{Wk)A = W^AWk for any A G Aj^{6). 

The role of the prefactor -| is to insure that Qe,N is trace preserving. One can easily 

check (see Appendix 15. 5|) that Z = N'^'^g^]\f{0) and that moreover 

Proposition 5.12 Q^^n is a completely positive trace preserving map and admits the 
following spectral representation on A^j^{0) 

g,,NA = J2 leNik^)akWk, (5.26) 

fcez2d 

where 

W0:=%^^, A := $^a,iy,. (5.27) 

Defining the subalgebra of observables orthogonal to identity 

A%{e) = {Ae A'^{e) : ao = 0} 
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and introducing the superoperator norm w.r.t. iJ*-norni on A%{0) 



?,,iv|| := sup \\Qe,NA\\HS 



one immediately gets 



||^,,Ar|| = max 7eAr(fc), 

which in particular means that the quantum noise operator acts as a strict contraction 

onA%{e). 

We end this section with a brief comment regarding the quantization of the noise 

operator in the case of Model II. In this case due to the lack of quasiperiodicity 

conditions and the uniform in v representation of the action of quantum translations 

on Al^ (see (|5.21|) ) the quantum noise operator acts in an isomorphic fashion to the 

classical one: i.e. we can introduce a continuous version of Kraus noise operator 



^eA 



g^{v)ad{T^)dv, 



which establishes uniform w.r.t. h G IR+ isometry between classical and quantum 
noise operators. 

5.4 Semiclassical analysis of quantum dissipation time 

For any canonical map of the torus $, the full noisy quantum dynamics T^tv is defined 
(in case of Model I) on A%{6) by the composition 7^ jy := ^eW$. We will also consider 
coarse-grained family of quantum operators defined as follows 

'^}!"J ■■= OM^Ge. (5.28) 

We can now introduce the notion of quantum dissipation time. 

Tg{e, N) := min{n G Z+ : ||T^;^^|| < e"^}, (5.29) 

fcie,N) := min{n G Z+ : ||tj;J|| < e"^} (5.30) 
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Similarly as in classical case the dissipation time provides an intermediate scale be- 
tween initial stage of the evolution (where conservative dynamics dominates due to 
the assumption of negligible contribution from the noise term) and final stage when 
the noise has already driven the system to its final equilibrium (maximally mixed) 
state. On the dissipation time scale the contributions from these competing terms 
are, roughly speaking, balanced. In the following sections we analyse the behavior of 
the dissipation time for fixed quantum system and for its scmiclassical limit. To avoid 
any confusion we will reserve the symbols T^, T^, Tc(e), ^^(e) for classical quantities. 
We note that in case on Model II, due to the above established isometry between 
quantum and classical propagators f Section l5. 2. 1|) and noise operators fSection l5.3p . 
quantum dissipation times for this model coincide with their classical counterparts 
(and do not depend on the value of the Planck constant). Thus in this case no further 
analysis is necessary. On the opposite, as we will see in the following section, in case 
of Model I, quantum dynamics differs considerably from its classical version and in 
order to recover similarities it will be necessary to consider appropriate semiclassical 
limit. 

The main goal of the present section is the semiclassical analysis of the dissipation time 
of noisy evolution of quantum toral maps. Since semiclassical analysis is meaningful 
only in the case of Model I, form now on all the considerations are restricted to this 
case. In order to better understand the need of semiclassical analysis in the case 
of Model I we start with some simple considerations regarding the behavior of the 
dissipation time for finite quantum systems with fixed A^. 

5.4.1 Dissipation time in the 'quantum limit' 

We denote by W<j, the quantum Koopman operator associated with a canonical map 
$ on the torus T^'^. Since dira Aj^{6) < oo there exists a unitary matrix (quantum 
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propagator) f/jy implementing W$ on Hn{G) (cf. Proposition 15 .111) . That is 

u^A = adiUN) = u*^aun, a G Aij{e). 

Since Un is unitary, there exists a basis of its eigenf unctions ipl G 7iAr(0). For each 
such function, V(q,\ipl. ){i'\. \ = \i'\. ){wk I- Moreover W$l = 1, and hence we have 

Proposition 5.13 The degeneracy of unity in the spectrum of any quantum Koop- 
man operator on Aj^{6) is at least of order N'^. For any fixed value of N the corre- 
sponding quantum system is nonergodic. 

In Section 12.41 we showed that the classical dissipation time behaves in a power-law 
fashion in e if the Koopman operator has a nontrivial eigenfunction (if it possesses 
a modicum of regularity). In finite dimensional quantum setting all observables are 
'smooth', since they are represented by a finite Fourier series. Thus one should expect 
that the existence of nontrivial pure point spectrum of the quantum propagator should 
lead to slow dissipation. The following Proposition formalizes this intuition. 

Proposition 5.14 Assume that the noise generating density g decays sufficiently fast 

at infinity: 37 > 2d s.t. \g{x)\ = 0{\x\~"') as \x\ — > 00. 

Then, for any N > and any 0, the quantum noise operator on An{0) satisfies 

||(l-a,^)||<C(eiV)^. (5.31) 

This hound is useful in the limit tN -^ 0. 

As a result, the quantum dissipation time associated with any quantized map U<s, is 
bounded from below as Tq{e,N) > C{eN)~^ , where C > is independent of the 
classical map $. Besides, in this regime the coarse-grained quantum dynamics does 
not undergo dissipation, meaning that fq{e,N) = 00. 
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Proof. We use the RHS of the exphcit expressions for the eigenvalues 7e,iv(fe) of Ge,N- 
From the decay assumptions on g, we see that in the hmit eA^ — > 0, 

The sum on the RHS converges because 7 > 2d. Therefore, we get •y^^N^k) = 1 + 
0(^{eNy^ uniformly w.r.to k G Z^. Since Qe,N is Hermitian, this yields the estimate 

(jESD). 

The lower bound for the quantum dissipation time then follows from the same con- 
siderations as in the proof of Theorem 12.121 ■ 

In view of the above Proposition, the only nontrivial information regarding chaoticity 
of quantum systems and asymptotics of the dissipation time can be retrieved in 
appropriate semi classical limit. We thus turn now to semiclassical analysis of quantum 
maps. 

Following the notation introduced in Section 15.1.31 we denote by IIi^ an orthogonal 
Galerkin-type projection of Lq(T^'*) onto its subspace Jn- It is easy to see that the 
map 

defines an isometric embedding of a finite dimensional quantum algebra of superop- 

erators B{A%{0)) into infinite dimensional classical one B{Ll(T^'^)). 

It has been shown in |1UUJ (see Lemma 1 and its proof there) that for any fixed 

e > 0, the operators (Tn(71^n), which are isometric to T^^n, converge uniformly (i.e. 

in the norm of B{Lq{T^'^))) as A^ — > 00 to the classical operator T^. This implies in 

particular that for any fixed e > 0, and n eN the sequence aNiX^N) converges to T" 

in B{Ll{T^'^)) as iV ^ 00. 

The above result, valid for arbitrary quantizable canonical toral map and any compact 

noise operator, implies that in appropriate semiclassical regime N —>■ 00 one recovers 

classical behavior of the dissipation time. More precisely we have 
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Proposition 5.15 For any quantizable, canonical map F on the torus and any noise 
generating function g, quantum dissipation time coincides asymptotically in suffi- 
ciently fast classical limit N = N{e) ^ oo as e ^ with its classical counterpart. 

Proof. We have \\TJ^j^\\ = ||a7v(7;»|| < lkjv(^e» - T^W + \\T^l Now for given 
e we set n = Tc{e) and choose A^ = N{e) such that \\aiy(T"j^) — T"|| < e. Thus 
ll^'Vll < e + e^^- On the other hand 

\\T, II < IIottvIT;^^) -TJI + ||a7v(T,_^)|| = \\aN[T^^N> - T, \\ + \\T,^n\\ 

Thus for n = Tc{e) — 1 and we have e~^ < e+ ||7^"yY||, which together with the previous 
inequahty establishes the desired result. ■ 

The above statement, despite its generality gives, no information about the actual 
behavior of Tq unless the behavior of the classical one is known. 
The behavior of classical dissipation time of general nonlinear maps has been analyzed 
in Part I of this work. In particular we can use here our result regarding the loga- 
rithmic asymptotics of the dissipation time of Anosov systems (see Theorem I3.24|) . 
Joining this result with the above Proposition we get 

Corollary 5.16 Let F be a volume preserving C^ Anosov diffeomorphisms on T'^ and 
let g he a C^ noise generating function with fast decay at infinity. Then there exist 
Ai, y42 > and C > such that 

I) Quantum dissipation time of the coarse-grained dynamics satisfies in sufficiently 
fast semiclassical and small noise limit e{N)N -^ oo the following estimate 

A,\n{e~\N)) -C< f,(e,iV) < A2He~\N)) + C, 

II) If in addition F has C^"^'"^ -regular foliations and g G C'^{M.'^) is compactly sup- 
ported, then in the above specified semiclassical limit the dissipation time of the noisy 
evolution satisfies for some C > 

A, He-\N)) -C< r,(e, N) < A^ He'^) + C 
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One has to note here however that the crucial from the semiclassical point of view 
question about the regime i.e. the relation between A^ and e for which the above 
asymptotics holds remains in this general setting open. In the next section we answer 
this question in full details in the case of linear Anosov systems. 

5.4.2 Classical limit for quantum toral symplectomorphisms 

In this section we analyze linear maps projected on the torus (generalized cat maps). 

In this case all computations can be carried out explicitly. 

To focus attention and avoid unnecessary notational and computational complications 

we restrict the considerations of this subsection to Gaussian noises (the assumption 

follows the made in classical case, where only a-stable laws were considered). Thus 

we set g(k) = e^''"' . 

It is easy to see that, as in classical case, for any A G A%{6) and any symplectic toral 

automorphism 

n 

This yields 

n n 

lir.'VIH max rr7,^(F-'fc)= max Jl-f^^iF^), (5.32) 

In the above computations we have used the fact that the dissipation time does not 
depend on the direction of time. Similarly in coarse grained case we get 

||tj;;)|| = max 7,^(fc)7,^(F"fc). (5.33) 

We are now in a position to state the main theorem of the present section. 

Theorem 5.17 For arbitrary symplectic, ergodic and diagonalizable F G SL{2d,7j) 
and Gaussian noise, one has the following estimates 
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I) 



r,(e, N) > r,(e), r,(e, N) > ^j^, r,(e, N) > f,(e), 



uniformly in N. 

II) There exists M > (of. \5.40{ l such that 

T^{e,N) ^ r,(e) ^ -^ln(e-i), e ^ , eiV > M, 

Moreover for any (3 satisfying 

d>'-.. 

h{F) 

one has 



InllFll 
(3 > ^^^^^ + 1 (5.34) 



f,ie,N) ^ f,(e)^^ln(e-i), 6^0,6^iV>l, 

h{F) 

where h{F) is a constant equal to minimal dimensionally averaged K-S entropy of F 
(cf ;52]). 

Ill) In the limit eN — > 0, the quantum coarse-grained dynamics does not undergo 
dissipation i.e. Tq{e,N) = oo. The noisy dynamics undergoes slow (power-law) dissi- 
pation. 

As a direct corollary of the above theorem we get the following relation between 
spatial (small e and small h) and time (dissipation and Ehrenfest) scales for classical- 
quantum correspondence of noisy quantum dynamics in case of toral automorphisms. 

Corollary 5.18 (Dissipation vs. Ehrenfest times) 

Under the assumptions of Theorem \5.17\ the following relations hold 

I) 

th'^ ^ oo ^ rg(e, A^) < te, 
II) There exists M > (see 5.40^ such that 



eh-^ r^O M ^ r,(e, A^) ~ te, 
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in) There exists Mq > d^^^'^ such that 

eh-^ ~ C < Mo, r,(e, A^) = oo, 

IV) //A^e(ln(e^^))^'^ <^ 1 and the classical dynamics has logarithmic dissipation time 
(i.e. F is ergodic) then 

Tg{e,N) >rc(e) > r^. 

Proof of Corollary. Statements I) and II) follow immediately from statement II) 
of the above theorem (and its proof). Statement III) is a direct consequence of the 
estimate 



7-W|| > p-d{eNf > 

■i ., AT W -^ C -^ 



-1 



-e,N 



which holds true (for all n) whenever eN < Mq := d ^/^. 

The last statement follows from the following estimates. First we note that 

^.(^'^)>^(]^»l^(^"')~^c(e) 

But we also have 

r,(e) ~ ln(e'^) > ln(e"^) - -lnln(e"^) > InA^ ~ te- 

This completes the proof of the corollary ■ 

To prove the theorem we will need the following estimate (for the proof see Appendix 

Lemma 5.19 

For any $, G M^'^ denote by $, the unique vector satisfying ^- G (—1/2, 1/2] for all 
j = 1, ...,2d and ^ = $ mod 1. There exists a constant C > such that for all a > 
and all | G R'^'^, 

ga{^) < gAl) < l.im <^ + Ce-l'^" < Ul) + Ce-i-' (5.35) 
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Proof of Theorem. 

We start with the proofs of statements I) and III). The proof that quantum dissipation 
time is never shorter than classical one follows from (J5.32J) (J5.33J) . Lemma 15.191 and 
the fact that ge{^) < 1- Indeed, we have 

n n 

||T,>||= max n%M^'fc)> sup n^^(^'^) = ll^"ll- 

Similar estimation holds in coarse-grained case. To show that the remaining assertion 
in I) and that III) hold we first consider coarse-grained version. Denoting $,q := N^-^k 
and ^„ := N^^F^k, we have for any k E Z^, 

Ilt!;;^|| > leNikhMF'^k) > gMio)9eN{L) 
= e-(^^)'(ICoP+l4P) > e-^(^^)'. 

And obviously e"'^^"^^ — > 1 as eA^ -^ 0. 

In noisy case, continuing the notation ^^ := N^^F^k^ we have 



\%y\\ > \{l.N{F^k) > UgMl) = e-(^^)^S-l^''l^ > e- 



1=1 1=1 

Thus 

2 1 



r,(e,iV)> 



Now we pass to the proof of the statement II). 

The lower bounds for both noisy and coarse-grained versions follow from the general 

estimate established in point I) and results obtained in classical setting. 

We turn now to upper bound computations. First we consider coarse-grained version. 

In view of fj5.33|) we have to estimate from above the following product 

\\f}''J\\= max 7,^(fc)7,^(F"fc). (5.36) 

Given j3 satisfying ()5.34|) we fix 5/3 > 0. According to the assumption that e^A^ > 1, 
for all < 5 < (5/3, all sufficiently small e > and all sufficiently big A^ there exists 
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ra G N such that 

1 



- ln(e-M <n< . ln(iV/2). (5.37) 

{l-5)h{F) {1- 6) h{F) In \\F\\ 

Now, given fco G Z^ we consider two cases 

a) feo generates classical orbit i.e. feo, -F"fco € Z^. Then in view of Theorem 13. 10^ for 
any 6' < 6, for sufficiently small e and any n satisfying (j5.37|) we have 

lk,\^ + |F-fco|' > min (|fc|2 + iF^fcH > e^d-'^')^^^)" > e^^^'^^^^^)" (5.38) 

Thus for any such k^, there exists / G {0, ra} such that 

Using ()5.35p we arrive at the following upper bound 

^NikhMF^'k,) < ^MF%) < e~^'^'^'-'^'^''"+Ce-'^^^^^' < e~-\ 

where the last inequality holds for sufficiently big eA^. 

b) fco generates non-classical orbit i.e. F^fco ^ Z^. In this case, |A;o| > e*^"'^^'^^^*^'^^". 
Indeed otherwise, in view of RHS of ()5.37p we would have 

\F"-hJ < II F"||p(l-^)'*(^)" < 

P "'01 — 11-' II t' ^ „ ; 

which would imply F"'k G Z^. 

Thus similarly as in previous case we have for sufficiently big eN, 

The above two cases exhaust all possible values of ko- We then conclude that 

which in view of the definition of dissipation time completes the proof in coarse- 
grained case. 
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Now we consider fully noisy case i.e. we need to estimate from above the following 
product 

n 

ll^!^^ll=o--.n7-(^'fc)- (5.39) 

Let C denote the constant of the RHS of ()5.H5|1 and define 



M := max|4e^(^)||F||, j'41n (^TT^^t;^) 1 ■ (5-40) 

We fix < 5 < 1/2. Using the fact that eA^ > M we obtain for every e > the 
existence of n G N such that 

■ln(e"^) <n < n + 1 < In ( 7777^777 ) • (5.41) 



il-6)h{F) {l-6)h{F) \2\\F\ 

Now, once again, given k^ G Z^ we consider two cases. 

a) ko generates classical orbit i.e. -F'fco G Z^ for / = 1, ..., n. Then in view of Theorem 
I3.1()| for any < 6' < 6, all sufficiently small e, and corresponding n satisfying 1)5.411) 
we have 



V|F^fcol'> min y|F^fc|2>e2a-'5'Wi^)(n+l)>g2(l-5)h(F)(n+l)_ r^,^^) 

1=1 ^ 1=1 

Thus for any such ko, there exists Iq G {1, ...,n} such that 

b) fco generates non-classical orbit i.e. there exists 1 < I < n such that -F'fco ^ Z^. 
Let lo be the largest exponent such that F^°ko e Z^ but F'o+^fco ^ '^n- Then 

|F'°fco|>e(^-')''(^)("+^). (5.43) 

Indeed, otherwise in view of the RHS of (|5.4H) one would have 

|F'»+ifco| < ||F|||F^ofco| < ||F||e(i-^)'^(^)("+i) < -, 
which would imply F'^+^fco G Zj^. 
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Thus in both cases, using (J5.35J) and (j5.4(jp . we get 

n 

1=1 

_ph{F) ^ _m£ _i 

< e "" +Ce 4 < e ^ 
This completes the proof of statement II) and the whole theorem. ■ 

5.5 Technical proofs 

Proof of Proposition 15.101 

In the proof we follow the approach presented in |lU8j . 

Since we already know that symplecticity is necessary for quantization it is enough 
to prove that for given h = N~'^ and 6 the condition ()5.20|) is satisfied iff Up is a 
*-automorphism of the algebra Aj^lO) which in view of ()5.19|) and assumed symplec- 
ticity is equivalent to 

N N 

—F'm V F'm + F'm A6 = —m y m + m AO mod 1. (5.44) 

4 4 ^ ^ 

If A^ is even then for all m, ^m \J m ^TL and hence the above condition is trivially 
satisfied for = 0. If A^ is odd we first note that due to the following identities 

F'mNF'm = {J_F^ f_)mJ+J_F^ Am 

= mJ^FJ^J+J^FU-m = mJ^FJ+FU-m 
F'm A6 = mA F0, 

condition ()5.44|) can be rewritten as 

— m\/m-\ mJ^FJ+F' J^m + m A FO = m A modi (5.45) 

4 4 

To see that ()5.44|) implies ()5.20p we take m = Bj {j E {1, ...,2c/}), where Cj denote 

the standard basis vectors, use the fact that Bj V e^ = and get 

A^ 

—BjJ^FJ+F^J^ej + Bj AF9 = Bj AO mod 1. (5.46) 
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Now we note that 



A^ 



^ e, J_FJ+Ft J_e, 



N 
4 -^ 



AD^-BC^ -AB^ + BA^ 

N IC-D 
A-B 



and 



N C-D\ N A-B 

Substituting the above expression in ()5.46|) immediately yields fj5.20p . 

Now we sketch the proof of the opposite implication. First we rearrange the LHS of 

^ 2d ^ 2d 2d 

— y^ rriimjeiJ^ej -\ 2_. fnirrijeiJ-F J^FU-Sj + y^rriiei A F6 



N 



2d 



i,j=i 



i=l 



— y, 'rnimjei{J_F J^F"^ J_ + J+)ej 



N 



2d 



2d 



+ — ^ m'^^eiJ^FJ+F^J^ei + J^ rriiei A FO 
Now we note that due to the symplecticity of F 

^ 2d ^ 2d 

— 2^ f^i'iTT'j^iJ-FJ^F^J^ej = — 2, minijeiJ^ej mod 1 



4 



4 



Thus it is enough to use the following obvious identity 

2d 2d 

y^mfeiJ_F Jj^F"^ J_ei = y^rriieiJ^FJj^F^ J_ei mod 2, 



4 = 1 



j=i 



to conclude in view of the first part of the proof the equivalence between fl5.44|) and 
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Proof of Proposition 15. Ill 

Denote by {bj} {j = 1,...,N) the standard basis of M^ and let Ejk := \ej){ek\, 
Fjk := T{Ejk). r is inner if there exists U G M.n such that Fj^ = U*\ej){ek\U = 
\U*ej){U*ek\ = \fj){fk\j where fj = U*ej. It is then enough to show that there 
exists an orthonormal basis {fj} such that Fjk = \fj){fk\- First we note that Fjj = 
r(E,,) = T{E*^) = F*^, F,,- = r(E,,) = T{E]^) = Ff^ and F,^F,, = TiE,,E,,) = 
SjkFjj, which implies that Fjj form a set of N mutually orthogonal rank-1 projections. 
Thus there exists an orthonormal basis {fj} such that Fjj = \fj){fj\- Moreover since 
Fjk = FjjFjkFkk one also finds that Fjk = \fj){fk\- ■• 

Proof of Proposition 15.121 

Since by definition Q^^^ is already given in Kraus form, it is completely positive [221 El- 
Trace preservation follows from the normalization. The value of the normalization 
constant can be found as follows 

Z=J2 UN-'k) = Y, 9eiN-'k) = N^' Yl 9Mk) = N''~gMO). 
Now using periodicity of ad{Wk, we get 



^N 



Qe,N can thus be diagonalized in the basis {Wk} 

ye,NVVko = 77^ > geN[k)e W kg = ~ — Wko, 

where in the last inequality we used Poisson summation formula and symmetricity of 

QeNik). 

Hence for any k G Z^'^, Ge,NWk = 7eAr(fc )W^fc, which yields spectral representation 
of Q.N- ■ 
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Proof of Lemma 15.191 

Let fei/2 = (1/2,0, ...,0) G M^"'. There exists a constant C depending only on d such 
that 

Tiez2d 
Now 

neZ n>0 n>0 nGZ neZ 

Thus 

id) < ^.d) + 2Ce-^ 5^ e-'^'l"l'=^(aO + 2Ci(0)e-^ 



and finally 



^a(O) ^a(O) 



To prove the other estimate we consider a splitting of the lattice Z^'^ into three pairwise 
disjoint sets Z^'^, {0}, Z^"' such that n G Z^'^ iff -n e Z^'^. With this notation we have 

TiGZ2d 



Now using the fact that e"^'^^^'"' + e^" ^^ > 2 we get 

TIGZ2'* 

Thus in particular 7eAr(A^^) > geN{^)- • 
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Appendix A 



The dynamics of cat maps 



For completeness we briefly recall in this appendix the most important facts regarding 
classical cat map dynamics on the plane and its canonical quantization. For simplicity 
we present the material in 2-dimensional setting. 

A.l Classical dynamics of cat maps 

In this section we describe both discrete and continuous-time classical cat map dynam- 
ics on the standard phase-plane M^. In continuous setting we derive the Hamiltonian, 
the Lagrangian, the Euler-Lagrange equations and the action of the cat map dynam- 
ics. This will constitute the basis for canonical quantization presented in the next 
section. 

Discrete dynamics 

The classical mechanics of a cat map with 1 degree of freedom is defined on 2 dimen- 
sional Euclidean phase space M^. Any point in the phase space is denoted by {q,p) 
and represents respectively position and momentum coordinates. The initial state of 
the system is denoted by {qo,Po)- The discrete time dynamics is generated by any 
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ergodic toral automorphism A G SL{2, Z) 



qi 


= A 


% 





a b 




% 


Vi_ 




Po_ 




c d 




Po_ 



In the case of 1 degree of freedom {d = 1), considered here, the ergodicity of A is 
equivalent to its hyperbohcity and is determined by the condition |TrA| > 2. For 
simphcity we will even assume that A is positive definite i.e. TtA> 2. 



Continuous dynamics 

Continuous version of the dynamics is defined by suspension over the discrete time 
sequence 



Pit) 



A' 



Po 



In order to provide an explicit formula for this dynamics we need to compute A* 



^tlnA 



We will use the following lemma 



Lemma A.l Let A G Ain and mA{z) = X]i=i('2 '~ ^j)™' be a polynomial which 
satisfies m^(A) = (itla can be e.g. the minimal or the characteristic polynomial 
of A). There exist matrices Mij G Ain such that for any sufficiently differentiable 
function f 

d mj — 1 
i=l j=o 

We will apply the lemma twice. First we apply it to the cat map A and the function 
In. We denote by A the largest eigenvalue of A and take m^ '■= {x — X){x — \~^). 
Thus for any continuous function f we have f{A) = f{X)M^ + f{\^^)M_. To find 
M± we take f± = x — X"^^ and get 

1 



A-A±i]l = (A^i-A±^)M^i 



M4 



A 



±1 



AT 



jiA-X^'l) 
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hence 






In particular 



HA) = J^^i2A-TrAl) = J^ 



a — d 26 
2c d — a 



(A.i: 



The prefactor can be written in many different ways, e.g., 

InA _ h{A) _ h{A) _ sinh-^(|v/(TrA)2-4) 

A-A-i ~ A-A-i ~ ^(TrA)2-4 ~ v^(TrA)2-4 

where /i(A) denotes the KS entropy of A. 

Now we apply the lemma to the matrix \nA, denoting by ;Ui = InA = h{A) and 
/i2 = — In A = —h{A) its eigenvalues. We have mi^^j^ = {x — h{A)){x + h{A)) and for 
any / and /(InA) = f{h{A))M+ + /(-/i(A))M_. To find M± in this case we take 
first f± = x± h{A) 

lnA± h{A) = ±2h{A)M± ^ M± = ^^ \nA + ^1, 



hence 



f(lnA) = 1 fihiA))-fi-hiA)) fihiA)) + fi-hiA)) 

■'^ ' h{A) 2 2 ■ 



In particular for any even function /, /(InA) = /(/i(A))IL (e.g. (InA)^ = /^^(A)!) 

(fe(A) 



and for any odd function /, /(InA) = L ... InA. Taking /(x) = e^ we finally get 



A* = — — sinh(/i( A)t) In A + cosh(/i(A)t) 1. (A.2) 

From the above formula we see that A* can be considered as a solution of a Cauchy 
problem for the following matrix-valued ODE 

X-/i2(A)X = 0, 

X(0) = 1, X(0)=lnA. 



168 



Appendix A. The dynamics of cat maps 



The fact that A^ indeed satisfies this equation can also be verified directly using the 
identity (InA)^ = h'^{A)l and the obvious fact that A* = (InA)^A*. This equation 
will appear later as the Euler-Lagrange equation for this dynamics. 
Combining ()A.2|) with ()A.1|) we arrive at 



A' 



smh{h{A)t) 
^{TiAy-A 



,-d 2b 
2c d-( 



cosh{h{A)t)l. 



Now we write down the explicit equations of continuous version of cat map dynamics 
smh{h{A)t) I fur A\ \\ /n\ oz, sinh(/i(A)t) 



Pit) 



2c 



v/(TrA)2-4 

smh{h{A)t) 
^(TrA)2-4 



[a-d) + cosh{h{A)t) g(0) + 2b " "^^^^'^^^^' ^p{0) 






v/(TW 



g(0) 



^^L^id-a)+cosmA)t)]piO). 



Hamiltonian and Lagrangian 



To find the Hamiltonian of the system we compute 



m 
Pit) 



In AA* 






In A 



Qit) 
Pit) 



In horizontal vector notation this can be rewritten as iq,p) = lnAiq,p). 

The Hamiltonian Hip,q) satisfies (g,p) = V"'"iJ(g,p) = iHpiq,p),—Hgiq,p)). Thus 

V-^if = lnA and 



H 



KA) 



v/(TrA)2-4 



which gives 



Hiq,p) 



HA) 



v/(TrA)2-4 
We note that the Hamiltonian is time independent. 



— c ^ia — d) 
^ia — d) b 



ibp^ — cq^ + (a — d)pq). 



(A.3) 
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The explicit form of Hamilton equations reads 

We can now find the Lagrangian using the Legendre transform 

L{q,q) =pq-H{q,p) (A.4) 

We note that since H is time independent, L will not depend explicitly on time either 
and we drop the time variable. The use of the first Hamilton equation, yields 



P=Yb[ ^^1(1) ^ 'i-(--d)i\- 



Now inserting this formula for p into ()A.4|) and performing necessary transformations 
and rearrangements one arrives at 

L{q, q) = ^ (V(TrA)2-4/i(A)g2 _ 2(a - d)qq + ^/ {Tr Af - 4{h{A))-' f'^ . 

Euler-Lagrange equation 

Now we find the Euler-Lagrange equation for cat map dynamics. The general form 

of the Euler-Lagrange equations can be written as 

d dL _dL 

dt dq dq 

Using the explicit formula for L one gets 

^ = ^^{y'{1,Ar-i{h(A))2q-2{a-d)q). 

which immediately yields the following Euler-Lagrange equation (Cauchy problem) 
for our system 

q-h\A)q = Q, 

g(0) = go, g(0) = vq, 
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where Vq denotes the initial velocity. 

The action 

Now we want to compute the action along the classical path of our dynamics. In 
terms of Lagrangian the action is defined as 

Siquqo) = Siqih),qito))= r L{q{t)A{t))dt. (A.5) 

J to 

Here we derive only equations for action and compute it for particular choice of time 
interval, i.e., to = and ti = 1, that is, we find the action associated with the discrete 
dynamics. To derive action equations we use once again the Legendre transform ()A.4|) 
in conjunction with Euler-Lagrange equations to get 

dL dL 

oq oq 



Now from ()A.5|) we get 

^g(ti) = ^^(g(ti),g(to)) = L{q{t{),q{t^)), 
^g(to) = —S{q{U),q{U)) = -L{q{to),q{to)). 

Dropping the explicit dependence on time and taking derivative w.r.t. q we obtain 

dS dL{q,,q) dS dL{qo,q) 

oqi oq ago oq 

Hence the action equations are 

dS dS 

Pi = 7^, Po 



dqi dqo 

Using these equations one easily checks that for our discrete dynamics 

S{qi, go) = ^("^0 - 2gogi + dqf). (A.6) 
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Indeed 



q' = T{dqi-qo) = T{d{aqo + bpo)-qo) 

{da — 1) 



b 



go + dpo = ego + dpo = pi, 



dS{qi,qo) I, ^ 1/ r N 

g = T(«go-gi) = Tiago -ago-opo) = -Po- 

A. 2 Canonical quantization of cat maps 

Notation 

Throughout the rest of the appendix the following notation will be used 

Zn = Z mod Ar = {0,l,...,iV-l}, 

1 iV-1 



[Z/N) mod 1 = <^0, — ,... 



N' ' N 
Continuous quantum Fourier transform will be denoted by J-'h, i.e. 

M^)iP) = j^j^me~'^'"^dq. (A.7) 

With this normalization, Dirac delta function (on R"') and periodic Dirac delta comb 
(on T"*) can be written as 

5^{q) = ^ f e^^-J-^^Prfp = / e2"^('?-")Prfp, 



For any periodic function ip with period 1 and any vector k G hZ'^ with h = 1/N and 
A^ G N we denote by ip{k) its quantum Fourier coefficients defined as 

V^(^) = 7^/^(9)^"'^^'^ ^9. (A.8) 

Then the inverse transform (i.e. the Fourier representation of ip) is given by 

i,{q) = h'"^ Y, V^(fc)e2"*^. (A.9) 
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The normalization in ()A.8|1 implies that the following relation between continuous 
and periodic Fourier Transforms holds 

Th{^){p) = h'' J2 ^ik)SM (A.IO) 

Parseval identity in this setting reads 

•^ ^ k€hZd 

For a sequence c = {csj^gQd we define its discrete Fourier transform c = {ck}keQ'^ ^^ 

(A.ll) 



Then 





± 




CsC 


-2mNs-k 


L-fc 


J\fd/2 






1 


1: 


4e' 


2niNk-s 


•-S 


p^d/2 





(A.12) 



In discrete case the Parseval identity takes the standard form 



{c,d) = ^ Csds = ^ Cfcdfc = {c,d). 

The relation between continuous and discrete Fourier transforms for periodic delta 
combs and h = 1/N is given by 



sez^/N J keZ'i/N 



Finally we note that 



J2 cJs = N'"^ Y. Cfce^^''^'^-" (A. 13) 

Quasi-periodic wave functions 

Now we are in a position to determine the space of all quasiperiodic wave functions. 
Let 6 = {6q, Op) G T^"^. The wave function ip is 0-quasiperiodic (cf. 113121]) if for all 
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m = {rrii, 7TI2) G Z^'^ 

The second condition gives 

Using the bijective property of the Fourier transform on iS'(M'^) we get 

J2-Ki{eq-q/h)-m2 — 1 

which gives 6q — q/h G Z'^. Hence the only possible values of q for which '?/' 7^ are 
determined by condition q G hU^ + hOq. The wave ip is then necessarily a delta comb 
of the form 

^(9) = X] as6s{q). 
Now 

ip{q + mi)= ^ as5s(q + mi)= ^ as+rrii5s(g) 

sehZ't+hOg sehZ-i-mi+hOq 

Thus quasiperiodicity of ■?/' implies that hU^ — mi = hU^ and ag+mi = e^'^*^^"*^as- 
The former condition implies the existence oi N E Z such that h = 1/N and hence 

^(9) = JjdJ^ Yl ^s^s+Oq/Niq), (A.14) 

where Cs = N'^ag+Oq/N is a quasi-periodic sequence of arbitrary numbers supported on 
U^ /N lattice satisfying Cs+mi = e^'^'^^^''^^ Cs- Throughout the appendix the following 
normalization will be used 

W2 = ^EN' = 1- (A-15) 

This enables us to identify the space of all admissible quasiperiodic wave functions 
with the Hilbert space C^. 
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Quantum propagator 

In this section, following jHSl and (HHl; we recall the explicit formula for the planar 
quantum cat map propagator. We denote by f/(9i, Qq) the kernel of the propagator 
i.e. the Green function for the Schrodinger equation for quantum cat map dynamics 
on M^. In order to find U{qi, Qq) we can either quantize the Hamiltonian (derived in 
Section IA.1|) and solve the corresponding Schrodinger equation, or use the fact that 
the Hamiltonian is quadratic which implies that the semiclassical approximation to 
the quantum propagator is exact [HH] and hence the kernel can be directly expressed 
in terms of the classical action given by formula ()A.6|) in section lA.ll Applying the 
second method we immediately get 

Uiq„qo) = (^^) exp(^l(aq^-2qog, + rfq?)) 

The quantum evolution ipi = Uipo of an initial wave function ipQ is then given by 

MQi) = U{q^,qo)MQo)dqo- 
In particular, for an initial wave in form of a 0-quasiperiodic delta comb 

V'o(q) = Yl ^^^s+0,/N{q), 

the evolution is given by 

MQi) = lu{qi,qo) Yl cA+0,/iv(go)c^9o (A.16) 

^^"^ sezd/N 

= (t) E c.exp(^iais + 0,/Nf-2is + d,/N)q, + dql)Y 

We will use the above formula in Appendix |Bl where we show that the evolution of the 
Wigner transform of ipo is given by the Frobenius- Perron operator associated with a 
classical cat map. 
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Appendix B 



Wigner transform 



This appendix is devoted to the Wigner transform. By means of distribution the- 
ory we derive an exphcit form of its discrete version for wave functions assuming 
quasiperiodic delta comb form. We show that the discrete Wigner function takes the 
form of a uniformly spaced Dirac "delta brush" , which is periodic and supported on 
the 2N X 2N rational grid on the unit torus. We derive several important proper- 
ties of the discrete Wigner transform and show that it provides a clear geometric 
interpretation of the quantum cat maps evolution f Proposition IB.1|) . We conclude 
the appendix with a simple proof of the equivalence between canonical and algebraic 
finite dimensional quantizations of the cat map dynamics. 

We start by recalling briefly standard definition of the Wigner transform. In order 
to simplify the notation we will use the symbol x = {q,p) to denote the phase 
space variables and the symbol k = (fci, ^2) to denote the conjugate variables in the 
frequency space. 
Let ip denote a wave function. The associated Wigner transform W^{q,p) is given by 



W^x) = W4q,p) = ^^J^y(^q + ^j^(^q-^]e-'^^-'^'dq' 
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where p{u,v) = \ip){ip\ = ip ® ip = ipiu)ip{v) is the density kernel associated with 
the wave function ip. 

In its most general form the Wigner transform is defined for any tempered distribution 
ip G 5'(]R") or density p G iS'(]R^"). Using this fact one can derive a useful integral 
version of the Wigner transform. To this end, let denote an arbitrary test function 
from tS(M'^) and apply the change of variables 

q' q' 

u = q + -, v = q--; 

to get 

/ Wp{x)<P{x)dx = ^/ p(q + ^,q-pje-'^^-'^'<P{x)dq'dx 

= -^ / p{u,v)5J^^^^]e~'^P<'''''^(P{x)dudvdx. 
h J^4d \ 2 / 

From above computations we see that the kernel of the Wigner transform is given by 
the continuous family of integral Fano operators Sx 

Thus the Wigner transform Wp{x) of p admits the following kernel representation 

1 If 

^pi^) = Td^^P^'^^Td p{u,v)Sx{v,u)dvdu 

Discrete Wigner function 

In this section we construct a discrete version of a Wigner transform. 

Using the kernel representation and applying it to a quasiperiodic distributional wave 

function ijj of the form ()A.14|) with h = N~'^ we get 

W4x)= V c,,c,,^.,+.. .Jq)e-^"^P-(-^--^). (B.l) 

^-^ 2 "'' JV 

si,S2€Z''/N 
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Our main goal is to show that ()B.1|) assumes a periodic "delta brush" structure. To 
this end we first change summation indexes in ()B.H) . 

Si + S2 Sl ~" •S2 

r = s = 

2 2 

Note that the above transformation does not merely undo the change of variables 
performed in deriving the kernel representation of W^ (the Jacobian this time is 
non-unital). 

r,seZ'i/2N 

(We use a convention that Cq = whenever q ^ Z/N.) Next we have 
= E c.+*c._,e-^-^^-* E ^^(2^(P - |)) 

Thus 

^^(^) = 7^ E E -.+*c._,e---^(^-^^)X^«,(g)5^^^(p) 

The above formulas provide a clear geometric interpretation of the Wigner transform 
of a quasiperiodic delta comb. A few remarks are in order here. First of all, the 
resulting Wigner function is strictly periodic (even though the original wave need not 
be) and is supported on the grid with a mesh spacing of the size of h/2 (half of the 
corresponding resolution for the wave). The support of the Wigner transform forms a 
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2N X 2N lattice centered at (shifted from the origin by) 6. The later property explains 
the interpretation of 6 as Bloch or Floquet "angles" (cf . Section 15.1 .41) and the term 
"quantum boundary conditions" coined in [TU] in the context of the quantization 
condition 15.201 (see Proposition 15.101 in Section IB. 2. Ij) . 

Now we want to find discrete Fourier coefficients of W^. Taking the Fourier transform 
in this case requires some care, since one has to adjust the value of the Planck constant 
adequately. Taking into account the support of W^, its discrete Fourier transform 
agrees with a distributional one if the latter is taken with the value of the Planck 
constant equal to h/2. Indeed, for any k G 1?'^ /2N , we have 

The above formula can be simplified as follows 
1 

W) 



W^{k) = -^ Yl c.+tC._,e-2-2^'=-(^+^) Y e-^-^^(^+^)-(*+^^) 



i2N)~ 



2jv ■'^V2JV 

— y 



c.-.,c.+.,e-2'^™-(^+i^). 



For further simplification one needs to note that half of the coefficients are zero, since 
Cr are supported on Z/A^. Moreover due to the quasiperiodicity of c^ the product 
Cj-Cr is periodic. We can thus apply the change of indices r = £ — ^2 to obtain 



^V-l^) = -n^d X] ^rCr+2k^er 



-2Tri2Nki{r+k2+j^) 



(2N¥ ^ CrCr+2k2(i 

_ __}^^-2m2Nki{k2 + ^) \^ „ y, ^^2TTiNki-r 

~ (2N¥ ^ <-r<-'r+2k2<i 

Thus we conclude that the discrete Fourier coefficients of W^ are given by the so-called 
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discrete ambiguity function. 

Wm = (^ E Cr.-fc.c.+..e-^-^^'=-(^+^) (B.3) 



rei. 



(2iV) 



^ $: c.c.+2..e-^-^^'=^(^+^^+^). (B.4) 



red, 



From the above formula we see that W{k) is 1-quasiperiodic in both variables and is 
supported on 1/ /2N lattice. 

Properties of the Discrete Wigner function 

It is well known that in the continuous setting the Wigner function can be viewed as a 
quantum counterpart of the joint phase space density for the position and momentum 
variables. For any ip G LF'iW^) satisfying normalization condition \\il)\\2 = 1 one has 

f W^{x)dp=\ij{q)\\ f W^{x)dq=\J'hij{p)\\ [ W4x)dx = 1, 

JR'i J^d jT^2d 

The Wigner function is always real but need not be nonnegative (one can easily note 
that W^{0) = — (I) ||V'(9)ll2 is negative for any odd tp). However, the following 
useful property (Parseval identity) holds 

Not all of these properties can be generalized in any obvious way to the whole iS'(]R) 
since the waves from outside of L^(]R) are no longer normalizable. Nevertheless, as 
we show in this section, all these properties are preserved for quasiperiodic waves. 
There are also some interesting characteristic properties of discrete Wigner function 
(not shared by the standard continuous version). Below we mention two of them 
I. W{x) is completely determined by {2N)'^'^ matrix of the strengths of its delta 
functions. 

(2NY ^ Cr+t^^r-td 
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II. In general all {2N)^'^ coefficients may be nonzero, but at the same time the coef- 
ficients attached to a vertices of a |-box may differ at most in sign. 






(2N¥ "^ r+t r t 



teQ 



d 



and 



— y 



r,s+^ (2iV)' 









g-.i2Art.e,^^^^-^_^g-2..2iV(s+^)t 



1 J2 (-l)^^*--c.^,c._,e---^(^-^^)*. 



^2]V 



Now it suffices to note that Cr+t are zero whenever r ■ Cj and t ■ Cj are of different 
parity which gives 






>,s+i 



Next we derive the properties regarding marginal distributions 
III. g-marginal projection. 



1 

Jmy 



seQjJV ''^V2JV -^V2JV 

- ^ ^ C ^,C , Y" ^-2.i2N(s+'j)t ^ - _ I ,2 



IV. p-marginal projection results in |csp (we omit the proof, which is in this case a 
bit more technically involved). 

We end this section by deriving the discrete version of the Parseval identity for the 
Wigner function. That is, we want to prove the following identity 
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V. 



j\rd Z^ '^'^""^ 



J2 W^AkmAk) = N' 

In the prove we use the Fourier transform of W^p derived in the previous section. 
J2 W^Akm,{k) 

kGQl% 

- ^ \^ \^ r r . ^-27Ti2Nkiir+k2 + ^)-l J , ^27rJ2iVfci-(s+fc2 + ^) 
~ (2NW ^ ^ CrCr+2k2e ^ dsds+2k2e 

- -*- \^ \^ r r , p-^^iNki-r-J j 2niNki-s 
^ l\(2d / J CrCr+k2d'sU's+k2 / j ^ 



k2,r,s<^Q% 



ki€Q% 



-^ ^ CrCr+k2dsds+k2N'^Sir - s) = N" 



JY2d 



fc2,r,sS 



^N 



]\rd 2-^ '^■^"'■' 



Cat map evolution of the Wigner function 

In this section we show the the Wigner function associated with a quasiperiodic delta 
comb evolves classically under the cat map dynamics. Indeed, we have the following 

Proposition B.l Let U denote the quantum propagator associated with a cat map 
F . For any initial -quasiperiodic delta wave ipQ, the Wigner transform of the evolved 
wave ipi = UipQ satisfies the property W^^ = W^^ o F~'^. 



Proof. The prove is obtained by a direct application of formula ()A.16|) (see Section 
^2]) which gives 



W^,{x) 



Ni\ 



'^''^— ^L(«l + ^)^-2(«l + ^)(q+^)+d(q+^ 



^^ E- 



^ y 



(i/2 



— j E^-^" 



^ ^a(s2 + ^)2-2{s2 + ^)(q-4)+4q- 



2; J ^-27TiNpq' ^^/ 



S2 
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From the above we can represent W^p-^ {x) by 



31,32 






n + ^)A , 2^[a'l^+a{s,-s,f-^-(s,-s,)q 



31,32 

Thus using the spectral resolution of the Dirac delta comb we get 

E2TTiN ( „( SX+S2 . Sq\ „\ / „ „ ^ 
2 ^ AT 
Si,S2 

= y^ Cs,CsJ s,+.2 I 9g (dq - bp)e 

S1,S2 



Sl,S2 



V" Cs,Cs2(5£i+£2 , eg ((ig - 6p)e 

■^-^ 2 "'" iV 

Sl,S2 



Sl,S2 

— c- / r 1 \ OiriKJl — r-n^r,n\l a, —an\ tt^ "^^ ^I 

2 ^ AT 
Sl,S2 



Wigner function, Weyl quantization and ambiguity function 

We end this appendix by recalling briefly the relation between the Wigner function, 
the Weyl quantization and the ambiguity function. We have 

Proposition B.2 For any classical observable f G L'^{M.'^'^) and any wave function 
if) G L'^{E?'^), the expectation of the Weyl quantization Op{f) of f satisfies 



(^,0p(/)^)= / f{x)W4x)dx. 

In particular, the expectation of any Weyl translation operator T^ (cf. ()5.2|) . Section 
\5.1.1]) is given by the ambiguity function A^{k) (the inverse Fourier transform of the 
Wigner function) associated with a wave ip 



{^P,Tk^) = A^{k) = [ W4x)e''^'''^^dx. 



Proof. Weyl translations generate the whole algebra of observables. It is thus enough 
to prove the second statement. Using the formula ()5.3|) for the explicit action of Weyl 
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translations on wave functions derived in section 15. 1.11 we get 

(^,TfcV^) = / i^iq)Tkijiq)dq = [ ii{q)e'^'^^<^-^^'^^i^{q - k^)dq 
= f ^{q-ki/2)'^{q + k,/2))e^'''-'^dq = A4k) 
On the other hand 
/ WJx)e^^''^dx = ^ f ^{q + q'/2)^{q-q'/2))e-^'^P-'''-''^-''+'''-P^dq'dx 

= [ ^{q-kj2)i,{q + ki/2))e^''^-'^dq = A^{k), 

JRd 

which completes the proof. ■ 

Remark B.3 The above proposition generalizes to the wave functions from S'{M.'^) 
(in this case the observables need to be taken from a smaller set). In particular, 
the Proposition holds for quasiperiodic delta combs and the observables from L^(T^'^), 
satisfying Y.k 1/(^)1 < ^■ 



The above proposition together with Proposition IB. II provides the proof that the 
algebraic quantization introduced in Section 15.2.11 coincides with a canonical one, 
derived in Section lA.2l of AppendixIXl Indeed, denoting by U the propagator obtained 
by means of the canonical quantization and by U the *-automorphism constructed in 
the algebraic approach, we get (using symplecticity of F) what follows 

= [ W^,{F-'x)e'-''-^dx= [ W^,{x)e'-''^dx = {^^,Tk^^) 

Thus the canonical quantum propagator U implements the algebraic *-automorphism 
U on the space of quasiperiodic delta waves. 
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